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BACKGROUND  AND  SIGNIFICANCE 

Contamination  of  soils,  sediments,  and  groundwater  with  non-aqueous  phase  liquids  (NAPLs)  is 
the  most  widespread  environmental  problem  facing  the  Air  Force.  In  situ  remediation  technologies 
are  the  only  viable  treatments  that  offer  the  potential  of  gaining  control  over  the  huge  costs  and 
extended  time  required  to  clean  up  contaminated  Air  Force  sites.  Cleanup  procedures  are  often 
instituted  and  operated  with  limited  knowledge  of  both  the  physical  and  chemical  processes  which 
may  govern  contaminant  vapor  and  liquid  mass  transfer  and  transport.  To  be  effective,  complete 
hy^nantic  control  of  the  subsurface  at  the  remediation  site  is  critical.  One  of  the  complicating 
factors  in  hydronamic  control  is  the  creation  of  preferential  flow  paths  generated  by  chaotic 
disturbances  in  the  flow  field.  These  phenomena  are  not  predicted  well  with  current  multiphase 
flow  process  models,  especially  when  gravity  plays  a  role  (Rimmer  et  al.,  1996),  which  dictates 
the  use  of  very  conservative  and  expensive  risk-analysis  criteria.  On  many  occasions,  this  has 
rendered  the  use  of  in  situ  treatment  remedies  ineffective  because  of  hydronamic  bypassing  of 
subsurface  pollutants. 

The  complexity  of  preferential  flow  of  water  and  non-aqueous  phase  liquids  is  largely  unappre¬ 
ciated  because  few  techniques  permit  accurate  measurement  of  water  and  NAPL  contents  in 
rapidly  changing  flow  fields.  Until  now,  rapid  measurements  of  fluid  concentrations  have  been 
limited  to  two-phase  systems  (i.e.,  systems  with  two  mobile  fluids),  and  the  first  part  of  this  work 
focussed  on  these  NAPL  and  water  systems.  However,  in  most  contamination  sites,  NAPLs, 
water,  and  air  coexist  within  the  soil  matrix  yielding  a  more  complex  three-phase  system.  For 
these  systems,  it  is  much  more  difficult  to  obtain  accurate  saturations,  as  two  independent 
measurements  are  needed  at  each  point.  Previous  measurement  techniques  have  been  limited  to 
static  or  steady  state  flow  conditions  as  the  measurement  times  are  long.  This  has  severely  limited 
the  study  of  instabilities  in  three-phase  systems. 

We  have  used  this  Air  Force  grant  to  improve  and  extend  into  three-phase  systems  two  separate 
visualization  techniques  which  can  rapidly  record  fluid  concentrations  in  our  soil  slabs:  light 
transmission  and  attenuation  of  synchrotron  radiation.  The  techniques  are  complementary  and 
both  provide  high  temporal  and  spatial  resolution  of  fluid  concentrations,  necessary  as  the 
properties  of  preferential  flow  can  ch^ge  on  short  time  scales.  Using  these  improved  techniques 
we  have  obtained  data  concerning  the  source  and  scope  of  preferential  flow  in  two-  and  three- 
phase  systems. 

The  specific  research  performed  during  this  grant  is  summarized  as  follows.  During  the  first  year 
of  the  grant,  we  began  testing  the  theoretical  predictions  for  two-phase  flow  by  using  both  tech¬ 
niques.  In  the  second  year,  we  extended  the  study  of  two-phase  flow  to  include  surfactants  and 
their  effects  on  preferential  flow.  The  light  transmission  technique  has  been  made  quantitative  for 
two-phase  oil  and  water  systems  and  we  have  extended  the  technique  to  three-phase  systems  (oil, 
water,  and  air)  (Damault  et  al.,  in  preparation).  Attenuation  of  synchrotron  radiation  has  been 
used  to  provide  highly  accurate  and  rapid  concentrations  within  preferential  two-phase  flow. 
Importantly,  rapid  measurements  of  three-phase  saturations  have  been  achieved  by  using  dual 
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energies  from  the  synchrotron  radiation.  Two  other  related  studies  have  been  carried  out.  The  first 
studied  oil  and  water  movement  through  wetted  and  non-wetted  tubes  to  understand  the  basics  of 
fluid  entrapment.  The  second  study  consisted  of  the  production  of  preferential  flow  as  a  function 
of  the  solid-liquid  contact  angle  as  in  water  repellant  or  NAPL  coated  soils. 

THEORETICAL  BACKGROUND 

The  theory  leading  to  the  instability  of  the  interface  between  fluids  in  a  porous  medium  has  been 
the  topic  of  continuous  research  in  fields  of  chemical  engineering,  fluid  mechanics,  and  petroleum 
engineering  since  the  discovery  of  interfacial  instability  by  Hill  (1952).  Hill  showed  that  the 
downward  displacement  of  fluid  1  into  fluid  2  is  unstable  if 

^^PrPz)  -  0a''(PrP2)  >  0  (1) 

Here  p,  and  pj  are  the  fluid  densities,  ft,  and  H2  are  the  viscosities,  u  is  the  interfacial  velocity, 
k  is  the  intrinsic  permeability  of  the  media,  0,  is  the  porosity  available  for  fluid  transmission,  and 
g  is  the  gravitational  acceleration.  This  instability  condition  is  satisfied  when  water  infiltrates  into 
a  sandy  soil  saturated  with  a  light  oil,  thus  producing  preferential  flow  paths. 

For  the  water-air  system,  previous  work  by  Parlange  and  Hill  (1976)  calculated  the  finger  size  (d) 
for  water  infiltrating  air  dry  soil  which  subsequent  experiments  found  to  be  valid  (Glass  et  al., 
1989).  Under  this  contract,  we  have  broadened  this  analysis  to  include  water  infiltrating  into  oil 
saturated  soil  with  the  assumption  that  the  oil  pressure  remains  constant  ahead  of  the  finger 
(Chandler  et  al.,  1997,  see  Appendix  A  for  more  detail).  The  resulting  finger  width  is  related  to 
the  air-water  width  by 


*oil 
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(2) 


where  o„a  and  are  the  water-oil  and  water-air  surface  tensions,  respectively,  is  the 
conductivity  of  the  tip  in  the  water  finger,  v  is  the  velocity  of  the  front  if  it  would  have  remained 
flat,  and  and  0,^,  are  the  initial  and  final  water  contents  in  the  finger,  respectively. 


EXPERIMENTS 


The  two  high-speed  data  acquisition  techniques  -  Light  Transmission  and  Synchrotron  X-Ray 
attenuation  -  each  have  their  unique  advantages.  The  light  transmission  technique  involves  the 
video  recording  of  light  and/or  color  transmission  through  a  soil  slab.  The  concentration  of  each 
fluid  can  be  directly  related  to  the  transmitted  intoisity  for  the  air-water  system  and  hue  (or  color) 
for  the  NAPL- water  system  in  which  the  water  is  colored  blue  (Chandler  et  al.,  1997;  Rimmer 
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et  al.,  1997,  see  Appendix  B).  Light  transmission  allows  a  spatial  field  of  250,000  points  to  be 
measured  and,  thus,  provides  measurements  of  the  flow  throughout  the  whole  chamber. 
Synchrotron  x-ray  attenuation  is  based  upon  the  same  technique  as  gamma-ray  attenuation,  the 
more  fluid  in  the  radiation  path,  the  more  the  radiation  is  absorbed  (Liu  et  al.,  1993).  The 
advantage  of  a  synchrotron  source  is  that  the  large  quantity  of  radiation  (10,(K)0  times  greater  than 
conventional  sources)  allows  data  to  be  taken  at  a  much  faster  rate  than  from  gamma-ray  sources. 
X-ray  attenuation  measures  fluid  concentration  at  one  position  at  any  time,  but  at  a  much  higher 
accuracy  and  with  much  greater  soil  choices  than  the  light  transmission  technique. 

In  all  of  the  experiments,  a  two-dimensional  polycarbonate  chamber  (thickness  =  1  cm)  was  filled 
with  sand  and  shaken  to  prevent  settling  during  the  fluid  injection.  Several  ports  were  strategically 
located  throughout  the  slab,  in  which  high-speed  tensiometers  were  placed  to  measure  the 
pressures  of  the  oil  and  water.  The  tensiometer  readings  were  fed  into  a  computer  for  direct 
comparison  with  fluid  concentrations.  For  two-phase  measurements,  soils  were  saturated  with 
Soltrol  220,  a  light  non-toxic  NAPL  somewhat  similar  to  kerosene  and/or  quadricyclane.  Water 
was  applied  to  the  top  surface,  to  simulate  rain  or  water  forcing  on  a  contaminated  soil.  Fingering 
was  found  to  occur  and  was  studied  as  a  function  of  soil  sieve  sizes,  water  application  rates,  and 
surfactant  concentration. 

Light  Transmission 

As  both  phases  in  the  oil-water  system  are  clear  liquids,  in  our  original  experiment  the  visual 
contrast  was  enhanced  by  dying  the  water  with  0. 1  mg/1  FD&C  blue  dye  and  the  oil  with  0.2  mg/1 
Sudan  IV  red  dye  (Aldrich  Chemical  Co.).  Rainfall  was  simulated  by  dripping  water  uniformly 
across  the  top,  and  Figure  1  depicts  time  slices  of  the  water  infiltration.  The  measured  width  of 
the  fingers,  5-8  cm,  compare  very  favorably  with  the  predicted  width  from  the  stability  analysis. 
When  surfactants  are  applied  to  the  water,  the  fingers  tend  to  break  up  and  wander  horizontally, 
and  the  water  flow  covers  more  of  the  chamber. 

Recently,  we  have  modified  the  video  imaging  technique  so  as  to  obtain  quantitative  fluid 
saturations.  The  modification  consists  of  dyeing  only  the  water  phase  and  analyzing  the  color  (or 
hue)  content  of  the  transmitted  light.  To  test  and  calibrate  this  technique,  a  chamber  with 
individual  cells  was  built  from  the  same  materials  and  the  same  thickness  as  the  experimental 
chambers.  Figure  2  is  a  picture  of  the  light  transmission  through  such  a  chamber,  with  each  cell 
packed  with  sand  and  different  concentrations  of  oil,  water,  and  air.  Cells  with  more  water 
tqypeared  bluer  to  the  eye,  and  cells  with  more  oil  appeared  brighter  to  the  eye.  By  trial  and  error 
with  different  dyes,  it  was  determined  that  when  a  concentration  of  0.025  mg/1  FD&C  blue  dye 
is  added  to  the  water  phase,  the  measured  hue  content  of  the  transmitted  light  is  proportional  to 
the  water  content  inside  the  porous  media.  This  relationship  is  shown  in  Figure  3.  The  linearity 
of  this  relationship  produces  higher  accuracy  fluid  concentrations  than  those  obtained  in  water-air 
systems  which  use  the  intensity  of  the  transmitted  light.  The  visualization  technique  for  oil  and 
water  has  been  accepted  for  publication  (Damault  et  al.,  1997,  see  Appendix  C). 
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Figure  1:  Water  fingering  in  oil  saturated  soil.  Each  line  represents  the  position  of  the  water 
fingers  at  equal  time  intervals.  The  gray  circles  are  the  positions  of  the  tensiometers  (from 
Chandler  et  al.,  1995,  see  Appendix  A). 

Using  this  light  transmission  apparatus,  fluid  concentrations  within  preferential  flow  paths  were 
measured.  Inside  a  water  finger  penetrating  an  oil  saturated  soil  there  exists  two  distinct  regions. 
The  lower  part  of  the  finger  has  a  high  concentration  of  water  (0.25  cmVcm^),  effectively  dis¬ 
placing  most  of  the  oil.  But,  behind  this  region  the  water  concentration  is  lower  (0. 17  cmVcm^), 
indicating  that  the  oil  is  re-entering  the  finger.  This  behavior  is  similar  to  that  observed  in  water 
and  air  systems.  The  lower  water  content  behind  the  finger  tip  also  inhibits  further  widening  of 
the  finger,  leaving  large  portions  of  the  oil  unmoved.  These  results  match  those  seen  in  the 
synchrotron  work  of  the  first  year,  but  the  light  transmission  technique  allows  for  a  study  of  the 
whole  flow  field,  at  least  for  homogenous  soils. 


Figure  2:  A  picture  of  the  calibration  chamber,  with  cells  filled  with  varying  amounts  of  oil, 
water,  and  air.  The  transmitted  intensity  and  hue  can  be  related  to  the  fluid  contents  inside  each 
cell. 


WATER  CONTENT  IN  SOIL-AIR-OIL-WATER  SYSTEM 
Plot  of  average  Hue  versus  water  content 


Figure  3:  Relationship  between  the  hue  of  the  transmitted  light  and  water  content  of  the  sand. 
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LIQUID  FLUIDS  CONTENTS  IN  SOIL-AIR-OIL-WATER  SYSTEM 
Plot  of  average  Intensity  versus  liquid  fluids  contents 
with  different  water  contents 


20  40  60  80  100  120  140  160  180  200  220  240 

Intensity 


Figure  4:  Relationship  betwewi  the  intensity  of  the  transmitted  light  and  the  total  fluid  content  of 
the  sand.  Different  calibration  curves  are  used  depending  on  the  hue  of  the  transmitted  light. 
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Recently,  we  extended  light  transmission  to  rapidly  measure  saturations  in  three-phase  flow.  When 
the  output  of  the  color  camera  is  converted  into  hue,  saturation,  and  intensity,  we  find  that,  as 
before,  the  hue  is  related  only  to  the  water  content,  and  for  known  water  content  the  intensity  is 
related  to  the  total  fluid  content.  Figure  4  shows  how  the  intensity  depends  on  the  total  fluid 
content  for  different  water  contents.  The  curves  are  identical  except  at  low  water  contents  where 
the  intensity  is  much  greater  for  a  given  fluid  content.  Still,  there  exists  a  one-to-one  correspon¬ 
dence  between  fluid  content  and  intensity  once  the  water  content  is  known  from  the  measured  hue. 
Thus,  we  can  convert  the  measured  color  and  intensity  at  every  pixel  of  the  camera  into  oil,  water, 
and  air  contents  inside  our  sand  filled  chamber.  This  is  the  first  technique  that  can  rapidly  measure 
three-phase  concentrations  over  a  whole  flow  field.  More  details  are  provided  in  Appendix  D 
(DiCarlo  et  al.,  1996). 

Synchrotron  X-ray  Attenuation 

During  the  first  year  of  funding,  we  used  synchrotron  x-ray  attenuation  to  measure  concentrations 
of  fluids  to  a  volumetric  accuracy  of  0.01  on  a  time  scale  of  one  second.  The  source  of  the  radia¬ 
tion  was  the  F-2  beam  line  at  the  Cornell  High  Energy  Synchrotron  Source  (CHESS).  Further 
details  of  the  apparatus  are  described  in  the  first  progress  report. 

In  each  experiment,  water  was  applied  at  a  single  point  in  the  top  at  varying  fluxes  (0.5  to  70.0 
cmVmin).  A  finger  began  to  form  within  10  cm  of  the  top,  which  methodically  moved  down  the 
chamber.  Measurements  of  the  oil  and  water  pressure  at  points  within  the  finger  were  obtained 
versus  time.  Water  concentrations  wctc  measured  approximately  30  cm  from  the  top  of  the  cham¬ 
ber  as  functions  of  time  and  horizontal  position  within  the  finger.  Figure  5  shows  an  oil-water 
hysteresis  curve  obtained  by  measuring  the  water  content  and  pressure  through  an  infiltration  and 
re-infiltration  cycle.  Hysteresis  curves  are  an  important  parameter  which  controls  the  persistence 
of  fingered  flow,  and  can  be  used  to  predict  the  success  of  different  cleanup  strategies. 

For  important  three-phase  saturations,  we  have  developed  a  technique  where  the  harmonics  of 
synchrotron  radiation  produce  the  necessary  intense  dual  energy  x-ray  beam.  The  water  and  oil 
phases  were  doped  such  that  each  phase  was  sensitive  to  each  energy  photon.  This  beam  passed 
through  our  sample  chamber  and  incident  and  transmitted  x-rays  were  measured  using  scintillation 
detectors  which  provided  the  necessary  energy  resolution.  The  chamber  was  placed  on  a  new 
movable  platform  which  had  a  range  of  50  x  30  cm.  A  schematic  of  the  beam  line  is  shown  in 
Figure  6.  Using  this  technique,  we  were  able  to  obtain  three-phase  saturations  with  an  accuracy 
of  .02  volumetric  over  time  scales  as  small  as  5  seconds. 

Our  initial  three-phase  studies  consisted  of  oil  infiltrating  into  dry  and  partially  water  saturated 
soil.  The  chamber  was  filled  with  water  from  below  to  a  hdght  in  the  chamber  of  40  cm,  and  sub¬ 
sequently  drained  until  the  water  fringe  was  at  a  height  of  12  cm.  Thus,  initially  the  soil  slab  had 
thiw  distinct  regions,  12  cm  of  water  saturated  sand  at  the  bottom,  followed  by  28  cm  of  partially 
water  saturated  sand,  and  finally  15  cm  of  air  dry  sand  on  top.  Oil  was  applied  through  a  point 
source  at  the  top  of  the  chamber  at  a  rate  of  4  ml/min.  A  finger  formed  which  moved  quickly 
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through  the  dry  sand,  and  slowly  through  the  partially  wet  sand,  before  finally  hitting  the  water 
table  and  spreading  laterally.  The  experiment  was  run  with  varying  sieve  size  sands,  and  the  fluid 
contents  were  measured  as  the  finger  passed  through  the  dry  and  partially  water  saturated  sands. 

Soltrol  220  -  Water  Dynamic 
Pressure-  Saturation  Relationship 


0.00  0.05  0.10  0.15  0.20  0.25  0.30  0.35 

9^  Water  Content  (cm^/cm^) 


Figure  5:  Hysteresis  of  the  primary  and  secondary  wetting  and  drying  curves  in  an  oil  and  water 
system  measured  with  synchrotron  x-rays. 
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Figure  6:  A  schematic  of  the  experimental  setup  at  the  F-2  hutch  of  CHESS.  The  monochroma¬ 
tor  transforms  the  white  x-ray  beam  into  a  multiple  energy  beam.  Scintillation  detectors  measure 
the  intensity  of  both  energy  x-rays  by  detecting  air  scattering  from  the  main  beam.  The  sample 
chamber  is  mounted  on  a  movable  platform  allowing  measurements  from  many  different  positions 
within  the  chamber. 

Figure  7  shows  the  measured  fluid  saturations  as  the  oil  finger  penetrated  the  partially  water 
saturated  30/40  sand,  one  of  many  data  sets.  Data  were  accumulated  over  3  second  intervals,  and 
the  error  bars  are  the  expected  variance  due  to  the  count  times  in  the  detector.  As  the  oil 
penetrated  the  sand,  it  drove  out  the  water. 

Other  Experiments 

In  a  separate  experiment  (Rimmer  et  al.,  1996,  see  Appendix  E),  we  compared  how  the  theory 
of  wetting  and  nonwetting  liquid  motion  in  tubes  was  related  to  the  size  of  ganglia  in  a  porous 
medium.  It  appeared  that  the  size  of  the  entrapped  ganglia,  when  gravity  alone  is  the  cause  of  the 
motion,  was  closely  ^elated  to  the  tube  size  in  which  the  wetting-nonwetting  fluid  configuration 
was  stable. 

In  studies  of  preferential  flow,  the  contact  angle  is  important.  A  special  case  is  wetting  and 
nonwetting  soils,  in  which  the  finger  width  depends  highly  on  the  degree  of  water  repellency. 


Volumetric  Content 
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Identical  soils  were  treated  with  OTS  (octadecyltrichlorasilane  Fisher  Scientific),  to  varying 
degrees  of  repellency  and,  thus,  contact  angle.  When  the  soil  had  a  positive  water  entry  pressure, 
fingering  took  place,  bypassing  most  of  the  matrix.  The  results  are  presented  in  Appendix  F 
(Bauters  et  al.,  1997). 
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Figure  7:  Water  and  oil  contents  versus  time  as  the  oil  finger  passes  through  partially  water 
saturated  30/40  sand.  The  oil  content  does  not  reach  saturation  and  water  content  is  driven  lower 
by  the  invading  oil. 
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Initial  experiments  have  shown  that  the  addition  of  surfactants  caused  the  water  density  in  the 
finger  to  be  less  uniform,  and  the  infiltrating  water  covered  much  more  of  the  experimental 
chamber,  albeit  with  smaller  fingers.  We  have  repeated  these  experiments  with  four  different 
surfactants,  using  synchrotron  x-rays.  We  found  that  if  we  used  ]^.  (2),  the  finger  width  was 
slightly  over  predicted,  because  the  theory  assumes  that  the  energy  loss  in  oil  is  negligible 
flowing  around  the  tip,  while  in  reality  there  is  some  resistance.  The  results  are  presented  in 
Appendix  G  (DiCarlo  et  al.,  1997).  (Note  that,  in  this  research,  we  examined  fingering  under 
gravity  with  two  fluids  of  about  equal  viscosity.  In  contrast,  viscous  fingering  is  caused  by  the 
difference  in  viscosity  of  the  two  fluids.)  The  advantage  of  these  studies  was  that  we  were  able 
to  investigate  the  effect  of  contact  angle  without  the  complication  influences  of  changing  the 
properties  of  the  oil  in  the  system. 

In  conclusion,  for  certain  soils,  fingering  can  control  the  movement  of  water  through  oil-contam¬ 
inated  soils  and,  thus,  affect  many  remediation  techniques.  From  theoretical  arguments  and  high 
speed  experimaits,  we  have  been  able  to  determine  properties  such  as  the  size  and  water  content 
in  the  water  fingers.  Early  experiments  suggest  that  surfactants  break  up  the  fingers,  yielding  a 
better  water-oil  mixture  with  positive  effects  for  biological  breakdown  of  oils. 
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5  ABSTRACT 

6 

7  The  fingered  flow  theory  for  air-water  systems  was  extended  to  predict  the  width  of  water 

8  fingers  in  oil-water  systems,  assuming  that  the  pressure  in  the  oil  phase  remains  constant.  It  is 

9  shown  that  the  scaled  diameters  of  oil-water  fingers  can  be  found  from  the  diameters  of  air-water 

10  fingers  using  the  ratio  of  the  surface  tension  and  the  difference  between  the  densities  of  the  oil 

1 1  and  water.  Finger  diameters  measured  for  water  infiltrating  into  light  oil  were  in  agreement 

12  with  finger  diameters  scaled  from  measurements  of  air-water  fingers. 
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1  INTRODUCTION 

2 

3  Interest  in  subsurface  transport  mechanisms  responsible  for  anomalies  from  conventional  models 

4  of  convective-dispersive  flow  has  increased  with  concern  over  groundwater  contamination  by 

5  pesticides  and  non-aqueous  phase  liquids.  One  particularly  important  mechanism  is  fingered 

6  flow  due  to  instabilities.  Studies  of  unstable  flow  of  water  through  unsaturated  porous  media 

7  and  fingering  of  two  liquids  in  Hele-Shaw  cells  have  led  to  various  descriptions  of  these 

8  phenomena,  both  physical  and  theoretical  (Saffman  and  Taylor,  1958).  Chuoke  et  al.  (1959) 

9  extended  the  Hele-Shaw  result  to  soils  using  an  undefined  "effective"  surface  tension. 

10 

1 1  Parlange  and  Hill  (1976)  accounted  for  the  capillarity  of  the  media  by  treating  the  wetting  front 

12  as  diffuse,  rather  than  discontinuous,  in  their  theory  of  wetting  front  instability.  The  moisture 

13  content  and  capillary  pressure  of  the  fingered  flow  resulting  from  such  an  instability  have  been 

14  well  described  for  air-water  systems  in  terms  of  the  initial  moisture  content  and  grain  size 

15  distribution  of  the  media  (Glass  et  al.,  1989a,b;  Baker  and  Hillel,  1990;  Liu  et  al.,  1993;  Selker 

16  et  al.,  1992b,c).  Glass  et  al.  (1991)  apparently  following  the  approach  of  Parlange  and  Hill 

17  developed  general  equations  for  finger  diameters  in  immiscible  fluids,  but  subsequently  Glass 

18  and  Nicholl  (1996)  found  these  equations  did  not  predict  the  observed  finger  behavior  in  oil- 

19  water  systems.  However,  Glass  et  al.  (1991)  assumed  inadvertently  that  the  front  was 

20  discontinuous.  In  this  paper  we  again  use  the  diffuse  front  behavior  of  Parlange  and  Hill  (1976) 

21  but  assume  that  the  oil  provides  little  resistance  to  the  movement  of  the  water  finger  (in 

22  agreement  with  the  conditions  of  our  experiment).  We  show  that  minimum  finger  diameters  for 


3 


1  oil-water  systems  are  thus  obtained  by  taking  into  account  the  oil  phase  density  and  the  oil-water 

2  surface  tension.  We  also  demonstrate  that  viscosity  is  not  a  factor  in  determining  the  minimum 

3  finger  size. 

4 

5  BACKGROUND 

6 

7  This  theory  of  unstable  displacement  of  a  light  oil  by  water  in  porous  media  follows  the 

8  approach  used  to  describe  wetting  front  instability  of  water  in  homogeneous  soils  by  Parlange 

9  and  Hill  (1976).  The  purpose  of  the  theory  is  to  estimate  the  stabilizing  effect  of  surface  tension 

10  (see  Appendix)  and  requires  an  analysis  of  the  transition  from  oil-phase  predominance  to  water- 

11  phase  predominance  across  the  front.  The  earlier  approach  of  Glass  et  al.  (1991)  in  which  the 

12  transition  takes  place  explicitly  across  a  discontinuous  front  is  inconsistent  with  the  physical 

13  process  of  stabilization  introduced  by  the  matric  potential  (surface  tension  effect). 

14 

15  To  formulate  analytically  an  equation  for  the  finger  diameter  in  an  oil-water  system,  we  make 

16  use  of  our  experimental  observations  that  the  oil  pressure,  Po,  within  the  chamber  remains  static 

17  (no  flow)  when  an  unstable  water-oil  front  moves  downward,  indicating  that  the  oil  has  no 

18  difficulty  escaping  ahead  of  the  front  as  the  instability  initially  develops.  This  assumption  is 

19  crucial  to  the  present  formulation.  Only  when  the  finger  penetrates  deeply  into  the  chamber 

20  does  the  oil  pressure  around  the  finger  tip  influence  finger  propagation  and  the  present  theory 

21  break  down.  Defining  Z  as  the  vertical  coordinate  with  the  origin  and  at  the  free  oil  surface 

22  then. 


A. 

PoS 
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(1) 


-  Z  =  0 

1  where 

2  Po  =  density  of  oil 

3  Po  =  oil  pressure 

4  Z  =  vertical  coordinate 

5  g  =  acceleration  of  gravity 

6  and  the  o  subscript  indicates  the  oil  phase. 

7 

8  The  hydraulic  potential  in  the  water  phase,  4>„,  is  defined  in  terms  of  water  head  as: 

4),,  =  —  -  Z;  Z  ^  Z.  +  z,  (2) 

PwS 

9  When  the  oil  pressure  is  static,  Eq.  (2)  can  be  rewritten  by  subtracting  pjp^g  -Zpjp^,  (which 

10  here  according  to  Eq.  (1)  is  zero  because  of  the  choice  of  the  reference  point  for  the  gravity 

11  potential)  as: 

(3) 

I 

12  where 

13  Pw  =  pressure  in  the  water  phase  (units  of  energy  per  unit  volume  of  water) 

14  =  density  of  water 

15  =  matric  potential,  defined  as  (Pw-Po)/pwg 
Zf  +  Zf  =  location  of  the  perturbed  water-oil  front 
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Zf  =  average  position  of  the  water-oil  front 
and  the  w  subscript  indicates  the  water  phase. 

A  water-oil  front  is  defined  as  the  boundary  between  the  oil  and  the  water  where  the  water 
content  is  equal  to  0  (Figure  1).  The  average  position  of  the  front  is  Z  =  Zf.  To  simplify 
the  notation  we  introduce  a  new  Lagrangian  coordinate,  z,  with  its  origin  at  the  average  position 

of  the  front,  thus,  z  =  Z  -  Zf. 


To  find  the  most  unstable  wavelength,  the  front  at  time  t,  is  perturbed  as  (Saffman  and  Taylor, 
1958;  Parlange  and  Hill,  1976): 

Zf  =  W 

For  two-dimensional  gravity  driven  flow,  the  distance  Zf  is  measured  vertically  down  from  the 
average  position  of  the  front,  z  =  0,  in  the  direction  of  the  gravity,  along  the  x  axis  (Figure  1); 
X  is  the  wave  number,  w  is  the  growth  rate  (Glass  et  al.,  1990).  At  the  front  there  must  be  a 
critical  pressure  difference  before  fingering  can  be  initiated.  Hillel  (1987)  referred  to  this 
pressure  difference  as  the  "water  entry"  potential,  In  the  oil-water  system,  h^^  can  be 
expressed  as: 


where 


=  oil  pressure  at  the  front 


(5) 
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water  pressure  at  the  front 


6 


1  and  the  f  subscript  refers  to  properties  measured  at  the  front.  By  definition  IVe  is  equal  to  the 

2  matric  potential  at  the  front,  As  we  will  see  later  the  exact  value  of  h,,^  is  not  important 

3  for  calculation  of  the  finger  width  because  it  is  being  offset  by  other  constants. 

4 

5  The  relationship  between  the  Darcy  velocity  of  the  water,  v,  in  the  z  direction  and  the  velocity 

6  of  the  curved  front  at  the  onset  of  the  instability,  u,  is  (Saffman  and  Taylor,  1958;  Parlange  and 

7  Hill,  1976): 

„  .  „  +  (6) 

dt 

8  Noting  that  the  Darcy  velocities  are  proportional  to  the  hydraulic  potential  gradient,  we  can  find, 

9  by  combining  Eqs.  (4)  and  (6),  that  at  the  front: 


10  where 


dZf 

- 1 - =  M  +  _Z 

dz  dt 


(7) 


11 

12 

13 

14 

15 

16 

17 


the  hydraulic  conductivity  of  the  water  entry  pressure,  which  is  usually 
assumed  to  be  the  saturated  hydraulic  conductivity 


0  =  water  content  at  and  behind  the  front 


0  =  initial  water  content,  which  could  be  equal  to  zero 


Following  Parlange  and  Hill  (1976)  we  expect  the  velocity  of  the  curved  front  and  the  velocity 
of  the  front,  v ,  if  it  had  remained  flat,  to  be  related  as: 
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Combining  Eqs.  (4),  (7),  and  (8)  and  differentiating  with  respect  to  time  and  x  we  find  at  the 
front  that: 

.  (10) 

0  -0  az  '  ' 

W| 

Note  that  the  hydraulic  potential,  behind  and  up  to  the  front  (Figure  1)  satisfies  Laplace’s 
equation  because  in  that  region  the  moisture  content  is  constant  and  equal  to  0^^, 

=0  (11) 
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A  general  solution  of  Eq.  (11)  where  <f>^  is  differentiated  with  respect  to  x  and  Z  and  that 
satisfies  Eq.  (10)  at  the  front  so  that  the  perturbation  decays  as  Z  -oo  is: 


<t>  = 


0...  -e...r 


w, 


v(Z-Z,)  +  ' 


Pwj 


(12) 


where  IXre-<(l-P<yPw)]  is  the  integration  constant  such  that  at  the  front  the  matric  potential  equals 
the  water  entry  potential  (Eq.  (5)).  Equation  (12)  may  be  rewritten  as: 


__^V_  -Pw  ■Pq  _  Pw 

Pw 

Evaluating  Eq.  (13)  at  the  front  where  z  =  Zf  and  noting  that  for  Zf  <  <  1  the  exponential  Xz 
becomes  1,  we  find  that: 


Pw  Pq 


(14) 


The  most  unstable  wave  number  (in  the  oil-water  system),  Xqu,  is  obtained  from  Eq.  (14)  by 
calculating  dw/dX  and  setting  it  equal  to  zero,  i.e.,  co/X^a  =  oiou^oa-  Equation  (14)  can  then  be 
rewritten  for  the  most  unstable  wave  number  as: 


ef\  v  oil  oil  y 

-d.  Ac 

uy  if  bf  bf 


(15) 


In  order  to  compare  the  fingers  in  the  oil-water  system  with  air-water  fingers,  v/e  take  as 
reference  the  X,ar  for  Ih®  air-water  system  which,  when  v  ^  0  (Parlange  and  Hill,  1976),  is  given 
by: 


9 


1  where  for  the  air-water  system  is  defined  for  the  same  fluid  contents  as  in  the  oil-water 

2  system: 


3 

4 

5 


Finally,  since  the  finger  width,  d,  is  inversely  proportional  to  X,  d  «  ir/X  (Parlange  and  Hill, 
1976,  Figure  1),  Eqs.  (15)  and  (16)  give  the  ratio  and  finger  diameters  in  the  oil-water  system. 


doa,  and  in  the  air- water  system,  d^,  as: 


^oU 


(18) 


6  If  the  assumption  of  oil  moving  freely  from  the  wetting  front  away  is  incorrect,  we  expect  the 

7  finger  to  be  wider  than  predicted.  Thus,  Eq.  (18)  predicts  the  narrowest  finger  expected  to  be 

8  observed  in  the  oil- water  system.  Finally,  if  the  term  with  v  is  negligible  (which  it  will  be  for 

9  our  case),  Eq.  (18)  differs  from  the  result  of  Glass  et  al.  (1991)  by  a  term  which 

10  appears  in  the  denominator  (see  their  Eq.  (17)). 

11 


12 
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1  MATERIALS  AND  METHODS 

2 

3  Two  sets  of  experiments  were  performed  with  oil-water  systems.  The  first  set  was  performed 

4  with  20-30  sand  in  a  relatively  large  chamber  and  gave  initial  estimates  of  finger  diameter.  To 

5  confirm  that  these  experiments  were  not  an  artifact,  a  second  set  of  experiments  were  performed 

6  at  a  later  date  with  12-20  sand  in  which  the  methods  by  which  the  water  was  added  and  the  oil 

7  was  removed  were  changed.  Moreover,  the  oil  pressure  was  measured  to  test  the  assumption 

8  in  the  model  of  static  pressure  in  the  displaced  phase. 

9 

10  In  experimental  set  I  a  glass  walled  chamber  with  interior  dimensions  of  80  cm  by  50  cm  by  1 

11  cm  was  used  (Figure  2).  For  each  trial,  the  chamber  was  loaded  with  white  silica  sand  which 

12  had  been  sorted  between  standard  #20  and  #30  sieves,  cleaned  in  the  manner  of  Glass  et  al. 

13  (1989a)  to  remove  organic  contaminants  and  air  dried.  Following  vibrational  packing  (Selker, 

14  1992a),  vacuum  degassed  oil  (Soltrol  220,  Philips  Petroleum  Co.)  was  introduced  slowly  from 

15  the  bottom  of  the  chamber  until  it  ponded  above  the  sand  surface.  Degassed,  distilled  water  was 

16  applied  to  the  oil-saturated  slab  by  two  techniques:  At  a  rate  of  0.2  ml/cm^/min,  using  a  cam- 

17  driven  oscillating  dripper  and  2.0  ml/min  via  a  stationary  dripper  located  near  the  center  of  the 

18  chamber  above  the  oil  surface.  Overflow  drains  at  both  sides  of  the  chamber  and  an  underdrain 

19  manifold  were  arranged  to  maintain  the  oil-free  surface  at  1  cm  above  the  sand. 

20 

21  In  experimental  set  n  a  chamber  was  used  with  interior  dimensions  of  55  cm  by  33  cm  by  1  cm. 

22  The  chamber  was  filled  and  saturated  with  oil  in  a  manner  similar  to  experimental  set  I.  From 


11 


1  five  ports  at  the  bottom  of  the  chamber  the  oil  was  pumped  at  a  rate  of  2  ml/min  into  a 

2  stoppered  Ehrlenmeyer  flask  with  water.  The  water  which  was  forced  out  of  the  flask  was 

3  applied  to  the  top  of  the  chamber  0.5  cm  below  the  oil  surface  via  a  small  glass  tube.  In  this 

4  case  no  free  water  surface  formed  on  top  of  the  sand  and  fingers  started  immediately.  Pressures 

5  in  the  oil  phase  were  measured  at  11,  19,  and  36  cm  below  the  oil-free  surface  with  1  cm 

6  diameter,  fast  responding  tensiometers  mounted  flush  with  the  chamber  wall  (Rimmer  et  al., 

7  1996).  Of  the  20  experiments  which  were  performed,  we  reported  finger  diameters  from  the 

8  six  experiments  for  which  the  chamber  loading  technique  was  adequately  uniform  and  oil 

9  pressure  could  be  measured  consistently. 

10 

11  Finger  dimensions  were  determined  by  coloring  the  water  with  50  mg/1  FD&C  blue  dye  #1  and 

12  the  oil  with  0.2  mg/1  Sudan  IV  red  dye  (Aldrich  Chemical  Co.).  Traces  were  made  from  video 

13  images  (experimental  set  I)  or  directly  with  transparent  sheets  on  the  chamber  wall  (experimental 

14  set  II).  Finger  diameters  were  measured  when  the  finger  tip  did  not  change  any  more  in  time 

15  (experimental  set  I)  or  was  near  or  at  the  bottom  of  the  chamber  (experimental  set  II). 

16 

17  RESULTS 

18 

19  Two  types  of  results  were  obtained:  Finger  diameters  for  experimental  sets  I  &  n  and  the  oil 

20  pressure  for  experimental  set  II. 

21 
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1  Examples  of  finger  traces  for  trials  from  each  of  the  experimental  sets  in  which  the  flow  rate 

2  was  near  2  cm’min  *  are  shown  in  Figures  3  and  4.  Finger  widths  for  both  experimental  sets 

3  are  shown  in  Figures  5  and  6  for  flow  rates  comparable  to  air-water  experiments  which  we  have 

4  performed  previously.  In  general,  the  average  finger  diameter  was  nearly  constant  for  the  first 

5  15  cm  and  started  to  increase  gradually  below  that  depth.  For  experimental  set  II  the  increase 

6  in  finger  width  was  even  more  notable  near  the  bottom  of  the  chamber.  In  experimental  set  I 

7  the  fingers  did  not  reach  the  bottom  of  the  chamber.  Instead,  new  fingers  formed  when  the 

8  fingers  were  near  30  cm  depth.  The  finger  width  for  experimental  set  I  with  the  20-30  sand  was 

9  4.8  cm  with  a  standard  deviation  of  1  cm.  For  experimental  set  II,  the  width  was  initially  4.7 

10  cm,  also  with  a  standard  deviation  of  1  cm  and  then  below  30  cm  increased  to  a  width  of  8.6 

11  cm  with  a  standard  deviation  of  3  cm. 

12 

13  An  example  of  the  oil  pressure  measured  in  experimental  set  n  for  tensiometers  at  11,  19,  and 

14  36  cm  below  the  oil-free  surface  is  given  in  Figure  7.  The  pressure  in  each  of  the  tensiometers 

15  remained  nearly  constant  once  the  sand  slab  had  been  saturated  with  oil.  A  slight  pressure 

16  increase  was  found  at  all  tensiometers  throughout  the  course  of  the  experiment,  reflecting 

17  displacement  of  oil  by  water  during  the  course  of  the  experiment.  More  importantly  the  oil 

18  pressure  difference  between  the  tensiometers  was  constant,  regardless  of  the  finger  having  passed 

19  by  the  tensiometers. 

20 


21 
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1  DISCUSSION 

2 

3  Fingers  grew  at  constant  but  small  velocities  for  each  experiment.  However,  in  contrast  to 

4  fingers  observed  in  similar  experiments  in  the  air-water  system  (Selker  et  al.,  1992b, c)  finger 

5  velocities  (although  internally  constant)  were  more  variable  for  different  application  rates  and 

6  chamber  regions.  For  stationary  water  application  rates  of  2  ml/min,  finger  velocities  ranged 

7  from  0.019  to  0.022  cm/sec.  For  linear  water  application  of  0.2  ml/min/cm^,  the  average  finger 

8  velocities  were  lower,  varying  between  0.005  and  0.009  cm/sec.  In  all  cases  the  influence  of 

9  the  term  with  v  in  Eq.  (18)  was  negligible. 

10 

11  A  critical  assumption  in  the  development  of  the  unstable  flow  theory  for  the  oil-water  system 

12  is  that  the  oil  pressure  remains  static  throughout  the  experiment.  This  assumption  is  validated 

13  by  the  static  difference  in  oil  pressure  between  any  two  tensiometers  (Figure  7),  e.g.:  The 

14  pressure  difference  between  tensiometers  3  and  4  remained  static  at  approximately  13.6  cm, 

15  which  is  equivalent  to  the  difference  in  oil  head  between  them,  17  cm,  multiplied  by  the  density 

16  of  Soltrol,  0.80  g/cm^.  Taking  into  account  the  specific  weight  of  oil  and  a  distance  of  2.5  cm 

17  between  them,  show  that  the  pressure  is  static.  Hence  the  static  oil  pressure  assumed  in  the 

18  model  was  valid. 

19 

20  The  predicted  finger  diameter  for  the  oil-water  system  is  based  on  the  physical  constants  and  the 

21  velocities  of  the  fingers  in  each  system,  assuming  that  the  diameter  and  velocity  of  an  individual 

22  finger  are  constant.  The  physical  constants  required  for  scaling  fingers  between  systems  (Eq. 
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1  (18))  are  the  fluid  densities  and  interfacial  tensions  for  each  of  the  two  fluids  within  the  porous 

2  medium,  the  hydraulic  conductivity  of  the  sand  and  the  porosity  of  the  media.  Water  and 

3  Soltrol-220,  have  densities  of  1.00  and  0.80  g*cm'^,  respectively,  at  20°  C  (Cary  et  al.,  1989). 

4  The  water-air  interfacial  tension  is  73  g*sec'^,  however,  interfacial  tension  is  neither  well 

5  defined  nor  constant  over  time  (Schroth  et  al.,  1995).  Schroth  et  al.  (1996)  found  the  Soltrol- 

6  water  interfacial  tension  between  40  and  46  g  •  sec'^.  Cary  et  al.  (1989)  also  reported  a  surface 

7  tension  for  Soltrol  and  water  of  45  g  •  sec'^.  We  repeatedly  measured  the  interfacial  tension  for 

8  the  Soltrol  and  water  used  by  us  (with  and  without  the  dyes)  to  be  32  to  33  g  •  sec'^  with  a  Fisher 

9  Scientific  Model  21  ring  tensiometer.  The  hydraulic  conductivities  for  the  12-20  and  20-30  sand 

10  were  30  cm/min  (Schroth  et  al.,  1996)  and  20  cm/min,  respectively.  The  porosities  are  0.36 

11  cmVcm^  (Selker  et  al.,  1992a;  Liu  et  al.,  1993).  Additional  parameters  required  for  calculating 

12  finger  diameter  by  Glass  et  al.  (1991)  are  the  viscosities  of  water  and  Soltrol  which  are  0.01  and 

13  0.047  g* sec-1  •cm  ‘,  respectively  (Cary  et  al.,  1989). 

14 

15  The  finger  diameters  for  the  top  15  cm  of  our  experiments  as  well  as  others  (Selker  et  al.,  1992; 

16  Liu  et  al.,  1994;  Damault  et  al.,  1996;  Hendrickx  and  Yao,  1996;  DiCarlo  et  al.,  1997;  Rimmer 

17  et  al.,  1996)  are  shown  in  Table  1.  The  difference  in  finger  widths  between  experimental  sets 

18  where  point  source  water  application  and  those  where  uniform  linear  water  applications  were 

19  used,  is  remarkably  robust.  Taking  all  finger  diameters  into  account  we  find  the  average  oil- 

20  water  finger  width  in  the  12-20  sands  to  be  2.9  times  that  of  air-water  fingers.  For  the  20-30 

21  sand  it  is  2.3  times  greater. 


22 
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1  To  calculate  the  theoretical  scaling  for  finger  diameters  in  air-water  and  oil-water  systems,  we 

2  neglect  the  velocity  term  because  the  finger  velocity  is  much  less  than  the  saturated  conductivity 

3  of  the  media.  We  will  calculate  the  scaling  factor  for  Eq.  (18)  and  compare  this  with  the  scaling 

4  of  Glass  et  al.  (1991)  and  Chuoke  et  al.  (1959).  The  scaling  calculated  for  Glass  et  al.  (1991) 

5  is  larger  by  a  factor  of  over  that  in  Eq.  (18).  Chuoke  et  al.’s  (1959)  scaling  (as 

6  formulated  by  Glass  et  al.,  1991)  is  similar  to  Eq.  (18)  but  to  the  power  0.5.  Because  of  the 

7  uncertainty  of  the  Soltrol- water  interfacial  tension,  the  calculated  scaling  factors  are  shown  in 

8  Figure  8  for  the  full  range  of  the  aforementioned  interfacial  tensions.  Also  plotted  are  the 

9  observed  ranges  of  scaling  between  the  oil- water  and  air  water  fingers  of  Table  1.  Whereas  that 

10  Eq.  (18)  predicts  the  scaling  in  minimum  finger  diameters  quite  well,  the  scaling  proposed  by 

11  Glass  et  al.  (1991)  overpredicts  and  that  of  Chuoke  et  al.  (1959)  slightly  underpredicts  finger 

12  width  for  oil-water  systems. 

13 

14  In  both  experimental  sets  the  finger  width  increased  below  15  cm  (Figures  5  and  6).  In  earlier 

15  experiments  with  air-water  systems  (Glass  et  al.,  1989a,  1990;  Selker  et  al.,  1992a,b;  Liu  et  al., 

16  1993,  1994)  this  widening  at  the  bottom  was  not  observed.  There  are  several  explanations 

17  possible  for  the  widening.  First,  the  nature  of  the  flow  field  produced  in  the  oil-water  system 

18  differed  from  that  in  the  air- water  system.  A  second  less  likely  explanation  (because  the  theory 

19  takes  this  into  account)  is  the  difference  in  pressure  profile:  The  air  pressure  remains 

20  atmospheric  while  in  the  oil  the  pressure  increases  with  depth.  One  would  expect  a  continuous 

21  increase  with  depth  for  the  oil- water  fingers  which  was  not  the  case.  A  final  explanation  for 

22  variability  in  finger  width  in  the  oil-water  system  is  that  the  greater  density  of  the  displaced  fluid 
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1  in  the  oil-water  system  as  compared  to  the  air-water  system  reduces  the  driving  force  and 

2  renders  the  finger  width  more  sensitive  to  anisotropy  in  the  sand  slab.  Currently,  we  are 

3  investigating  the  density  differences  in  the  chambers  with  x-rays  at  the  Cornell  High  Energy 

4  Synchrotron  Source. 

5 

6  CONCLUSIONS 

7 

8  Unstable  wetting  front  theory  was  developed  for  oil-water  systems.  The  theory  is  different  from 

9  air- water  systems  (where  the  air  pressure  remains  atmospheric)  in  that  the  pressure  in  the  oil  is 

10  increasing  with  depth.  Nevertheless,  if  the  oil  pressure  does  not  change  as  the  finger 

1 1  approaches,  then  the  oil  has  no  difficulty  escaping  ahead  of  the  wetting  front.  According  to  the 

12  presented  theory,  the  minimum  finger  widths  in  oil-water  systems  may  be  scaled  to  air-water 

13  fingers  using  the  densities  and  interfacial  tensions  of  the  fluids,  whereas  fluid  viscosities  need 

14  not  be  considered  for  scaling  widths  of  gravity  driven  fingers.  Experiments  showed  that  fingers 

15  in  oil-water  systems  are  consistent  with  the  theoretical  assumptions  of  constant  finger  velocity 

16  and  width  made  for  fingers  in  air-water  systems.  Moreover,  the  finger  diameter  observations 

17  agreed  well  with  the  theoretically  predicted  behavior.  However  oil- water  fingers  are  apparently 

18  more  sensitive  to  packing  heterogeneities  in  the  porous  media  than  are  air- water  fingers.  The 

19  increase  of  pressure  head  with  depth  in  the  oil- water  system  may  also  contribute  to  variations 

20  in  finger  dimensions.  A  more  complex  theoretical  model  is  required  to  describe  this  case. 

21  Hence,  further  research  is  needed  to  describe  the  forces  associated  with  flow  of  oil  and  the  effect 

22  on  fingered  flow. 
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1  APPENDIX 


To  find  ttou*  we  introduce  water  in  oil-saturated  soil  at  a  constant  rate,  Q,  per  unit  length  of  a 
vertical  line  source.  Then,  for  any  distance,  r,  from  the  line  source  assuming  that  oil  moves 
freely  ahead  of  the  front: 


6  Integration  of  Eq.  (Al)  yields: 


Q  =  -2nrK 


=  Ik,. 


V  e, 


since  r;  (where  0=0,  Figure  6)  and  rf  (where  0=0  )  do  not  differ  by  much.  Equation 

w^  Hy 

(A2)  can  be  rewritten  by  expanding  the  logarithmic  term  and  neglecting  higher  order  terms: 


r2it  7 

r.  =  r  +  — —  f  K  —dQ 
'  *  Q  I  dQ 


9  From  conservation  of  mass: 


-  Kyj  ”  Q<  <■ 

10  Then,  following  Parlange  and  Hill  (1976)  we  obtain  (differentiating  Eq.  (A3)  with  time): 


u  =  V  - 


1 1  where  u  =  dr,/dt  is  the  velocity  of  the  curved  front  and  v  »  dr/dt  or  the  velocity  in  the  absence 
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of  surface  tension  effects: 
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J  ’"0-0 
Ow,  “y  *"f 


1  assuming  the  contact  angles  are  the  same  then  for  the  oil-water  and  the  air-water  system 

2  scales  as  we  can  write  Eq.  (A6)  as: 


= 


r,2  ^Od 

^air - 


where  S.;,  is  the  sorptivity  in  the  air-water  system.  Equation  (8)  is  identical  to  Eq.  (A5), 
replacing  rf  by  the  curvature  of  the  front. 
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Table  1 .  Summary  of  observed  two-dimensional  finger  widths  for  oil-water  and  air-water  fluid 
pairs  in  12-20  and  20-30  sand. 


Sand  type:  12-20  20-30 


Fluids  (initial-invading):  oil- water 

air-water 

oil-water 

air-water 

FINGER  DIAMETER  (cm) 

This  study  4.7  (1.0)' 

4.8  (1.0) 

Rimmer  et  al.  (1996) 

— 

4.6 

— 

Damault  et  al.  (1996)  5.3 

— 

— 

— 

DiCarlo  et  al.  (1997)  — 

1.7  (0.2) 

— 

2.3  (0.3) 

Hendrickx  &  Yao  (1996)^  — 

1.6 

— 

2.6 

Liu  et  al.  (1994) 

— 

— 

2.25 

Selker  et  al.  (1992a)  — 

— 

— 

1.7 

OBSERVED  SCALING 

BETWEEN  FLUID  PAIRS 

2.8 

2.1 

‘  The  number  in  parentheses  is  the  standard  deviation  when  known. 

^  Three-dimensional  finger  width  converted  to  a  two-dimensional  width  by  multiplying  by  ir/4.8. 
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1  LIST  OF  nCURES 

2 

3  Figure  1:  Perturbed  water-oil  front. 

4 

5  Figure  2:  Schematic  diagram  of  the  light  chamber. 

6 

7  Figure  3:  Traces  of  three  fingers  initiated  by  point  source  application  of  2  ml/min.  Each  trace 

8  is  labeled  by  time,  in  minutes  (experimental  set  I). 

9 

10  Figure  4:  Traces  of  three  fingers  initiated  by  point  source  application  of  2  ml/min.  Each  trace 

11  is  labeled  by  time,  in  minutes  (experimental  set  II). 

12 

13  Figure  5:  Observed  finger  width  with  depth  for  experimental  set  I.  Each  symbol  represents  one 

14  finger.  The  solid  line  is  the  running  average. 

15 

16  Figure  6:  Observed  finger  width  with  depth  for  experimental  set  II.  The  solid  line  is  the 

17  running  average. 

18 

19  Figure  7:  Oil  pressure  readings  for  tensiometers  located  at  11,  19,  and  36  cm  below  the  surface 

20  (experimental  set  II). 

21 
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1  Figure  8:  Calculated  and  observed  scaling  of  finger  diameters  in  air-water  and  oil-water 

2  systems. 

3 


4  Figure  9:  Water-oil  front  for  a  vertical  line  source. 
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1  ABSTRACT 

2 

3  The  complexity  of  simultaneous  flow  of  water  and  non-aqueous  phase  liquids  is  largely 

4  unappreciated  because  few  techniques  permit  accurate  quantitative  measurement  of  water  and  oil 

5  contents  in  rapidly  changing  flow  fields.  High  intensity  x-rays  were  used  at  the  Cornell  High 

6  Energy  Synchrotron  Source  (CHESS)  to  obtain  rapid,  accurate,  and  non-destructive  quantitative 

7  measurements  of  the  changing  fluid  contents  in  a  porous  media  during  infiltration  events. 

8  Concomitant  temporal  pressure  measurements  were  obtained  for  each  liquid  phase  using  rapid 

9  responding  tensiometers.  The  system  was  used  for  measuring  temporal  volumetric  fluid  content 

10  changes  during  a  water  finger  infiltration  into  sand  saturated  with  a  NAPL  (Soltrol-220)  in  a 

11  two-dimensional  soil  slab.  The  fluid  content  distribution  of  a  finger  in  the  oil-water  system  was 

12  found  to  be  similar  to  air-water  systems.  The  hysteretic  constitutive  relationship  between 

13  pressure  and  the  content  was  developed  from  the  data.  The  relationship  was  used  to  explain  why 

14  the  finger  did  not  widen  behind  the  tip,  and  why,  upon  re-infiltration,  water  followed  the 

15  previously  established  path.  These  findings  are  relevant  for  cleanup  of  oil  contaminated  sites 

16  because  it  aids  in  the  understanding  of  hydrologic  control  which  is  an  essential  component  of 

17  cost-effective  in  situ  remediation. 
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1  INTRODUCTION 

2 

3  Groundwater  contaminants  are  introduced  at  or  near  the  soil  surface  via  spills,  leakage  from  fuel 

4  tanks,  or  overfilling  of  tanks.  Many  of  these  contaminants  are  organic  compounds  of  low  water 

5  solubility,  existing  as  separate  liquid  phases  in  soil  pore  spaces  (Dane  et  al.,  1992).  The 

6  problem  is  severe.  For  example,  the  volume  of  contaminants  at  130  Department  of  Energy 

7  (DOE)  facilities  alone  exceeds  2,500  billion  liters  (Department  of  Energy,  1995;  Wilkins  et  al., 

8  1995).  The  anticipated  life-cycle  cleanup  costs  are  expected  to  be  greater  than  300  billion 

9  dollars  and  will  take  at  least  100  years  using  existing  technologies. 

10 

11  In  situ  remediation  technologies  are  the  only  viable  approaches  with  potential  for  controlling  the 

12  huge  costs  and  extended  times  required  to  clean  up  these  sites.  However,  cleanup  procedures 

13  are  often  instituted  and  operated  with  limited  knowledge  of  both  the  physical  and  chemical 

14  processes  governing  contaminant  transport  (Wilkens  et  al.,  1995).  To  be  effective  complete 

15  hydrologic  control  of  the  subsurface  at  the  remediation  site  is  critical.  Hydrologic  control  at 

16  NAPL  contaminated  sites  can  be  improved  by  accurately  measuring  water  and  oil  contents  in 

17  rapidly  changing  flow  fields  (Mercer  and  Cohen,  1990). 

18 

19  Early  observations  of  oil  and  water  flow  were  made  in  Hele-Shaw  flow  cells  by  Saffman  and 

20  Taylor  (1958)  and  Chuoke  et  al.  (1959).  Currently,  most  real-time,  non-intrusive  measurements 

21  of  oil  and  water  contents  in  porous  media  are  made  using  gamma-  and  x-ray  attenuation 

22  techniques.  For  example,  a  dual  energy  gamma  radiation  method  was  described  by  Hopmans 
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1  and  Dane  (1986),  and  was  applied  to  study  saturation  changes  during  monotonic  liquid  drainage 

2  (Lenhard  et  al.,  1988;  Dane  et  al.,  1992).  The  x-ray  absorption  (film-based  radiography) 

3  technique  described  by  Tidwell  and  Glass  (1994)  can  be  used  to  obtain  full  field  oil  and  water 

4  contents.  However,  the  conventional  x-  and  gamma-ray  methods  are  limited  by  the  relatively 

5  low  intensity  of  Am  241,  Cs  137,  and  x-ray  tube  sources,  which  do  not  permit  accurate 

6  monitoring  of  flow  phenomena  occurring  in  porous  media  on  temporal  scales  of  seconds. 

7 

8  Other  methods  for  determining  oil  and  water  fluid  contents  and  that  permit  measurement  of 

9  rapidly  changing  flow  fields  is  the  full-field  visualization  techniques,  described  by  Van  Geel  and 

10  Sykes  (1994)  and  Damault  et  al.  (1996)  which  involve  dying  the  fluids  and  analyzing  the  visual 

11  information  with  computer  graphics  techniques.  Van  Geel  and  Sykes  (1994)  dyed  the  oil  red, 

12  and  analyzed  the  reflected  light.  This  determined  the  fluid  contents  near  the  glass  wall  which 

13  may  not  represent  that  of  the  chamber.  Darnault  et  al.  (1996)  colored  the  water  with  0.(X)5% 

14  blue  dye  and  recorded  the  transmitted  light  from  a  bank  of  high  frequency  fluorescent  lights 

15  behind  the  slab  chamber.  The  color  (or  hue)  of  the  images  was  related  to  moisture  content 

16  values.  However,  scattering  near  the  edges  made  it  more  difficult  to  see  sharp  boundaries. 

17  Moreover,  the  sands  need  to  be  translucent  (Tidwell  and  Glass,  1994).  While  light  transmission 

18  techniques  allow  monitoring  of  many  array  points  in  a  relatively  large  two-dimensional  section 

19  of  a  soil  slab,  x-ray  attenuation  permits  more  accurate  water  content  measurements  in  a  small 

20  section  of  any  soil  type. 


21 
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1  A  final  destructive  method  for  visualizing  oil  and  water  contents  is  by  using  oils  that  solidify. 

2  After  solidification,  the  columns  are  sectioned  and  the  distribution  of  oil  studied  (Conrad  et  al., 

3  1992;  Powers  et  al.,  1992). 

4 

5  The  objective  of  this  paper  is  to  present  a  method  using  synchrotron  x-rays  that  allows  rapid 

6  quantitative,  and  precise  measurements  of  fluid  contents  in  oil-water  systems.  Synchrotron 

7  radiation  can  be  up  to  10,000  times  as  intense  than  conventional  x-ray  sources,  which  decreases 

8  errors  and  counting  times.  Synchrotron  beam  lines  are  rapidly  increasing  and  are  available  to 

9  the  general  scientific  community. 

10 

11  To  illustrate  the  application  of  the  synchrotron  x-ray  technique,  we  experimented  with  fingered 

12  flow  through  a  soil  initially  saturated  with  a  light  oil  (Soltrol  220).  Soltrol  has  negligible  water 

13  solubility,  very  low  volatility  at  room  temperatures,  as  well  as  a  low  health  hazard.  We 

14  measured  fluid  contents  and  pressures  in  both  phases  using  phase  sensitive  tensiometers.  This 

15  study  is  relevant  because  fingered  flow  generated  by  chaotic  disturbances  in  a  flow  field  is  one 

16  of  the  factors  complicating  hydrologic  control  during  in  situ  cleanup,  and  can  render  the  use  of 

17  in  situ  treatment  remedies  ineffective  because  of  hydrodynamic  bypassing  of  subsurface 

18  pollutants.  With  the  result  of  the  synchrotron  experiments,  the  forming  and  persistence  of  these 

19  bypass  flow  paths  are  explained.  It  is  shown  that  because  of  hysteresis  in  the  constitutive 

20  relationships,  these  paths  are  in  equilibrium  with  the  surrounding  oil  fluid.  Finally,  the 

21  synchrotron  and  other  measurements  are  used  to  scale  the  finger  diameter.  By  taking  soil  and 

22  fluid  properties  into  account  we  calculate  the  dimensions  of  the  bypass  flow  paths. 
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BACKGROUND 


Unstable  Flow  Hienoniena 


Instabilities  involving  two  immiscible  fluids  (known  as  bypass  flow)  have  been  observed  in  the 
oil  industry  for  some  time  (Craft  and  Hawkins,  1959).  General  stability  criteria  for  two  liquids 
have  been  given  by  Saffman  and  Taylor  (1958)  and  Chuoke  et  al.  (1959)  and  specific  criteria 
for  air-water  systems  in  soil  by  Parlange  and  Hill  (1976),  Raats  (1973),  Philip  (1975),  and  Hillel 
and  Baker  (1988).  In  general,  the  flow  will  be  unstable  if  the  pressure  difference  between  the 
two  liquids  is  greatest  where  the  fluid  has  intruded  the  farthest.  This  condition  can  be  met  when 
a  less  viscous  fluid  intrudes  into  a  more  viscous  fluid,  and/or  when  a  denser  fluid  intrudes  with 
gravity  into  a  less  dense  fluid. 


Diameters  of  fingers  in  air-water  systems  were  calculated  by  Parlange  and  Hill  (1976),  Glass 
et  al.  (1989a, b),  and  Liu  et  al.  (1994a).  The  equation  derived  from  Parlange  and  Hill  (1976) 
and  by  Liu  et  al.  (1994a)  is  the  simplest  and  applies  to  well  sorted  sands  (i.e.,  those  with  narrow 
particle  size  distributions)  where  the  slope  of  the  wetting  branch  of  the  soil  characteristic  curve, 
dh^^/d0^,  is  straight  over  the  portion  that  K(0)  is  decreases  rapidly  with  moisture  content  and 
is  written  as  by  Liu  et  al.  (1994a)  as: 


2  K  0. 


dh 


afw 


^atw 


dQ^ 


n  +  1.5 


(1) 
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where  #c  is  equal  to  x  for  two-dimensional  fingers  and  4.8  for  three-dimensional  fingers  (Glass 
et  al.,  1990)  and  6f  is  the  moisture  content  at  the  finger  tip.  17  may  be  found  as  (Brooks  and 
Corey,  1964): 

r,  =  1  +  3  (2) 

A 

where  X  is  a  positive  index  related  to  the  pore  size  distribution  of  soil  with  values  ranging  from 
2  to  7  resulting  in  values  between  3.3  and  4. 


Based  on  Hele-Shaw  type  of  experiments,  finger  width  for  two-phase  immiscible  liquid  systems 
was  given  by  Chuoke  et  al.  (1959)  who  assumed  that  a  relationship  analogous  to  the  Laplace- 
Young  relationship  exists  between  the  macroscopic  frontal  curvature  and  the  pressure  jump 
across  the  front  as  (Glass  et  al.,  1991b): 


3a* 


10.5 


«(Pw  -  Pop  - 


(3) 


A  difficulty  with  Eq.  (5)  is  the  "effective"  surface  tension  a*  which  has  not  been  defined, 
although  Homsy  (1987)  associated  it  with  a  modified  capillary  number,  g  is  the  gravity,  p„  is 
specific  weight  of  water,  is  specific  weight  of  oil,  q  is  flux,  and  K^f  is  the  conductivity  in  the 
finger. 


Chandler  et  al.  (1996)  assumed  little  resistance  to  moving  of  the  oil  and  related  finger  widths 
for  air-water  systems,  d^„  and  oil-water  systems  do,„  as: 
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Equation  (3),  assuming  that  a  scales  like  the  surface  tension,  and  Eq.  (5)  predict  different  finger 
scaling  behavior  with  fluid  density,  fluid  interfacial  tension,  and  the  coarseness  of  the  porous 
media.  In  the  following  analysis  we  will  compare  widths  of  slow  moving  fingers  and,  thus, 
minimize  the  effects  of  the  flux  dependency.  Noting  that  matric  potential  scales  as  the  inverse 
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1  of  the  microscopic  length  (Miller  and  Miller,  1956)  we  can  also  write  the  finger  scaling  in  the 

2  form: 


d  =  d*c  (6) 

G  f,* 

where  the  superscript  *  is  the  "reference  soil  with  reference  fluids",  G  =  (l-pn«/Pw)  accounts 
for  different  fluid  densities,  and  D50  is  mean  particle  diameter.  The  subscript  w  is  for  wetting 
and  subscript  nw  for  nonwetting.  The  coefficient,  b,  equals  1  for  scaling  according  to  Eq.  (4) 
and  b  is  0.5  using  Eq.  (3)  of  Chuoke  et  al.  (1959).  The  factor  c  is  a  correction  factor  for 
difference  in  moisture  content  in  the  finger  tip  and  type  of  chamber  used  (2  or  3  dimensional). 
By  introducing  scaling  factors  for  the  particle  size,  fp,  for  gravity,  fo,  and  surface  tension,  f„, 


Eq.  (6)  can  be  rewritten  as: 


10  where 


f  ^ 

Jd  rx 


Jo  ~  . 
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C  =  CgC^ 


(11) 


In  Eq.  (11)  the  k  is  the  coefficient  introduced  by  Glass  et  al.  (1990),  for  2-d  (k=x)  and  for  3-d 
(ic=4.8).  The  correction  for  moisture  content  in  the  finger  tip,  c*  is  based  on  fluid  content 
measurements.  The  effect  of  the  moisture  content  on  the  finger  is  not  specified  in  Chuoke’s  Eq. 
(3).  The  correction  factor  tends  to  be  small  and,  therefore,  the  high  precision  is  not  necessary. 


Oil  and  Water  Content  Measurement  Using  Synchrotron  X-Rays 


Bierck  et  al.  (1988)  described  a  method  using  synchrotron  radiation  to  monitor  consolidation  of 
clays  under  seepage  stress.  The  high  intensity  of  synchrotron  radiation  permitted  acquisition  of 
transmitted  x-ray  intensities  through  relatively  thick  regions  (on  the  scale  of  centimeters)  during 
short  time  intervals  (in  a  second  or  less).  This  method  was  adapted  as  described  by  Liu  et  al. 
(1993)  to  measure  water  content  changes  during  rapid  water  finger  movement  through  sand. 


In  a  two-fluid  system,  here  oil  and  water,  the  water  content  can  be  directly  related  to  the 
attenuation  of  the  x-rays,  ln(I/lo),  where  I  and  Iq  are  the  transmitted  and  incident  radiation 
intensities,  respectively.  Straightforward  calculations  (Liu  et  al.,  1993;  Lenhard  et  al.,  1988) 
yield  the  water  content  in  terms  of  the  measured  x-ray  attenuation,  as: 
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(12) 


Here  6^  is  the  volumetric  water  content  (cm^/cm^),  I„t  is  the  transmitted  intensity  through  the 
initially  oil-saturated  porous  media,  x  is  the  thickness  of  the  chamber  (cm),  and  a*  and  oc^  are 
the  attenuation  coefficients  (cm  *)  for  the  water  and  oil.  The  volumetric  oil  content  is  given 
by: 

,  (13) 

where  n  is  the  total  pore  volume. 


MATERIALS  AND  METHODS 


Experiments  were  carried  out  at  the  F-2  beam  line  of  the  Cornell  High  Energy  Synchrotron 
Source  (CHESS),  which  provides  an  intense,  high  energy  x-ray  beam  (Batterman  and  Ashcroft, 
1979).  Using  a  Si(311)  monochromator,  the  beam’s  energy  was  tuned  to  33.2  KeV.  This 
energy  level  was  chosen  because  it  was  slightly  above  the  iodine  absorption  edge  and  increased 
the  contrast  between  oil  and  water  doped  with  30  g/1  of  potassium  iodine.  Subsequently,  the 
beam  passed  through  a  set  of  slits,  a  sealed  Xe  ionization  chamber,  the  sample  chamber,  and 
another  sealed  Xe  ionization  chamber.  Ionization  chambers  produce  a  current  directly 
proportional  to  the  x-ray  intensity  that  passes  through  them.  The  slits  defined  a  beam  cross- 
section  of  0. 1  cm  high  and  0.6  cm  wide.  For  each  ion  chamber,  the  current  was  converted  into 
a  frequency,  and  acquired  by  a  computer  which  automatically  recorded  the  relative  incident  and 
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1  transmitted  radiation.  The  sample  chamber  was  mounted  on  a  computer  controlled  translation 

2  stage,  allowing  monitoring  of  fluid  content  changes  at  different  positions  within  the  chamber. 

3 

4  The  soil  slab,  illustrated  in  Figure  1,  was  a  chamber  54  cm  high,  23  cm  wide,  and  0.95  cm 

5  thick  with  0.95  cm  polycarbonate  walls,  equipped  with  four  water  and  four  oil  tensiometer  ports. 

6  The  chamber  was  positioned  vertically  such  that  the  x-rays  struck  (and,  thus,  monitored  the  fluid 

7  contents)  2  cm  below  the  lower  oil  tensiometers  and  30  cm  below  the  top  of  the  chamber.  This 

8  initial  position  is  indicated  by  a  dotted  line  in  Figure  1 . 

9 

10  Two  types  of  sieved  silica  sand  were  used  to  pack  the  chamber;  either  "20/30  sand"  or  "12/20 

11  sand"  with  grain  diameters  between  0.59-0.85  mm  and  0.85-1.5  mm,  respectively.  Soltrol  220, 

12  a  mixture  of  alkanes  C13  through  C17  (Phillips  66  Company),  provided  the  oil  phase  for  the 

13  experiments.  Distilled  water  doped  with  potassium  iodine  at  30  g/L,  and  colored  blue  with 

14  FD&C  Blue  provided  the  aqueous  phase.  The  potassium  iodine  increased  the  x-ray  attenua- 

15  tion  of  the  aqueous  phase,  and  the  blue  food  coloring  allowed  a  visual  inspection  of  water 

16  position. 

17 

18  Oil  and  water  specific  tensiometers  were  developed  by  modifying  high-speed  tensiometers 

19  (Selker  et  al.,  1992a)  to  measure  oil  and  water  (^w)  pressures.  The  matric  potential  is  the 

20  difference  of  the  oil  and  water  pressure  measurements.  The  interiors  of  the  oil  tensiometers 

21  were  filled  with  oil  and  a  hydrophobic  plastic  (hydrophilic  glass)  porous  plate  provided  the 

22  connection  to  the  soil.  The  water  tensiometers  used  a  porous  glass  plate  with  relatively  large 
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1  pores.  The  oil  tensiometers  and  their  connections  are  illustrated  in  Figure  2.  A  three-way  valve 

2  connected  pressure  transducers  to  either  a  tensiometer  or  the  atmosphere,  enabling  readings  of 

3  reference  pressures  during  experiments.  Each  transducer  output  was  amplified,  converted  into 

4  a  frequency,  and  captured  by  a  data  acquisition  system.  Calibration  of  the  readings  into 

5  centimeters  of  water  was  achieved  by  recording  the  transducer  output  under  a  sequence  of 

6  measured  water  column  heights.  During  the  experiments,  the  pressure  data  was  recorded 

7  simultaneously  with  the  x-ray  intensity  measurements. 

8 

9  EXPERIMENTAL  PROCEDURE 

10 

11  Initially  the  constants  in  Eqs.  (12)  and  (13)  (x,  tj,  a„,  and  «„),  were  measured.  The  thickness 

12  of  the  chamber,  x,  was  determined  using  a  dial  caliper,  and  the  porosity  was  calculated  from  the 

13  volume  of  oil  required  to  saturate  a  given  volume  of  the  packed  soil  slab. 

14 

15  The  attenuation  constants  and  were  determined  using  the  method  described  by  Lenhard 

16  et  al.  (1988).  A  polycarbonate  standard  cell  with  five  sequential,  0.5  cm  thick  compartments 

17  was  employed.  X-ray  attenuation  was  measured  for  the  (empty)  standard  cell  after  each 

18  compartment  was  sequentially  filled  with  either  water  or  oil.  A  linear  regression  of  the 

19  attenuations  yielded  the  attenuation  constants  for  the  two  fluids.  This  procedure  was  repeated 

20  several  times  throughout  the  experimental  period  to  verify  the  stability  of  the  standards  and  the 

21  beam. 
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1  Six  separate  finger  experiments  were  conducted  using  two  sieved  sand  types  (12/20  and  20/30; 

2  the  numbers  refer  to  the  sieve  numbers)  and  three  water  infiltration  rates  (Table  1).  As  we  were 

3  interested  in  testing  both  the  changing  fluid  contents  versus  time  within  the  finger,  and  the  fluid 

4  contents  in  a  horizontal  cross  section  of  the  finger,  x-ray  attenuation  measurements  were 

5  collected  in  two  different  fashions.  The  first  consisted  of  measuring  the  attenuation  at  one 

6  particular  position  versus  time,  the  second  consisted  of  measuring  the  attenuation  versus  the 

7  horizontal  position  of  the  chamber.  In  a  typical  experiment,  most  of  the  time  was  spent 

8  observing  fluid  changes  within  the  finger  boundaries,  with  periodic  checks  of  how  the  finger 

9  cross  section  had  changed.  The  tensiometer  data  was  taken  at  all  times. 

10 

11  Using  a  peristaltic  pump,  the  sand  was  saturated  from  the  bottom  with  oil.  Once  the  chamber 

12  was  filled  with  oil,  water  was  infiltrated  through  a  point  source  at  the  top  of  the  chamber.  A 

13  closed  drainage  system,  described  below,  maintained  a  relatively  constant  oil  level  and  pressure 

14  even  as  the  water  displaced  the  oil.  A  nalgene  tube  connected  the  bottom  of  the  chamber  to  the 

15  side  of  an  oil  filled  drainage  bottle  (Figure  1).  Another  tube  connected  the  top  of  the  drainage 

16  bottle  to  a  T-fitting  set  at  the  height  of  the  chamber  oil  level.  As  water  displaced  the  oil,  the 

17  excess  oil  drained  through  the  bottom  of  the  chamber  and  out  through  the  T-fitting,  keeping  the 

18  oil  level  constant.  The  experiment  ended  when  water  displaced  all  the  oil  in  the  collecting 

19  bottle,  thus  filling  the  exit  tube  with  water  and  raising  the  oil  pressure. 

20 

21  The  applied  water  initially  spread  both  horizontally  and  vertically  while  at  the  top  of  the 

22  chamber,  but  soon  a  water  finger  formed  and  moved  slowly  down  the  chamber.  In  the  coarse 
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1  12/20  sand  the  finger  moved  straight  down,  while  in  the  finer  20/30  sand  its  path  tended  to 

2  meander  slightly  sideways.  During  re-infiltration  episodes  the  water  followed  its  previous  path 

3  down  through  the  chamber. 

4 

5  Relative  pressure  measurements  in  the  synchrotron  were  found  to  have  an  accuracy  of  0. 1  cm 

6  of  water,  but  absolute  pressures  were  found  to  be  much  less  accurate  (approximately  5  cm  of 

7  water).  In  order  to  improve  the  overall  accuracy,  the  experiments  were  repeated  without  using 

8  the  x-ray  technique  under  the  same  conditions  to  determine  absolute  pressures  in  the  fingers. 

9  These  experiments  showed  identical  relative  pressure  changes,  but  also  provided  absolute 

10  pressures  accurate  to  1.0  cm  of  water.  Thus,  pressure  data  taken  simultaneously  with  the  x-ray 

11  data  was  offset  to  match  the  absolute  pressures  obtained  from  the  concomitant  experiments.  In 

12  these  experiments  we  also  measured  size  and  shape  of  the  fingers. 

13 

14  RESULTS  AND  DISCUSSION 

15 

16  The  observations  of  the  distribution  of  the  fluid  contents  and  pressure  in  the  fingers  is  shown 

17  first  because  it  allows  us  to  illustrate  the  utility  of  the  synchrotron  x-rays  for  documenting  the 

18  relatively  fast  moving  unstable  flows.  Then  we  will  show  how  these  measurements  can  be  used 

19  to  obtain  the  hysteretic  constitutive  relationships  that  explain  why  fingers  are  persistent  and  how 

20  fingers  scale  for  Miller  and  Miller  similar  soils. 

21 


22 
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1  Observations 

2 

3  Table  1  shows  the  diameter  and  moisture  content  in  the  finger  tips  for  the  six  synchrotron 

4  experiments.  The  porosity  of  the  sand  was  0.36.  The  oil  and  water  contents  in  the  tip  of  the 

5  fingers  for  the  12/20  sand  obtained  with  the  synchrotron  experiment  compared  well  with  those 

6  obtained  with  the  light  transmission  technique  of  Damault  et  al.  (1996).  To  our  knowledge  there 

7  are  no  other  measurements  available  for  comparison  purposes  of  fluid  contents  in  rapidly  moving 

8  oil-water  fingers. 

9 

10  Although  in  each  of  the  six  synchrotron  experiments  accurate  measurements  of  the  oil  and  water 

11  contents  in  the  finger  were  obtained,  both  sets  of  tensiometers  were  contacted  by  the  water 

12  finger  in  only  one  of  the  experiments  (#5).  The  12/20  sand  was  used  and  the  water  infiltration 

13  and  (re-infiltration)  rate  was  2.3  ml  min  ^  A  finger,  3-4  cm  in  width,  formed  and  moved  down 

14  the  chamber,  contacting  the  center  sets  of  tensiometers.  Figure  3  shows  the  temporal  changes 

15  in  water  content  and  capillary  potential  (i/'*  -  ^o)  during  the  first  infiltration.  The  passage  of  the 

16  finger  can  be  seen  in  the  time  delays  between  the  response  of  the  top  and  center  tensiometers. 

17  Also,  the  water  content  readings  were  obtained  about  50  sec  later  than  the  center  tensiometer 

18  reading.  The  pressures  and  water  content  exhibit  behavior  typical  of  fingers  and  unstable  flow 

19  as  observed  in  water-air  systems  by  Selker  et  al.  (1992b);  the  water  content  rises  abruptly,  levels 

20  off  temporarily,  and  then  drops  to  a  constant  value  when  the  finger  hits  the  bottom  of  the 

21  chamber  at  around  1000  sec.  When  the  water  was  turned  off  at  approximately  1500  sec,  the 

22  water  content  underwent  another  abrupt  drop  as  the  water  drains  to  the  bottom  and  is  replaced 
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1  by  oil.  The  tensiometer  readings  also  show  the  typical  behavior  of  unstable  flow  in  which  the 

2  matric  potential  decreases  behind  the  front.  Also,  after  the  water  is  turned  off  the  matric 

3  potential  drops  even  further. 

4 

5  Figure  4  shows  the  temporal  changes  in  water  content  during  a  re-infiltration  episode.  In 

6  contrast  to  the  initial  infiltration,  neither  the  water  content  nor  the  pressure  difference  decreased 

7  until  the  water  was  turned  off.  The  flow  within  the  original  finger  boundaries  does  not  show 

8  the  characteristic  drop  of  pressure  behind  the  finger  tip  and  has  characteristics  similar  to  a 

9  "standard"  wetting  front  (i.e.,  when  water  infiltrates  in  a  soil  with  a  horizontal  wetting  front) 

10  in  which  the  matric  potential  increases  at  the  front  and  then  remains  relatively  constant  behind 

11  the  wetting  front.  Also,  the  time  delays  between  the  tensiometer  and  water  readings  is  much 

12  less  than  was  observed  in  the  initial  infiltration  (Figure  3),  showing  a  much  quicker  passage  of 

13  the  pulse.  All  re-infiltration  episodes  exhibited  identical  behavior  although  they  are  not  reported 

14  here. 

15 

16  Figure  5  shows  the  horizontal  profile  at  selected  times  during  the  infiltration  and  re-infiltration 

17  cycles.  The  profile  extends  from  about  the  middle  of  the  finger  to  outside  of  the  finger, 

18  depicting  the  wide  variation  of  water  contents  at  the  same  height.  The  measured  horizontal 

19  boundary  of  the  finger  is  about  8-10  mm  wide,  but  the  actual  boundary  is  likely  much  smaller 

20  than  this,  as  much  of  the  measured  width  is  due  to  the  fact  that  the  x-ray  beam  was  6  mm  wide. 

21  In  addition,  the  boundary  (and,  thus,  the  width  of  the  finger)  did  not  change  measurably  during 


18 


1  re-infiltration  episodes.  Note  that  this  follows  the  air-water  finger  observations  in  which  during 

2  re-infiltration  the  finger  width  remains  the  same  and  the  moisture  content  is  lower. 

3 

4  Persistence 

5 

6  Most  theory  of  unstable  flow  has  been  focussed  on  the  instability  itself.  Persistence  of  fingers 

7  over  time  especially  in  oil-water  systems  has  received  only  little  attention  as  noted  by  Glass  and 

8  Nicholl  (1996).  In  air- water  systems  the  hysteretic  nature  of  the  constitutive  relationship  has 

9  been  crucial  to  understanding  the  finger  persistence  (Glass  et  al.,  1989c;  Liu  et  al.,  1994b).  We 

10  will  show  here  that  for  oil-water  systems  the  hysteresis  plays  much  the  same  role. 

11 

12  The  synchrotron  data  presented  above  allows  us  to  construct  for  oil- water  systems  the  constitu- 

13  tive  relationship  between  the  fluid  pressure  difference  and  water  content.  Taking  into  account 

14  the  travel  time  between  pressure  and  fluid  content  measurement  location,  the  data  in  Figures  3 

15  and  4  are  replotted  to  give  the  primary  and  secondary  wetting  loops  in  Figure  6.  In  this  figure, 

16  the  letter  "A"  refers  to  the  primary  wetting  branch  beginning  at  Al  (oil  saturated  soil)  via  A2 

17  (finger  front  starts  affecting  the  readings)  to  A3  (just  behind  the  finger  tip  front).  Shortly  after 

18  the  finger  tip  passes  the  tensiometers,  the  pressure  starts  decreasing  at  B1  which  is  the  starting 

19  point  of  the  primary  drying  branch,  "B".  Between  points  B1  and  B2,  the  matric  potential 

20  decreases  further  while  the  finger  remains  at  the  high  moisture  content.  This  represents  the 

21  wetted  water  bulb  similar  to  the  air-water  systems  (Selker  et  al.,  1992b).  Then  while  the  finger 

22  tip  is  moving  farther  down,  both  moisture  content  and  matric  potential  are  decreasing  until  point 
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1  B3  where  the  water  is  turned  off.  At  point  B4  no  water  is  draining  out  of  the  chamber  anymore. 

2  Branches  "C"  and  "D"  are  the  wetting  and  drying  curves  during  re-infiltration.  Both  water 

3  content  and  matric  potential  are  lower  in  the  finger  tip  during  re-infiltration  than  the  initial 

4  infiltration  resulting  in  secondary  wetting  and  drying  branches  within  the  primary  loops.  All  the 

5  re-infiltrations  followed  the  second  smaller  hysteresis  loop  (only  one  is  shown  here). 

6 

7  The  data  in  Figure  6  make  it  possible  to  understand  why  the  finger  during  the  initial  infiltration 

8  into  the  oil  saturated  soil  does  expand  only  in  the  lowest  portion  of  the  tip  and  why  the  finger 

9  during  re-infiltration  does  not  go  outside  its  confines  set  during  the  initial  experiments. 

10  Sideways  movement  of  the  finger  is  governed  by  the  sorptivity.  Only  when  the  sorptivity  is  zero 

11  the  finger  will  not  expand  sideways.  The  sorptivity,  S,  can  be  defined  as  (Parlange,  1975; 

12  Brutsaert,  1976): 

♦/ 

//§;  ir,(i|>)  dill  (W) 

13  where  the  conductivity,  K,v(^)  and  the  moisture  content,  are  functions  of  matric  potential  on 

14  the  wetting  branch.  The  integration  constants  are  ypf  which  is  the  matric  potential  inside  the 

15  finger  and  is  potential  outside  the  finger.  Although  the  conductivity  function  is  not  known 

16  it  is  obvious  that  the  is  zero  when  the  water  content  is  zero.  Consequently,  if  the  moisture 

17  content  of  the  wetting  curve  between  the  integration  constants  is  zero,  we  can  safely  assume  that 

18  the  sorptivity  is  zero  too. 

19 
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1  To  examine  whether  the  finger  is  expanding  it  is  necessary  to  investigate  the  relationships  of 

2  fluid  contents  in  and  outside  the  finger.  Earlier  we  found  that  the  water  content  matric  potential 

3  relationship  for  the  soil  outside  the  finger  is  on  the  primary  wetting  loop  (indicated  by  "A”)  and 

4  the  soil  inside  the  finger  is  on  curve  "B"  during  the  initial  infiltration  experiment  or  curves  "C" 

5  or  "D"  for  each  of  the  re-infiltration  experiments  (Figure  6). 

6 

7  The  matric  potential  decreases  with  increasing  distance  above  the  finger  tip.  When  the  matric 

8  potential,  \J/(,  decreases  below  the  matric  potential  of  point  A2  (-0.75  cm  in  Figure  6),  the 

9  moisture  content  between  and  is  zero  and,  hence,  the  sorptivity  calculated  in  Eq.  (8) 

10  becomes  zero.  Accordingly,  anywhere  in  the  finger  with  a  matric  potential  less  than  -0.75  cm, 

11  the  finger  will  not  expand.  Only  in  the  region  near  the  finger  front  the  finger  is  expanding. 

12 

13  During  re-infiltration  the  matric  potential  inside  the  finger  varies  from  -2  to  -4  cm.  At  these 

14  tensions  the  water  contents  on  the  primary  wetting  curve  remain  zero.  Thus,  the  conductivity 

15  and  sorptivity  are  also  zero,  and  the  finger  cannot  cross  the  original  boundary  of  the  earlier 

16  established  finger.  We  can  now  also  understand  why  the  finger  has  a  stable-like  v/etting  front 

17  inside  the  "old"  (i.e.,  previously  formed)  finger.  Within  the  "old"  finger  the  infiltration  and 

1 8  sideways  movement  are  governed  by  the  secondary  hysteresis  loops,  with  a  finite  sorptivity  value 

19  at  all  pressures  behind  the  finger  tip. 

20 

21  These  findings  are  important  because  it  shows  that  knowledge  concerning  the  constitutive 

22  relationships  are  required  for  hydrodynamic  control  in  remediation  sites  with  immiscible  liquids. 
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1  It  also  shows  that  when  models  do  not  include  these  hysteretic  constitutive  relationships,  no 

2  realistic  planning  of  cleanup  procedures  can  be  achieved.  The  difficulty  is,  however,  that  the 

3  hysteretic  constitutive  relations  for  oil-water  systems  are  generally  not  known.  We  will 

4  investigate  next  if  Miller  and  Miller  (1956)  scaling  can  be  used  to  obtain  the  hysteretic 

5  constitutive  relations  as  well  as  the  finger  diameter. 

6 

7  Scaling  of  Fingers 

8 

9  Before  we  compare  the  observed  and  predicted  width  by  scaling  we  will  justify  the  scaling 

10  parameters  used. 

11 

12  Glass  et  al.  (1991a),  Parlange  et  al.  (1990),  and  Schroth  et  al.  (1996)  found  that  the  sands  in 

13  Table  2  are  Miller  and  Miller  scale  similar  and  can  be  scaled  according  to  Eq.  (8)  with  the  mean 

14  particle  sizes  given  in  Table  2. 

15 

16  The  scaling  factor  for  interfacial  tension,  f„,  is  somewhat  uncertain  due  to  unknown  contact  angle 

17  (Chandler  et  al.,  1996)  and  the  change  of  interfacial  tension  of  oil  and  water  in  time  (Schroth 

18  et  al.,  1995).  Here  we  derive  the  scaling  for  interfacial  tension,  f„,  by  matching  the  soil 

19  characteristic  curves  of  oil-water  and  air-water  systems.  This  has  the  advantage  that 

20  includes  the  unknown  contact  angle.  Schroth  et  al.  (1996)  found  f,,  values  between  1 .6  and  1 .8 

21  by  matching  drainage  curves  in  oil-water  and  air-water  systems  for  different  sand  sizes.  In 

22  Figure  7,  we  compared  the  scaled  12/20  air-water  drainage  curve  of  Schroth  et  al.  (1996)  using 
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1  f„  =  1.76  and  the  soil  oil-water  characteristic  curve  (Figure  5).  The  scaled  soil  characteristic 

2  air-water  curve  for  20/30  sand  (Liu  et  al. ,  1994a)  is  also  depicted.  All  three  drainage  curves 

3  are  close  and  certainly  within  the  experimental  error  of  the  tensiometers.  The  slight  difference 

4  in  slope  between  Schroth  et  al.  (1996)  and  both  Liu  et  al.  (1994a)  and  ours  is  likely  caused  by 

5  the  calculation  procedure  followed  by  Schroth  et  al.  (1996)  to  correct  for  the  difference  in 

6  moisture  content  within  the  tempi  cell.  The  wetting  curves  of  Liu  et  al.  (1994a)  and  ours  are 

7  also  very  close.  The  slight  offset  between  the  two  curves  is  either  caused  by  the  oil-water  curve 

8  not  going  all  the  way  to  saturation  or  slight  differences  in  reading  the  absolute  pressures.  A 

9  value  of  f„  =  1.76  or  an  interfacial  tension  of  41  g*sec'^  for  oil-water  and  72  g‘sec'^  for  air- 

10  water  will  be  used  here. 

11 

12  The  scaling  factor  for  2-d  and  3-d  fingers,  fn,  for  air-water  systems  is  equal  to  4.8/x  (Glass  et 

13  al.,  1990).  It  is  the  same  for  oil-water  systems.  The  scaling  factor,  fo,  is  either  equal  to  1  or 

14  0.2  depending  on  the  reference  fluid. 

15 

16  Table  2  gives  the  finger  diameters  for  a  number  of  experiments  reported  in  the  literature.  For 

17  2-d  fingers  the  average  finger  width  was  taken.  For  3-d  column  experiments  by  Yao  and 

18  Hendrickx  (1996)  coalescence  of  fingers  is  a  problem  in  measuring  the  diameter  and  because 

19  most  fingers  were  small  (see  their  Figure  6).  The  minimum  finger  diameter  was  more 

20  appropriate  and  is  reported  in  Table  2. 
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1  Using  the  above  scaling  factors  and  the  20/30  sand  oil-water  as  the  reference  soil  and  liquid  with 

2  an  average  diameter  of  6  cm,  the  finger  diameters  for  the  other  experiment  in  Table  2  were 

3  calculated  with  Eq.  (7)  for  b  =  1  and  b  =  0.5.  The  predicted  scaled  finger  diameters  are 

4  compared  with  the  observed  finger  diameters  in  Figure  8.  Note  that  not  all  moisture  contents 

5  in  the  tips  were  measured  and  some  were  estimated  based  on  similar  experiments.  In  general, 

6  there  was  a  satisfactory  agreement  between  the  scaled  finger  size  and  the  observed  value.  A 

7  value  of  b  =  1  gives  a  slightly  better  estimate  than  b  =  0.5. 

8 

9  CONCLUSIONS 

10 

11  Unstable  flow  experiments  were  carried  out  in  slab  chambers  with  Soltrol  and  water.  High- 

12  intensity  synchrotron  radiation  coupled  with  pressure  measurements  was  an  effective  method  of 

13  studying  flow  instabilities  in  oil-saturated  porous  media. 

14 

15  Hysteresis  in  the  constitutive  relations  was  the  governing  mechanism  for  the  temporal  persistence 

16  of  fingered  flow.  Hysteresis  provided  the  soil  with  a  memory  of  previously  formed  finger  paths. 

17 

18  The  results  show  that  fingered  flow  into  an  oil  saturated  porous  media  has  similar  properties  to 

19  fingered  flow  into  dry  porous  media,  but  with  different  length  scales,  pressure  heads,  and  fluid 

20  contents.  Scaling  of  unstable  flow  paths  might  overcome  some  of  the  difficulties  of  measuring 

21  the  constitutive  relations  in  oil- water  systems  for  different  soil. 
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1  These  results  will  ultimately  be  incorporated  in  theories  of  preferential  flow  and  lead  to  better 

2  hydrological  control  of  in  situ  remediation  in  which  immiscible  fluids  are  a  problem. 
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1  Table  1:  Synchrotron  experiments. 
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1  Table  2:  Overview  of  finger  experiments  and  parameters  needed  for  scaling. 
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DiCarlo  et  al..  1996 

37 

72 

1.1 
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1  nOURE  CAPTIONS 

2 

3  Figure  1:  Front  and  side  view  of  the  experimental  chamber.  There  was  one  set  of  top 

4  tensiometers  and  three  sets  of  middle  tensiometers.  Water  contents  were  measured 

5  at  the  scanning  elevation,  directly  below  the  middle  tensiometers.  The  bottle  and  the 

6  outflow  tube  were  used  to  keep  the  oil  pressure  constant  during  water  infiltration. 

7 

8  Figure  2:  Close-up  view  of  an  oil  tensiometer  inserted  into  the  chamber.  Water  tensiometers 

9  were  completely  filled  with  water  and  used  a  ceramic  porous  plate. 

10 

11  Figure  3:  Water  content  and  pressure  measurements  as  the  initial  finger  moved  down  the 

12  chamber  through  oil-saturated  12/20  sand.  Both  water  content  and  pressure  rose 

13  abruptly  at  the  front  of  the  finger,  and  fell  behind  the  finger  tip,  indicative  of 

14  unstable  flow.  The  water  was  turned  on  at  t  =  0  sec.  When  the  water  was  turned 

15  off  at  t  =  1500  sec,  both  readings  dropped  again. 

16 

17  Figure  4: 

18 

19 

20 

21  Figure  5:  Horizontal  profile  of  the  water  content  within  the  finger.  There  is  no  increase  in 

22  finger  size  for  re-infiltration  events. 


Water  content  and  pressure  measurements  during  a  re-infiltration  of  the  finger.  No 
drops  in  content  or  pressure  are  observed  while  the  water  is  on,  indicating  that  the 
re-infiltration  is  now  in  a  stable  flow  pattern. 
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Figure  6:  Pressure-water  content  relationship  from  the  finger  in  the  12/20  sand.  A  and  B  are 
the  initial  wetting  and  drying  curves,  and  C  and  D  are  the  secondary  wetting  and 
drying  curves. 

Figure  7:  Observed  (obs)  and  scaled  (sc’d)  oil-water  characteristic  curves  for  Soltrol  220  and 
water  (o/w).  The  scaled  curves  were  derived  from  Schroth  et  al.  (1996)  and  Liu  et 
al.  (1994a). 

Figure  8:  Comparison  of  scaled  finger  diameter  (Eqs.  (7-11))  according  to  theory  originally 
developed  by  Parlange  and  Hill  (1976)  (b  =  1)  and  Chuoke  et  al.  (1959)  (b  =  0.5). 
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1  ABSTRACT 

2 

3  The  difficulty  of  determining  transient  fluid  contents  in  a  soil-oil-water  system  is  hampering  an 

4  understanding  of  the  system’s  flow  characteristics.  In  this  paper,  we  describe  a  light 

5  transmission  method  (LTM)  which  can  rapidly  obtain  oil  and  water  contents  throughout  a  large 

6  two-dimensional  flow  field.  By  appropriately  coloring  the  water  with  0.005%  FD&C  blue  #1, 

7  the  hue  of  the  transmitted  light  is  found  to  be  directly  related  to  the  water  content  within  the 

8  porous  media.  Then,  the  hue  provides  a  high  resolution  measurement  of  the  water  and  oil 

9  contents  in  transient  flow  fields  (such  as  unstable  flow).  Evaluation  of  the  reliability  of  LTM 

10  was  assessed  by  checking  mass  balance  for  a  known  water  injection,  and  its  utility  in  visualizing 

11  a  whole  flow  field  was  exemplified  for  an  unstable  fingered  flow  by  comparing  fluid  contents 

12  to  those  obtained  with  synchrotron  x-ray  radiation. 


3 


1  INTRODUCTION 

2 

3  The  complexity  of  simultaneous  flow  of  water  and  non-aqueous  phase  liquids  (NAPLs)  is  little 

4  understood  because  few  techniques  permit  accurate  measurement  of  water  and  oil  contents  in 

5  transient  flow  fields.  Mercer  and  Cohen  (1990)  noted  that  in  situ  fluid  content  measurements 

6  are  needed  to  understand  immiscible  fluid  flows  in  porous  media  and  to  perform  effective 

7  subsurface  remediation. 

8 

9  Most  techniques  available  for  visualization  of  oil  and  water  fluid  contents  can  measure  only  near 

10  steady  state  flow  conditions.  They  include  conventional  x-ray  and  gamma  ray  radiation 

11  techniques  used  by,  among  others,  Hopmans  and  Dane  (1986),  Le.nhr.rd  et  al.  (1988),  Dane  et 

12  al.  (1992),  and  Illangasekare  et  al.  (1995).  These  measurements  are  slow  because  of  the 

13  relatively  long  counting  times  and  only  one  point  can  be  measured  at  one  time.  Synchrotron  x- 

14  rays  (Liu  et  al.,  1993;  Rimmer  et  al.,  1996)  allow  shorter  counting  times  (in  the  order  of 

15  seconds)  but  only  regions  of  less  than  0.5  cm^  can  be  measured  at  one  time.  Tidwell  and  Glass 

16  (1994)  have  performed  full  field  measurements  using  x-ray  film,  but  recording  times  are  too  long 

17  for  measuring  transient  flow  phenomena. 

18 

19  Transient  visualizations  (but  not  a  direct  measurement  of  fluid  contents)  have  been  made  in  Hele- 

20  Shaw  cells  with  smooth  walls  (Saffman  and  Taylor,  1958;  Chuoke  et  al.,  1959)  or  with  imprints 

21  of  porous  media  on  glass  (Schwille,  1988).  The  visualization  method  of  Van  Geel  and  Sykes 

22  (1994)  allows  measurement  of  liquid  content  near  the  wall  of  the  slab  chambers  by  relating  the 
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1  gray  level  of  the  reflected  light  to  water  content.  However,  as  they  noted,  the  fluid  content  near 

2  the  wall  might  not  necessarily  be  the  same  as  within  the  porous  media. 

3 

4  Thus,  there  is  a  need  for  a  method  that  allows  full  field  moisture  content  visualization  in  two- 

5  phase  liquid  systems.  This  paper  describes  the  development  of  such  a  technique  by  using  light 

6  transmission  through  a  two-dimensional  slab  chamber.  By  coloring  the  water  phase  with  0.005  % 

I  FD&C  blue  #1  the  hue  content  of  the  transmitted  light  is  directly  related  to  the  water  content. 

8  The  technique  is  a  new  application  of  the  method  of  Glass  et  al.  (1989)  for  air-water  systems 

9  in  which  the  light  intensity  was  directly  related  to  water  content.  Unlike  air  and  water,  in  the 

10  case  of  oil-water  systems  the  refractive  index  of  the  two  fluids  is  almost  identical,  making  it 

I I  impossible  to  measure  liquid  contents  of  each  of  the  fluids  from  the  light  intensity  measurements. 

12  A  blue  dye  was  chosen  because  it  produced  a  wide  color  spectrum  going  from  yellow  for  the 

13  oil  saturated  sand  via  green  to  blue  for  water  saturated  soil. 

14 

15  VISUALIZATION 

16 

17  Quantification  of  color  differences  can  be  complicated,  because  there  are  various  ways  to  specify 

18  "color".  The  color  for  video  cameras,  color  monitors  and  computer  graphics  is  defined  in  terms 

19  of  a  vector  with  the  three  components  of  red,  green,  and  blue  intensities  (RGB).  Another  system 

20  to  specify  the  color  vector  involves  hue,  saturation,  and  intensity  (HSI),  which  is  used,  for 

21  instance,  by  the  Munsell  color  chart  for  soil  classification.  Hue  is  the  attribute  that  describes 

22  the  pure  color  and  is  what  we  are  typically  referring  to  when  we  use  the  term  "color". 
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Saturation  is  the  attribute  that  describes  the  degree  to  which  the  color  is  diluted  with  white. 
Intensity  is  the  attribute  which  corresponds  to  the  gray  level  (black  and  white)  of  the  color 
image.  The  advantages  of  the  HSI  format  is  that  it  treats  color  roughly  the  same  way  that 
humans  perceive  and  interpret  color.  Therefore,  because  color  differences  are  apparent  to  the 
human  eye  for  different  water  and  oil  contents,  we  expect  that  water  content  is  directly  related 
with  hue  rather  than  intensity  or  saturation. 


Hue,  H,  in  HSI  format  can  be  obtained  from  the  RGB  vector  (R=red,  G=green,  and  B=blue 
intensity,  respectively,  with  values  for  each  vector  component  ranging  from  0  to  255)  following 
(Wilson,  1988): 


H  =  255 


1 

360 


Y-ARCTAN~-~^~^ 

y/3(G-B) 


(1) 


Y  =  90  for  G>B, 

Y  =  270  for  G<B. 

where  Y  is  a  constant  depending  on  the  green  and  blue  intensity.  This  equation  shows  that  if, 
as  we  expect,  H  is  the  best  variable  to  measure  water  content,  each  of  the  components  of  the 
RGB  vector  is  affected  by  the  water  content  and  the  combination  (2R-G-B)/(G-B)  would  be  the 
relevant  parameter  in  the  RGB  system  to  relate  to  water  content. 


MATERIALS  AND  METHODS 


For  fluid  measurements  in  porous  media,  light  transmission  techniques  involve  placing  a  two- 
dimensional  experimental  chamber  in  front  of  an  uniform  light  source  and  recording  the 
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1  transmitted  light  (Glass  and  Nicholl,  1995).  We  used  a  light  source  composed  of  a  bank  of  24 

2  florescent  light  bulbs  in  front  of  a  white  background.  Experimental  chambers  had  1  cm  thick 

3  polycarbonate  walls  with  a  1  cm  thick  interior.  The  transmitted  light  was  recorded  with  a  Cohu 

4  Solid  State  Color  Camera  located  1  m  in  front  of  the  chamber,  with  constant  settings  (zoom  = 

5  3  -  5  m  and  aperture  =  f  1.8)  and  stored  on  standard  VHS  tape  in  RGB  format.  Recorded 

6  images  were  converted  from  RGB  to  HSI  format  with  an  IBM  compatible  computer  equipped 

7  with  a  DT2871  Color  Frame  Grabber  and  analyzed  with  a  Color  Tutorial  (SP0228)  VI. (X) 

8  software,  both  from  Data  Translation,  Inc.,  -  07839A  1988. 

9 

10  To  obtain  a  calibration  curve  between  hue  value  and  oil-water  content,  a  two-dimensional 

1 1  calibration  chamber,  consisting  of  compartments  with  known  quantities  of  oil  and  water,  was 

12  constructed  (Figure  1).  The  oil  used  was  Soltrol  220,  an  isoparaffine  solvent  composed  of  a 

13  mixture  of  alkanes  C13  through  C17  from  the  Phillips  66  Co.  The  water  was  colored  with 

14  FD&C  blue  #1,  at  a  concentration  of  0.005%  by  weight.  By  trial  and  error,  the  color  blue  and 

15  the  specific  concentration  of  0.005%  were  found  to  give  the  widest  range  of  hue  values.  The 

16  calibration  chamber  had  dimensions  of  32.5  cm  high,  26  cm  wide,  and  1  cm  thick,  and  was 

17  divided  into  six  cells  of  3  cm  high,  26  cm  wide,  and  1  cm  thick.  The  walls  consisted  of  a  1  cm 

18  thick  Hyzod  polycarbonate  sheet,  made  by  Sheffield  Plastic,  Inc.  Cell  walls  were  covered  with 

19  black  tape  to  avoid  edge  effects  of  transmitted  light  through  the  plastic  walls.  The  cells  were 

20  packed  with  12/20  sieve  size  sand  and  water  contents  of  0,  20,  40,  60,  80,  and  100%  saturation 

21  (volumetric  moisture  content  of  0.(X),  0.07,  0.14,  0.21,  0.28,  and  0.35  cmVcm^,  respectively). 
The  remainder  of  the  pore  space  was  filled  with  oil.  The  coarse  12/20  sand  was  industrial 
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1  quartz  (Unimin  Corp.)  with  particle  sizes  between  0,85  and  1.5  mm,  which  was  packed  to  a 

2  porosity  of  35%.  The  0  and  100%  water  saturations  were  obtained  by  dripping  oil  and  water 

3  into  the  sand- filled  cells,  respectively.  The  20,  40,  60,  and  80%  water  saturations  were  made 

4  by  mixing  the  required  water  and  oil  quantities  in  a  beaker  with  the  sand,  stirring  the  mixture, 

5  loading  the  mixture  in  the  cells,  and  then  packing  the  cells  by  vibration. 

6 

7  For  each  cell,  the  hue  was  taken  as  the  mean  values  present  in  a  centered  line  of  100  pixels.  A 

8  histogram  of  all  the  hue  values  located  on  the  line  was  made  for  each  cell.  The  mean  hue  values 

9  were  plotted  versus  the  corresponding  water  content  to  obtain  a  calibration  curve. 

10 

11  Two  types  of  experiments  were  conducted  to  illustrate  the  reliability  and  applicability  of  the 

12  method.  The  first  one  aimed  to  test  the  accuracy  of  the  technique  for  a  simple  flow  experiment 

13  from  a  point  source  to  check  the  mass  balance,  while  the  second  one  illustrated  the  technique 

14  for  a  more  complex  flow  produced  by  fingering,  similar  to  the  experiment  of  Rimmer  et  al. 

15  (1996)  with  synchrotron  x-rays.  These  experiments  were  carried  out  in  a  two-dimensional  slab 

16  chamber  (Figure  2),  The  interior  dimensions  of  the  experimental  chamber  were  57  cm  high, 

17  51  cm  wide,  and  1  cm  thick.  Two  metal  bars  were  fixed  at  the  top  and  middle  of  the  chamber 

18  to  prevent  it  from  expanding  while  filling  with  oil.  The  walls  consisted  of  1  cm  thick  Hyzod 

19  polycarbonate  sheets  made  by  Sheffield  Plastic,  Inc,  The  experimental  chamber  was 

20  continuously  filled  with  industrial  quartz  12/20  sieve  size  sand  through  a  funnel-extension- 

21  randomizer  assembly  to  minimize  segregation  and  heterogeneity  (Glass  et  al.,  1989).  Extra  sand 

22  was  added  to  have  a  10  cm  sand  layer  in  the  funnel  to  achieve  a  natural  packing  in  order  to 
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1  avoid  the  movement  of  sand  particles  during  the  experiment.  After  packing,  the  top  2  cm  of  the 

2  sand  was  removed  by  sweeping  to  bring  the  sand  height  to  55  cm.  The  experimental  chamber 

3  was  saturated  with  oil  to  a  height  of  55  cm  by  adding  Soltrol  from  the  bottom  at  a  constant  rate 

4  of  8.4  ml/min.  The  oil  level  was  kept  constant  and  excess  oil  was  drained  through  the  bottom 

5  of  the  chamber  via  an  overflow  at  the  same  level  as  the  Soltrol  upper  level  in  the  chamber.  The 

6  experiments  were  performed  at  a  constant  temperature  of  20'’C. 

7 

8  The  first  test  consisted  of  checking  the  mass  balance.  Water  was  injected  with  a  flow  rate  of 

9  2.8  ml/min  at  a  point  source  located  3  cm  below  the  oil-saturated  12/20  sand  surface  in  the 

10  experimental  chamber.  The  hue  values  were  taken  for  all  pixels  constituting  the  area 

11  surrounding  the  finger  and  converted  to  water  content. 

12 

13  For  the  second  test,  we  measured  fluid  contents  for  unstable  flow  of  water  into  oil-saturated 

14  sand.  A  set  of  three  experiments  were  performed  in  the  experimental  chamber.  The  water  was 

15  uniformly  applied  through  a  needle,  attached  to  a  cam  driver  which  went  back  and  forth  at  a 

16  repetition  rate  of  8  seconds,  along  the  top,  1  cm  above  the  oil-saturated  12/20  sand  surface. 

17  Figure  2  depicts  the  experimental  setup.  Water  fingers  were  observed  under  three  different 

18  infiltration  flow  rates:  4.5,  7.5,  and  10.4  ml/min. 

19 


20 
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1  RESULTS 

2 

3  Calibration  of  the  Light  Transmission  Method  (LTM) 

4 

5  Color  differences  were  visible  between  the  different  cells  of  the  calibration  chamber:  the  oil- 

6  saturated  cell  appeared  yellow,  the  water-saturated  cell  appeared  blue  and  the  oil-water  saturated 

7  cells  appeared  varying  shades  of  green.  In  each  cell,  the  color  was  not  completely  uniform, 

8  slight  spatial  differences  existed  due  to  oil  and  water  ganglia  formation  (Figure  1). 

9 

10  Inside  each  cell,  the  color  attributes  of  pixels  were  correlated  with  water  content  in  HSI  and 

11  RGB  formats.  Using  RGB  format,  even  though  the  color  went  from  yellow  to  green  to  blue, 

12  only  the  red  component  of  the  RGB  vector  showed  any  trend  with  water  content.  This  trend, 

13  however,  was  not  sufficient  to  give  an  unique  relationship  with  water  content.  The  difficulty 

14  in  using  RGB  is  the  interdependence  of  color  saturation  with  values  of  the  components  of  the 

15  color  vector  (Wilson,  1988).  In  our  case,  the  bright  yellow  (0%  water)  had  approximately  the 

16  same  value  for  the  blue  component  as  the  100%  water  which  was  pale  blue.  Thus,  even  if  the 

17  eye  can  see  a  difference  in  color,  RGB  is  unable  to  pick  out  the  colors  in  a  simple  predictable 

18  manner.  In  HSI  format,  the  hue  (representing  the  color  for  the  human  eye)  is  independent  of 

19  the  color  saturation  (Wilson,  1988)  and,  indeed,  the  histograms  of  hue  values  for  the  water 

20  content  in  each  cell  were  found  to  be  unique,  whereas  neither  saturation  nor  intensity  had 

2 1  significant  correlation  to  water  content.  The  distribution  of  hue  values  for  the  six  different  water 

22  contents  is  presented  in  Figure  3.  The  yellow  oil-saturated  cell  had  a  hue  value  of  about  43,  and 
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1  the  blue-green  water-saturated  cell  had  a  hue  value  of  about  86.  The  oil-saturated  and  water- 

2  saturated  cells  had  the  most  uniform  color,  and  the  narrowest  hue  histograms.  Only  two  of  the 

3  histograms  had  overlapping  hue  values  because  of  formation  of  some  oil  and  water  ganglia,  but 

4  there  was  a  difference  in  peak  positions.  The  hue  value  increased  with  the  water  content  and 

5  a  linear  regression  provided  a  good  fit  with  a  of  97.8  %  (Figure  4).  The  hue-water  content 

6  equation  is; 

0  =  -0.9395  +  0.0222/f  (2) 

7  where  6  is  the  volumetric  water  content  and  H  is  the  hue  value  calculated  with  Eqn.  (1)  from 

8  the  recorded  RGB  color  attributes.  Note  that  the  constants  in  Eqn.  (2)  are  specific  for  our 

9  experimental  setup.  The  system  needs  to  be  calibrated  when  different  video  cameras  or  particle 

10  sizes  are  used. 

11 

12  Accuracy,  Utility,  and  Reliability  of  LTM 

13 

14  The  LTM  was  evaluated  with  a  mass  water  balance  and  synchrotron  x-ray  technique.  The  mass 

15  balance  test  consisted  of  measuring  a  known  water  volume  in  oil-saturated  sand.  The  injected 

16  20  ml  of  water  generated  a  finger  of  dimensions  approximately  10  cm  long,  6  cm  wide,  and  1 

17  cm  thick  in  the  experimental  chamber.  Hue  values  were  obtained  by  analyzing  all  the  14400 

18  pixels  present  in  a  surrounding  area  that  contains  the  finger,  and  were  converted  to  water  content 

19  using  Eqn.  (2).  The  total  amount  of  water  calculated  was  19.85  ml,  which  was  within  1%  of 

20  the  amount  of  water  applied. 
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1  As  an  application  and  to  check  the  accuracy  of  LTM  under  transient  flow  conditions,  water 

2  infiltrations  into  12/20  oil-saturated  sand  were  performed  for  three  flow  rates.  An  unstable 

3  wetting  front  formed  with  one  or  more  fingers.  An  example  of  the  fingered  flow  of  water  for 

4  the  flow  rate  of  4.5  ml/min  is  presented  in  Figure  5.  Before  the  first  finger  went  down,  a 

5  horizontal  wetting  front  of  about  2  cm  deep  formed.  The  finger  was  the  bluest  (highest  value 

6  for  hue)  at  the  tip,  indicating  the  highest  water  content.  Figure  6  shows  the  finger  water  content 

7  profile  after  15  and  20  minutes.  LTM  measurement  indicated  that  after  15  minutes  the  water 

8  content  at  the  finger  tip  was  between  0.27  and  0.30  cm^/cm\  and  the  finger  width  ranged  from 

9  6  to  8  cm.  At  20  minutes,  the  finger  tip  water  content  was  between  0.26  and  0.28  cm^/cm^. 

10  Approximately  25  cm  above  the  finger  tip,  the  moisture  content  decreased  to  an  average  of  0. 17 

11  cm^/cm^.  The  horizontal  water  content  profile  at  a  depth  of  48  cm  at  20  minutes  through  the 

12  tip  of  the  finger  is  shown  in  Figure  7.  Both  infiltration  flow  rates  of  7.5  and  10.5  ml/min 

13  produced  a  wetting  front  instability  that  generated  2  and  3  fingers,  respectively,  with 

14  approximately  the  same  characteristics  as  the  one  generated  with  4.5  ml/min.  The  only 

15  difference  was  that  the  water  content  was  slightly  higher  and  the  velocities  sometimes  slower 

16  (Table  1). 

17 

18  The  above  results  were  compared  with  available  data  from  synchrotron  x-ray  measurements  by 

19  Rimmer  et  al.  (1996)  with  12/20  oil-saturated  sand  and  a  flow  rate  of  2.3  ml/min  through  a  point 

20  source  (Figure  6).  For  the  synchrotron  measurements,  unlike  LTM,  there  is  a  fundamental 

21  relationship  between  the  attenuation  of  x-rays  and  moisture  content.  Although  slight  differences 

22  existed  in  observed  water  contents,  which  have  been  caused  by  slight  different  experimental  flow 
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1  rates,  both  methods  follow  the  same  finger  behavior.  The  highest  water  content  was  measured 

2  at  the  finger  tips:  0.265  cm^/cm^  for  the  x-ray  radiation  and  0.26  to  0.28  cm^/cm^  for  the  LTM. 

3  Also,  in  the  transition  zone  behind  the  finger  tip,  the  decrease  in  water  content  is  similar  over 

4  the  zone  that  the  synchrotron  was  available.  In  Figure  7,  the  moisture  content  in  a  horizontal 

5  profile  at  comparable  locations  is  also  very  similar  between  the  two  methods. 

6 

7  DISCUSSION 

8 

9  The  Light  Transmission  Method  (LTM)  was  found  to  yield  reliable  fluid  contents  in  transient 

10  oil-water  flow  fields  as  checked  by  the  mass  conservation  experiment.  Unstable  fingered  flow 

11  phenomena  measured  with  LTM  gave  similar  results  as  those  with  the  synchrotron  x-ray  source. 

12  The  advantage  of  the  LTM  is  that  it  records  the  fluid  contents  of  the  whole  flow  in  less  than 

13  0.05  seconds  (the  time  it  takes  for  a  video  camera  to  take  an  image).  Synchrotron  x-ray 

14  measurements  can  be  taken  also  within  one  second  but  only  at  one  location  at  a  time.  The 

15  disadvantage  of  LTM  is  that  silica  sand  has  to  be  used  as  porous  media  while  synchrotron  x-ray 

16  measurements  are  independent  of  the  translucence  of  the  porous  media.  The  LTM  when  coupled 

17  with  simultaneous  pressure  measurements  will  be  useful  in  validating  1-D  and  2-D  computer 

18  codes  for  transient  oil  and  water  flow. 
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1  CONCLUSIONS 

2 

3  A  method  was  developed  that  can  measure  full  field  average  fluid  contents  for  transient  flows 

4  in  soil-oil-water  systems.  Currently,  there  is  a  lack  of  such  measuring  techniques,  hampering 

5  the  development  of  theory  and  testing  of  two-phase  models. 

6 

7  We  found  that  by  dying  the  water  light  blue  and  analyzing  the  flow  field  for  this  blue  color  in 

8  HSI  format,  we  could  determine  the  water  content  in  the  slab  chamber  filled  with  sand.  RGB, 

9  which  is  used  in  recording  color  with  video  cameras,  could  not  easily  quantify  the  obvious  color 

10  difference  between  yellow  and  blue.  By  using  the  HSI  format,  an  unique  relationship  existed 

11  between  the  hue  and  water  content.  Comparison  with  moisture  contents,  determined  with 

12  synchrotron  x-rays,  showed  that  the  LTM  can  be  successfully  used  for  full  field  moisture 

13  contents  in  porous  media  consisting  of  silica  sand. 

14 
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High-speed  measurements  of  three-phase  flow 
using  synchrotron  X  rays 
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Abstract.  Accurate  monitoring  of  flow  instabilities,  which  can  occur  when  nonaqueous 
phase  liquids  (NAPLs)  flow  through  porous  media,  is  an  important  component  of 
predicting  the  transport  and  fate  of  these  compounds  in  the  subsurface.  In  particular,  flow 
situations  in  which  three  mobile  phases  (such  as  water,  NAPL,  and  air)  exist  in  the  porous 
media  are  inherently  complex.  Unfortunately,  the  relatively  low  source  intensities  and  the 
nontunable  source  energies  make  traditional  dual  gamma  techniques  unsuitable  to  study 
flow  instabilities  which  can  change  within  seconds.  We  present  an  alternate  technique, 
which  uses  synchrotron  X  rays  from  the  Cornell  High  Energy  Synchrotron  Source 
(CHESS)  to  measure  three-phase  fluid  saturations  on  the  time  scale  of  seconds.  Using  the 
harmonic  content  resulting  from  X  ray  diffraction,  we  obtained  a  high-intensity  X  ray 
beam  consisting  of  distinct  tunable  energies.  Three-phase  saturations  were  measured  on 
5-s  timescales  daring  fingering  of  light  NAPL  into  regions  of  dry  and  water  wet  sandy  soil. 
In  the  water  wet  soil  the  oil  finger  was  less  saturated,  slower,  and  wider  than  the  same 
finger  in  the  dry  soil.  The  results  yield  insights  into  the  nature  of  three-phase  preferential 
flow. 


Introduction 

Accurate  measurements  of  multiphase  flow  in  porous  media 
are  needed  to  understand  and  model  subsurface  transport  of 
contaminants  such  as  nonaqueous  phase  liquids  (NAPLs).  In 
particular,  flow  instabilities  and  the  subsequent  preferential 
flow  patterns  can  dramatically  affect  the  fate  of  NAPLs  and  the 
subsequent  attempts  at  remediation.  For  quantitative  studies 
of  preferential  flow  it  is  imperative  to  obtain  fast,  accurate,  and 
nondestructive  measurements  of  various  fluid  saturations 
within  the  porous  media.  Nondestructive  measurements  of 
multiphase  saturations  can  be  achieved  using  dual-energy 
gamma  radiation  systems,  but  these  techniques  have  been  lim¬ 
ited  in  speed  [Oostrom  et  a/.,  1995]. 

Attenuation  of  high-energy  electromagnetic  radiation,  such 
as  X  and  gamma  rays,  is  directly  related  to  the  amount  of  fluid 
or  solid  in  the  path  of  the  radiation,  thus  providing  a  means  to 
obtain  nondestructive  fluid  saturation  measurements.  Radia¬ 
tion  sources  and  techniques  have  progressed  over  the  years, 
providing  quicker  and  more  accurate  saturation  measure¬ 
ments.  Decays  of  radioactive  elements  such  as  cesium  (Cs)  and 
americium  (Am)  produce  monochromatic  gamma  rays,  and 
the  elements  have  been  used  as  portable  radiation  sources  for 
field  and  laboratory  measurements  {HiUeU  1980].  Combining 
and  collimating  Cs  and  Am  sources  allows  them  to  be  used  for 
three-phase  measurements  [Nopiger  and  Swartzendruber,  1974; 
Hopmans  and  Dane,  1986;  F errand  et  al.,  1986;  Lenhard  et  aL, 
1988;  lllangasekare  et  aL,  1995].  Recently,  McBride  and  Miller 
[1994]  have  developed  a  technique  that  obtains  dual  energies 
from  an  X  ray  tube.  Both  sources  can  yield  accurate  three- 
phase  saturations  but  only  over  long  timescales  as  enough 
photon  counts  are  needed  to  overcome  the  random  errors 
associated  with  counting  statistics  [Gardner  et  aL,  1972; 
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Oostrom  et  aL,  1995].  Recently,  synchrotron  radiation  sources 
have  also  been  used  for  nondestructive  measurements.  The 
greater  than  10,000-fold  intensity  increase  over  tube  sources 
and  radioactive  elements  produces  fast  two-phase  saturation 
measurements  and  allows  preferential  flow  to  be  precisely 
monitored  in  real  time  [Liu  et  aL  1993]. 

In  the  first  part  of  the  paper,  we  present  a  general  procedure 
for  obtaining  and  using  dual  energies  from  a  synchrotron  X  ray 
source.  This  procedure  has  been  used  to  simultaneously  mea¬ 
sure  water,  NAPL,  and  air  saturations  in  a  sandy  soil.  For 
saturation  data  taken  in  5-s  time  intervals,  the  measured  vari¬ 
ance  in  the  volumetric  content  is  approximately  0.01  cm^/cm^, 
which  is  comparable  to  the  destructive  oven  drying  technique. 
This  technique  can  be  applied  to  produce  fast,  precise  moni¬ 
toring  of  any  three  phases  in  porous  media.  In  the  second  part 
of  the  paper,  we  present  and  analyze  the  saturations  obtained 
when  a  light  NAPL,  hereafter  referred  to  as  oil,  infiltrates 
through  dry  and  partially  water  wetted  soil. 

Theory 

For  monochromatic  radiation  the  attenuation  of  high- 
energy  radiation  is  linearly  related  to  the  material  within  the 
radiation’s  path  [Cullity,  1956], 

A  =  \a  {T/D  = (1) 

i 

where  A  is  the  measured  attenuation  defined  as  the  natural 
logarithm  of  the  transmitted  intensity  T,  divided  by  the  inci¬ 
dent  intensity  /.  The  sum  is  over  each  material  in  the  radia¬ 
tion’s  path,  where  U  is  the  attenuation  constant  and  x  is  the 
cumulative  thickness  of  each  material.  The  attenuation  con¬ 
stants  depend  on  the  atomic  constituents  of  each  material  and 
the  energy  of  the  radiation.  In  general,  the  attenuation  de¬ 
creases  with  increasing  radiation  energy  and  decreasing  atomic 
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number,  with  the  exception  of  abrupt  attenuation  increases  at 
atomic  K  and  L  electron  energy  levels. 

Substituting  in  the  relevant  materials  for  multiphase  mea¬ 
surements  in  a  porous  media  of  total  thickness  x, 

A  =  -  UoO^x  *•  U,e,x  ~  U,  +  D.  (2) 

Here  x  is  the  thickness  of  the  experimental  chamber,  6  are  the 
volumetric  saturations,  and  the  subscripts  w,  o,  and  s  denote 
the  water,  oil  (or  NAPL).  and  soil  phases,  respectively.  Atten¬ 
uation  by  the  chamber  wails  and  the  air  between  the  detectors 
is  given  by  the  constant  The  term  D  is  a  constant  offset  that 
occurs  if  relative  rather  than  absolute  intensities  are  measured. 
Attenuation  by  the  air  inside  the  chamber  can  be  considered 
negligible.  In  two-phase  flow  situations  where  there  is  only  one 
independent  quantity,  for  example,  water  and  air  in  a  partially 
saturated  soil,  the  equation  can  be  inverted  to  yield  the  fluid 
saturations  in  terms  of  the  measured  attenuation. 


-^0 


(3) 


where  Aois  the  attenuation  measured  for  the  dry  soil.  Notice 
that  only  changes  in  the  attenuation  enter  into  the  fluid  deter¬ 
mination,  and  knowledge  of  the  individual  contributions  of  D 
and  Uf.  to  the  attenuation  is  unnecessary.  In  three-phase  flow 
situations  where  three  saturations  are  variable  (water,  oil  and 
air,  or  water,  soil  and  air),  another  independent  measurement 
is  needed  to  characterize  the  system.  This  can  be  provided  by 
measuring  the  attenuation  of  radiation  of  a  different  energy.  In 
a  sample  in  which  saturations  of  water,  oil,  and  air  are  variable, 
the  fluid  saturations  can  be  determined  from  the  attenuation  of 


^  1  /  UoH  ^Ol\  I  _ _ ^ ^ 

Var  (0H<)  t  \  Ri  Rh /\i^wL^oH  ^  ^wh^ol)x) 

(7) 

,  .  /  1 _ y 

Var  Bo  t  \  Rl  Rh  )\(^yvL^0H  ^  ^wh^ol)x) 


(8) 

Variances  can  be  decreased  by  either  increasing  the  counting 
times,  the  countable  photon  flux,  or  the  denominator  of  the 
third  term.  For  conventional  dual  gamma  sources,  long  count¬ 
ing  times  are  typically  used  for  increasing  precision.  The  high 
flux  and  tunable  energies  available  at  synchrotron  sources  al¬ 
low  the  second  two  factors  to  be  increased  dramatically,  thus 
lowering  the  variance,  and  decreasing  the  necessary  counting 
times. 

For  dual-energy  attenuation  measurements,  the  advantages 
of  large  photon  fluxes  are  more  subtle  than  a  brute  increase  in 
the  counting  rates.  Typically,  large  photon  fluxes  allow  ion 
chambers  to  be  used  as  detectors  for  attenuation  measure¬ 
ments.  Ion  chambers  detect  ionization  currents  produced  when 
X  rays  are  passed  through  an  inert  gas,  and  when  used  with 
synchrotron  source  fluxes  are  free  from  counting  statistics  as¬ 
sociated  with  individual  photon  detectors.  Unfonunately,  the 
lack  of  energy  resolution  precludes  ion  chambers  from  being 
used  for  dual  energy  attenuation  measurements.  Scintillation 
or  solid  state  detectors  are  still  needed,  both  of  which  give 
linear  results  up  to  a  maximum  photon  flux,  but  at  a  synchro¬ 
tron  source  the  detectors  can  be  configured  to  run  at  their 
maximum  rates  for  a  wealth  of  experimental  situations. 


the  high  (A^)  and  low  {Ai^)  energy  radiation, 

OH  ^  ^  WH^Ol}^ 

_  ^ A  U wl(‘^h  ^  wo) 

(UwL^OH~  ^WH^Ol)^ 

where  the  subscripts  H  and  L  denote  the  high-  and  low-energy 
attenuation  constants  andyl/^o  and/l^o  are  the  high  and  low 
attenuations  for  the  dry  sand,  respectively.  These  are  the  sat¬ 
uration  equations  for  any  dual-energy  attenuation  apparatus. 

The  advantages  of  synchrotron  radiation  can  be  demon¬ 
strated  by  studying  the  sources  of  variance  in  the  above  equa¬ 
tions.  The  variance  of  the  fluid  saturations  consists  of  system¬ 
atic  variances  due  to  uncertainties  in  attenuation  constants,  dry 
attenuations,  and  chamber  thickness,  and  temporal  variances 
due  to  the  finite  counting  times  for  the  detectors.  The  system¬ 
atic  variances  can  be  minimized  by  careful  calibration,  and  the 
variance  is  typically  dominated  by  the  temporal  variance 
[Oostrom  et  al,  1995].  Using  Poisson  statistics  for  the  detectors, 
photon  intensities  of  and  Ri^  (photons  per  second),  and  a 
counting  time  of  t  (seconds),  the  variance  of  the  attenuations 
are 

VarM„).i.  («) 


The  advantages  of  tunable  energies  and  attenuations  are 
much  greater  than  the  pure  increase  in  intensity.  The  attenu¬ 
ation  of  each  fluid  can  be  increased  in  a  number  of  ways; 
doping  the  fluids  with  heavy  elements,  increasing  the  path 
length,  or  using  lower  energy  X  rays.  For  conventional  radia¬ 
tion  sources,  fluid  doping  is  often  used,  but  any  method  of 
increasing  the  attenuation  is  limited  as  it  cuts  into  the  detected 
count  rate.  For  synchrotron  sources,  increasing  the  attenuation 
is  less  limited  as  the  large  incident  intensities  assure  plenty  of 
transmitted  photons.  Also,  in  some  situations  the  thickness  of 
the  chamber  and  fluid  contents  are  fixed.  In  this  case,  the 
attenuations  can  be  optimized  by  altering  the  energy  of  the  X 
rays.  This  is  impossible  for  conventional  sources,  but  the  white 
light  emitted  by  synchrotron  sources  allows  a  choice  from  a 
continuum  of  X  ray  energies. 

Most  importantly,  precision  can  be  increased  by  selecting 
the  X  ray  energies  such  that  the  specific  oil  and  water  attenu¬ 
ations  combine  to  yield  a  large  value  for  the  denominator  of 
the  third  term.  This  result  is  achieved  if  one  fluid  preferentially 
attenuates  the  low-energy  X  ray  and  the  other  fluid  preferen¬ 
tially  attenuates  the  high-energy  X  ray  such  that  » 

UwH^oL'  situation  occurs  only  when  an  element  within 
one  of  the  fluids  has  an  absorption  edge  between  the  two 
working  energies.  The  selection  of  energies  above  and  below 
absorption  edges  is  again  trivial  for  synchrotron  sources. 


For  synchrotron  sources  the  variance  is  typically  twice  this 
amount  as  counting  uncertainties  apply  to  both  the  incident 
and  transmitted  intensities.  Using  error  propagation  formulas, 
the  variance  in  the  fluid  saturations  can  be  calculated  [Oostrom 
et  aL  1995], 


Materials  and  Methods 

Figure  1  depicts  a  schematic  diagram  of  the  apparatus  used 
on  the  F-2  beam  line  at  the  Cornell  High  Energy  Synchrotron 
Source  (CHESS).  CHESS  is  a  National  Science  Foundation 
funded  synchrotron  facility  open  to  outside  users  through  peer 
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reviewed  proposals  (CHESS  Facility  Description,  available  at 
http://www.chess.comell.edu,  1995).  The  initial  white  X  ray 
beam  is  reflected  off  a  double-bounce  Si(220)  monochromator, 
producing  a  beam  of  two  distinct  energies.  The  multiple- 
energy  beam  entered  the  experimental  hutch  where  it  passed 
through  a  set  of  tantalum  slits  which  deflned  a  beam  size  of  1 
mm  vertical  X  8  mm  horizontal.  Subsequently,  the  beam 
passed  through  a  nitrogen-filled  ion  chamber,  70  cm  of  open 
air,  the  sample  chamber,  and  another  70  cm  of  open  air  before 
hitting  the  beam  stop  at  the  end  of  the  hutch.  X  ray  intensities 
were  measured  with  scintillation  detectors  pointing  down  at 
the  X  ray  beam.  The  incident  detector  was  placed  33  cm  before 
the  sample  chamber,  and  the  transmitted  deteaor  was  placed 
40  cm  after  the  sample  chamber.  The  individual  components  of 
the  apparatus  are  discussed  in  detail  below. 

Synchrotron  Radiation 

When  relativistic  charged  particles  are  bent  around  a  syn¬ 
chrotron  ring  they  emit  well-collimated  X  rays  of  a  broad 
energy  (white)  spectrum.  This  is  in  contrast  to  radioactive 
decay  sources  which  emit  X  rays  with  a  well-defined  fixed 
energy,  and  X  ray  vacuum  tube  sources  which  emit  both  X  rays 
with  a  well-defined  energy  (corresponding  to  the  characteristic 
emission  lines  of  the  target  metal)  and  a  low-intensity,  broad 
energy,  X  ray  background  (bremsstrahlung).  At  the  F-2  line  at 
CHESS  the  radiation  intensity  is  roughly  constant  from  9  to  40 
keV  and  only  drops  by  a  factor  of  2  at  60  keV  (CHESS  Facility 
Description,  1995).  For  absorption  measurements  the  white 
radiation  can  be  made  monochromatic  by  diffraction  from  a 
pair  of  perfect  silicon  crystals.  The  wavelength  A  of  the  X  rays 
diffracted  by  the  first  crystal  is  dependent  on  the  distance 
between  the  lattice  planes  d  and  the  incident  X  ray  angle  0, 
through  Bragg’s  law  [Cwl/fry,  1956], 

nA  =  Id  sin  B  (9) 

where  n  is  an  integer.  This  wavelength  is  related  to  the  X  ray 
energy, 
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Figure  1.  A  schematic  of  the  experimental  set  up  at  the  F-2 
hutch  of  CHESS.  The  monochromator  transforms  the  white  X 
ray  beam  into  a  multiple-energy  beam.  Scintillation  detectors 
measure  the  intensity  of  both  energy  X  rays  by  detecting  air 
scattering  from  the  main  beam.  The  sample  chamber  is 
mounted  on  a  movable  platform  allowing  measurements  from 
many  different  positions  within  the  chamber. 


robust  energy  setting.  When  tuned  to  20  keV,  the  low-  and 
high-energy  X  rays  were  20  and  40  keV,  respectively.  When 
tuned  to  12.5  keV,  the  fundamental  energy  was  eliminated  by 
placing  a  0.3-cm-thick  aluminum  plate  in  the  beam  path,  and 
the  low-  and  high-energy  X  rays  were  25  and  37.5  keV,  respec¬ 
tively.  In  either  case,  the  energies  were  chosen  to  be  above  and 
below  Iodine’s  absorption  edge  at  33.2  keV.  As  measured  by 
the  ion  chambers,  typical  initial  intensities  were  10' '  photons/s 
for  the  low-energy  beam,  and  10'°  photons/s  for  the  high- 
energy  beam. 

X  ray  Detection 

As  the  X  ray  intensity  is  not  constant  for  synchrotron 
sources,  the  incident  and  transmitted  intensities  of  both  ener¬ 
gies  must  be  accurately  measured.  Detection  methods  must 
have  adequate  energy  resolution,  enough  counts  for  good  sta¬ 
tistics,  and  not  destroy  the  incident  beam.  Also  to  reduce 
systematic  errors,  it  is  useful  if  the  incident  and  transmitted 


E  =  hcf\  (10) 

where  h  is  Planck’s  constant  and  c  is  the  speed  of  light.  The 
energies  of  the  X  ray  beam  can  be  varied  by  varying  the  angle 
of  the  monochromator.  The  second  silicon  crystal  is  placed 
below  and  downstream  of  the  first  crystal  and  set  to  an  iden¬ 
tical  angle.  This  standard  dual  monochromator  setup  allows 
the  energies  to  be  varied  without  varying  the  spatial  position  of 
the  beam. 

The  fundamental  energy  corresponding  to  n  =  1  is  the 
strongest  reflection  and  typically  the  only  one  used.  The  higher 
harmonic  energies,  corresponding  to  n  =  2,  3,  4,  •  •  • ,  are  of 
lower  intensity  due  to  the  smaller  band  pass  of  the  monochro¬ 
mator  crystal  at  higher-order  crystal  reflections  (the  band  pass 
scales  as  approximately  w”^),  and  from  the  lower  flux  at 
higher  energies  in  the  incident  beam.  Harmonic  energies  are 
usually  only  observed  at  synchrotron  facilities  due  to  the  broad 
energy  spectra  of  the  incident  beam.  In  fact,  harmonic  energies 
are  typically  a  nuisance  as  they  destroy  the  monochromaticity 
of  the  beam,  and  various  techniques  can  be  used  to  suppress 
them.  For  dual-energy  attenuation  measurements  the  harmon¬ 
ics  are  extremely  useful,  as  they  can  provide  a  simultaneous 
second  (and  possibly  third!)  attenuation  measurement. 

In  our  experiment  we  used  a  fundamental  energy  of  20  keV, 
but  12.5  keV  was  also  tested  in  order  to  determine  the  most 


detectors  are  of  the  same  type. 

We  have  used  Na  I(Tl)  scintillation  detectors,  which  count 
individual  X  ray  photons  at  rates  up  to  10^  counts/s  with  an 
energy  resolution  of  40%  \Tlwmpson^  1986].  Since  the  intensity 
of  the  beam  was  so  great,  the  detectors  were  placed  above  the 
beam  and  pointed  down  at  the  passing  beam.  In  this  fashion, 
the  detectors  counted  the  X  ray  photons  which  were  vertically 
scattered  from  the  beams  by  the  air  in  the  hutch.  Scattering  by 
the  air  consists  of  elastic  (Rayleigh)  and  inelastic  (Compton) 
processes  [Cullity,  1956].  At  the  working  photon  energies  the 
energy  loss  due  to  90°  Compton  scattering  is  less  than  1  keV, 
much  less  than  the  difference  in  photon  energies  and  the  en¬ 
ergy  resolution  of  the  detectors.  More  importantly,  both  pro¬ 
cesses  are  linear  with  the  X  ray  intensity.  To  control  the  mag¬ 
nitude  of  the  scattered  flux  into  the  detector  a  lead  aperture 
and  an  aluminum  attenuator  were  placed  in  front  of  each 
detector.  The  low-  and  high-energy  photon  flux  scattered  into 
the  detectors  are  given  by 

Idl  =  {o’(e.  El)£1  exp  (11) 

Io„  =  {o-(e,  £„)ft  exp  (12) 

where  and  are  the  low  and  high  X  ray  intensities  into 
the  detector,  7^^  and  are  the  low  and  high  intensities  of  the 
beam  below  the  detector,  (t  is  the  air  scattering  cross  section 
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and  is  a  function  of  scattering  angle  0  (90°  in  this  case)  and  X 
ray  energy  £,  ft  is  the  solid  angle  of  ^e  aperture,  is  the 
attenuation  constant  (001“^)  of  the  attenuator  in  front  of  the 
deteaor  and  y  is  the  thickness  (centimeters)  of  the  attenuator. 
Since  the  aluminum  preferentially  attenuates  the  low-energy  X 
ray,  by  adjusting  the  attenuator  thickness  y  and  the  aperture 
size  ft,  we  were  able  to  obtain  roughly  15,000  counts/s  of  high- 
and  low-energy  photons  scattered  into  each  detector.  To  elim¬ 
inate  stray  radiation  from  entering  the  detector  a  circular  lead 
flight  tube  was  attached  to  the  front  of  each  detector.  The 
incident  detector  was  placed  10  cm  above  the  beam,  had  a  3-cm 
flight  tube,  a  1.2-cm  circular  aperture,  and  a  0.6-cm-thick  alu¬ 
minum  attenuator.  The  transmined  detector  was  placed  8.5  cm 
above  the  beam,  had  a  7-cm  flight  tube,  a  2.2-cm  circular 
aperture,  and  no  attenuator.  Throughout  each  experiment  the 
configuration  of  the  detectors  was  not  altered  as  this  would 
change  the  proponionality  constants  and  destroy  any  previous 
calibration. 

The  output  of  each  scintillation  detector  was  fed  through  an 
attached  preamplifier  followed  by  a  variable  gain  spectroscopy 
amplifier.  The  gains  of  the  spectroscopy  amplifiers  were  set 
such  that  the  electronic  pulses  of  average  heights  3.4  and  6.8  V 
were  produced  for  the  20  and  40  keV  photons,  respectively. 
The  output  from  each  amplifier  was  split  and  fed  into  two 
separate  single  channel  analyzers,  one  analyzer  recorded 
pulses  between  2.5  and  4.5  V  as  20  keV  photons,  the  other 
recorded  pulses  between  5.2  and  7.7  V  as  40  keV  photons. 

For  any  dual  gamma  setup  there  is  a  finite  chance  that  a 
high-energy  photon  may  elicit  a  detector  response  which  will 
be  counted  in  the  low-energy  photon  bin,  and  vice  versa.  The 
sources  of  this  “cross-binning”  are  different  for  synchrotron  X 
rays  than  for  traditional  dual  gamma  apparatus  and  are  worth 
discussing.  In  a  traditional  dual  gamma  apparatus  using  Am 
and  Cs  sources,  the  “cross-binning”  is  due  to  Compton  scat¬ 
tering.  When  a  660  keV  enters  the  scintillation  detector,  it  may 
be  Compton  scattered  out  of  the  detector,  leaving  behind  a 
Compton  electron  with  energy  ranging  from  0  to  477  keV.  The 
Compton  electron  subsequently  produces  scintillation  light, 
resulting  in  a  broad  low-energy  tail  in  the  detected  energy 
spectra.  If  the  Compton  electron  happens  to  acquire  60  keV  of 
energy,  it  will  be  detected  as  a  60  keV  photon  and  will  be 
mistaken  for  a  Cs  photon.  For  the  lower  photon  energies  used 
here,  Compton  scattering  inside  the  detector  does  not  cause 
errant  readings.  The  maximum  Compton  electron  energy  cre¬ 
ated  by  a  high-energy  40  keV  photon  is  5  keV,  which  will  not 
be  counted  as  a  low-energy  20  keV  photon.  Instead  the  more 
likely  source  of  cross-binning  at  these  energies  is  due  to  the 
statistical  spread  in  the  small  amount  of  primary  electrons 
produced  for  each  X  ray,  resulting  in  an  energy  full  width  at 
half  maximum  (FWHM)  of  40%.  Thus  there  still  exists  a  finite 
chance  that  a  40  keV  photon  will  produce  a  response  in  the 
detector  that  would  be  fit  into  the  bin  reserved  for  counting  20 
keV  photons.  For  our  experimental  parameters  the  expected 
chance  of  a  high-energy  photon  being  counted  as  a  low-energy 
photon,  and  vice  versa  is  approximately  1%.  As  shown  in  the 
appendix,  the  cross-binning  creates  an  additional  variance  of 
1%  in  the  fluid  concentrations. 

Sample  Chamber  and  Fluids 

The  sample  chamber  had  0.94-cm-thick  polycarbonate  walls 
and  interior  dimensions  of  30  cm  wide  {x  direction),  0.94  cm 
thick  (y  direction),  and  55  cm  tall  (2  direction).  Attached  to 
the  bottom  end  was  a  manifold  of  five  fluid  ports,  and  the  top 


Table  1.  Measured  Attenuation  Constants  of  the  Water 
and  Oil  Phases 


Attenuation,  cm"^ 

Phase 

20  keV 

40keV 

Water 

2.56 

0.512 

Oil 

1.339 

1.102 

end  was  left  open  to  the  atmosphere.  The  chamber  was  filled 
with  clean,  dry,  sieved,  quartz  sand  (Unimin  Corporation)  us¬ 
ing  a  dual  randomizing  funnel.  Grain  sizes  used  were  12/20, 
20/30,  and  30/40.  where  the  numbers  correspond  to  the  sieve 
sizes  the  sand  passes  through  and  does  not  pass  through,  re¬ 
spectively.  To  measure  the  fluid  saturations  at  many  points,  the 
chamber  was  mounted  on  a  movable  platform  which  had  a 
horizontal  and  vertical  range  of  25  and  50  cm,  respectively,  and 
a  position  repeatability  of  0.01  cm. 

Distilled  water  doped  with  50  g/L  NaBr  and  FD&C  blue  dye 
was  used  for  the  aqueous  phase,  and  Soltrol  220  (Phillips 
Company)  doped  with  77  g/L  iodoheptane  (Aldrich  Chemical 
Company)  and  a  red  Sudan  IV  dye  was  used  for  the  oil  phase. 
Soltrol  220  has  been  widely  used  to  study  NAPL  behavior 
[Schroth  et  aL,  1995].  The  addition  of  the  dopants  raised  the 
specific  gravities  of  the  fluids  slightly  to  0.825  and  1.040  g/cm^, 
for  the  oil  and  water  phases,  respectively.  Measured  changes  in 
the  surface  tensions  were  less  than  1  dyne/cm  (du  Nouy  ring 
method),  with  surface  tensions  of  28  and  72  dynes/cm  for  the 
oil  and  water  phases,  respectively,  and  an  interfacial  tension  of 
32  dynes/cm.  Thus  the  oil  phase  was  a  spreading  NAPL.  Huid 
attenuation  constants  were  obtained  using  a  polycarbonate 
chamber  with  five  sequential  cells  of  thickness  0.56  cm  [Lcn- 
hard  et  aL,  1988].  The  attenuation  was  measured  for  the  empty 
chamber  and  after  each  cell  was  filled  with  either  fluid.  The 
measured  attenuations  were  linear  with  fluid  thickness,  except 
at  the  highest  thickness  (corresponding  to  2.8  cm  of  fluid,  or  a 
volumetric  saturation  of  3.00  cm^/cm^)  where  the  transmitted 
count  rates  were  smallest.  Attenuation  constants  were  calcu¬ 
lated  from  the  slope  of  the  attenuation  versus  thickness  curve. 
Table  1  lists  the  measured  attenuation  constants  and  their 
precision. 

To  test  the  precision  and  speed  of  the  technique  and  to 
observe  unstable  three-phase  flow,  the  saturations  within  an  oil 
finger  penetrating  partially  water  saturated  soil  were  mea¬ 
sured.  Initial  attenuations  of  the  dry  sand-packed  chamber 
were  measured  at  various  positions  throughout  the  chamber. 
The  chamber  was  then  filled  from  below  with  water  at  a  rate  of 
30  mL/min  to  a  height  of  2  =  40  cm,  and  subsequently  drained 
at  a  rate  of  20  mL/min  until  the  water  fringe  was  at  a  height  of 
2  =  12  cm.  The  chamber  had  three  distinct  regions,  below  2  = 
12  cm  the  sand  pack  was  water  saturated  (6^  =  0.37  cm^/ 
cm^),  between  2  =  12  cm  and  2  =  40  cm  the  sand  was  drained 
to  the  residual  water  saturation  (0/?  =  0.05  cm^/cm^),  and 
above  2  =  40  cm  the  sand  was  dry.  Oil  was  applied  through  a 
point  source  at  the  top  of  the  chamber  at  a  rate  of  4  mL/min. 
A  finger  formed  which  moved  quickly  through  the  dry  sand, 
and  slowly  through  the  partially  wet  sand  before  finally  hitting 
the  water  table  and  spreading  laterally.  Figure  2  presents  a 
sketch  of  the  chamber  during  the  finger  infiltration.  Temporal 
fluid  saturations  were  measured  in  the  water  unsaturated  and 
partially  water  saturated  regions  at  positions  2  =  45  cm  and 
z  =  30  cm,  respectively,  as  the  finger  traversed  these  regions. 
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Figure  2.  Front  view  of  the  experimental  chamber  during  an 
oil  infiltration.  The  oil  was  applied  from  a  point  source  at  the 
top  and  the  water  from  the  manifold  of  ports  at  the  bottom. 
Initially,  below  z  =  12  cm  the  sand  pack  was  water  saturated 
(0^  =  0.37  cm^/cm^),  between z  =  12  cm  andz  =  40  cm  the 
sand  was  drained  to  the  residual  water  saturation  (0/?  =  0.05 
cm^/cm-^).  and  above  z  =  40  cm  the  sand  was  dry.  The  initial 
approximate  water  content  as  a  function  of  depth  is  sketched 
on  the  right.  The  finger  became  wider  and  slower  when  it 
entered  the  partially  water  saturated  sand.  Oil  and  water  sat¬ 
urations  were  observed  as  a  function  of  time  at  z  =  45  cm  and 
z  =  30  cm. 


Horizontal  profiles  of  the  finger  were  taken  at  z  ==  30  cm  when 
the  finger  saturations  had  reached  steady  state. 

Results  and  Discussion 

Time  Resolution 

Figures  3  and  4  show  the  measured  fluid  saturations  as  the 
oil  finger  penetrated  the  initially  dry  and  partially  water  satu¬ 
rated  20/30  sand,  respectively.  To  resolve  the  tip  of  the  finger 
passing  through  the  observed  region,  attenuation  data  were 
accumulated  over  5-s  intervals.  The  predicted  variances  in  the 
fluid  contents  due  to  the  finite  counting  times  can  be  calculated 
from  (8)  and  (9).  Even  for  the  5-s  counting  times  the  variances 
are  only  0.01  cm^/cm^  and  are  depicted  by  the  error  bars  in 
Figures  3  and  4.  The  scatter  in  the  data  roughly  matches  the 
predicted  variance,  implying  that  the  expected  precision  has 
been  achieved,  and  that  the  counting  statistics  are  the  precision 
limiting  factor. 

Figures  5  and  6  show  the  measured  fluid  saturations  as  the 
oil  finger  penetrated  the  initially  dry  and  partially  water  satu¬ 
rated  30/40  sand,  respectively.  Data  were  accumulated  over  3-s 
intervals,  yielding  slightly  higher  variances,  but  greater  time 
resolution  than  for  the  results  in  Figures  3  and  4.  The  scatter 
in  the  data  is  slightly  larger  than  expected  from  counting  sta¬ 
tistics,  and  was  most  likely  caused  by  temporary  cooling  prob¬ 
lems  with  the  monochromator  crystal. 

It  is  instructive  to  compare  the  precision  and  speed  of  the 
above  synchrotron  X  ray  technique  and  other  dual  energy 
techniques.  Using  Am  and  Cs  sources,  Nofidger  and  Swanzen- 
druder  [1974]  obtained  a  water  precision  of  0.01  cm^/cm^  when 
counting  for  5  minutes,  and  0.04  cm^/cm^  when  counting  for 
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Figure  3.  Water  and  oil  contents  versus  time  as  the  oil  finger 
passes  through  initially  dry  20/30  sand.  Error  bars  depict  the 
expected  variance  due  to  counting  statistics  as  calculated  from 
(8)  and  (9),  and  they  roughly  match  the  scatter  in  the  data.  The 
rapid  rise  of  the  oil  content  can  be  easily  discriminated,  and  the 
oil  level  approaches  saturation. 


5  s.  Oostrom  ei  aL  [1995]  obtained  water  and  oil  precisions  of 
better  than  0.01  cm^/cm^  when  counting  for  60  s,  with  slightly 
worse  precisions  for  a  counting  time  of  30  s.  Studies  [Okuda  et 
aL,  1996;  Imhoff  et  aL,  1996]  using  the  X  ray  vacuum  tube 
apparatus  described  by  McBride  and  Miller  [1994]  have  re- 
poned  precisions  of  0.0024  cm^/cm^  for  porosity  and  0.005 
cm^/cm^  for  PCE  saturation  when  counted  for  6  s  three  sepa¬ 
rate  tinies.  It  is  unclear  if  this  reported  accuracy  is  from  a 
dual-energy  measurement  or  a  single  energy  measurement,  the 
latter  of  which  is  inherently  more  accurate,  as  the  attenuations 
do  not  have  to  be  decoupled.  These  comparisons  are  only 
rough  as  all  the  studies  have  used  different  fluids,  dopants,  and 
chamber  thickness. 

One  of  the  strengths  of  the  synchrotron  technique  can  be 
seen  from  a  comparison  of  Figures  4  and  6.  In  three-phase 
studies  it  can  be  important  to  know  which  phase  is  forced  out 
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Figure  4.  Water  and  oil  contents  versus  time  as  the  oil  finger 
passes  through  partially  water  saturated  20/30  sand.  The  oil 
content  rises  less  rapidly  and  to  a  lower  saturation  than  in  the 
dry  sand.  The  water  content  is  unaffected  by  the  oil  finger. 
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Time  (sec) 

Figure  5.  Water  and  oil  contents  versus  time  as  the  oil  finger 
passes  through  initially  dry  30/40  sand.  Again,  the  oil  rises 
rapidly  to  a  value  near  saturation. 


Figure  7.  A  horizontal  scan  through  the  oil  finger  in  the 
partially  water  saturated  30/40  sand.  The  lines  are  moving 
averages  of  the  saturations  and  are  a  guide  for  the  eye.  The 
extra  scatter  in  the  data  is  most  likely  due  to  local  variations  in 
the  sand  density  due  to  the  coarseness  of  the  sand. 


of  a  region  when  a  third  phase  enters  the  region  {Kao  and 
Hunt.  1996].  For  the  20/30  sand  the  water  content  remains 
constant  as  the  oil  finger  passes  implying  that  only  air  leaves 
the  volume,  but  in  30/40  sand  the  water  content  is  driven  down 
slightly  by  the  arrival  of  the  oil  finger.  These  small  and  fast 
saturation  changes  are  easily  observed. 

Figure  7  shows  the  measured  horizontal  profile  of  the  oil 
finger  in  30/40  sand.  In  this  case,  the  measured  variance  was 
much  greater  than  expected  for  counting  statistics.  This  vari¬ 
ance  was  even  greater  in  the  coarser  12/20  and  20/30  sands. 
The  extra  scatter  was  most  likely  due  to  the  variations  in  the 
soil  attenuation  y4  Q.  The  soil  volume  sampled  by  the  X  rays  was 
only  0.08  cm*^,  which  was  not  much  greater  than  the  average 
grain  size  volumes  for  the  12/20  and  20/30  sand.  Thus  the  sand 
density  and  porosity  fluctuates  on  this  length  scale,  causing 
fluctuations  in  the  dry  soil  attenuation.  If  the  sand  settles  and 


700  800  900  1000  1100  1200  1300 

Time  (sec) 


shifts  slightly  when  wetted,  the  sand  attenuation  at  each  posi¬ 
tion  can  change,  thereby  creating  more  scatter  in  the  data. 

Three-Phase  Unstable  Flow 

Quantitative  determinations  of  oil  finger  properties  can  be 
made  using  this  technique.  Finger  widths  were  measured  by 
taking  the  full  width  at  half  maximum  (FWHM)  in  the  hori¬ 
zontal  profiles.  Finger  saturations  were  measured  as  the  aver¬ 
age  near  the  center  of  the  profile.  These  results  are  summa¬ 
rized  in  Table  2.  In  the  initially  dry  region  the  oil  saturation 
was  close  to  full  saturation  (0.30  cm^/cm^)  and  the  finger  was 
2.5-3.0  cm  wide  (FWHM).  In  the  partially  water  saturated 
region  the  oil  saturation  was  0.15  cm^/cm^.  roughly  one  half  of 
the  completely  saturated  value,  the  finger  was  9.0-10.0  cm 
wide  (FWHM).  and  the  water  saturation  has  a  value  of  0.05 
cm^/cm^.  So  not  only  does  the  water  remain  when  the  oil  finger 
passes,  the  initial  water  content  creates  a  much  lower  oil  sat¬ 
uration  (and  thus  higher  air  situation)  within  the  oil  finger. 

Comparisons  can  be  made  between  this  three-phase  exper¬ 
iment  and  a  similar  two-phase  experiment  performed  by  Liu  et 
aL  [1991].  In  Liu’s  experiment,  finger  properties  were  mea¬ 
sured  as  a  water  finger  passed  from  initially  dry  sand  into 
partially  water  wetted  sand.  Similar  qualitative  effects  were 
observed:  in  the  initially  dry  sand  the  finger  was  2.3  cm  wide 
(FWHM)  with  a  saturation  of  0.35  cm^/cm^,  while  in  the  par¬ 
tially  wetted  sand  the  finger  was  5.0  cm  wide  (FWHM)  with  a 
saturation  of  0.24  cm^/cm^. 

For  either  invading  fluid  the  increase  in  the  finger  widths  can 
be  understood  through  the  following  argument.  In  the  partially 


Table  2.  Summary  of  Oil  Finger  Properties  in  the  20/30 
and  30/40  Sands  for  Both  Regions 

Finger  Width,  Oil  Saturation,  Water  Saturation, 
Sand  Region  cm  cm^/cm^  cm*Vcm’^ 


Figure  6.  Water  and  oil  contents  versus  time  as  the  oil  finger 
passes  through  partially  water  saturated  30/40  sand.  The  oil 
content  rises  to  a  lower  saturation  than  in  the  dry  sand,  similar 
to  that  seen  in  20/30  sand.  As  opposed  to  20/30  sand,  in  this 
case  the  water  content  is  driven  lower  by  the  oil  finger. 


Dry  20/30 

2.5 

0.30 

0.00 

Dry  30/40 

3.0 

0.30 

0.00 

Water  wetted  20/30 

9.0 

0.15 
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Water  welted  30/40 

10.0 

0.15 
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wetted  sand  the  water  forms  menisci  which  in  turn  creates  three-phase  flow,  where  the  real  time  saturation  measurements 


isolated  islands  of  the  water  or  air  phases.  The  invading  fluid 
has  to  either  bypass  these  islands  or  force  the  water  out  of  the 
pores,  resulting  in  a  decrease  of  the  fluid  conductivity.  The 
initial  water  film  covering  the  solid  surfaces  also  enhances  the 
lateral  spreading  of  water  and/or  oil,  the  latter  of  which  is  a 
spreading  fluid  when  in  contact  with  water.  The  simultaneous 
decrease  in  conductivity  (/C)  and  increase  in  sorptivity  (5) 
combine  to  increase  the  finger  width.  Experimentally,  it  is 
observed  that  the  ratio  of  the  finger  widths  between  the  dry 
and  prewetted  regions  is  greater  for  the  oil  finger  (a  factor  of 
3),  than  it  is  for  the  water  finger  (a  factor  of  2).  This  difference 
in  ratios  may  be  attributed  to  the  fact  that  the  water  finger  is 
miscible  with  the  initial  water  in  the  soil,  while  the  oil  is  im¬ 
miscible.  Thus  the  water  flow  is  blocked  or  slowed  by  pores 
with  entrapped  air,  while  the  oil  flow  is  blocked  by  pores  with 
entrapped  air  or  entrapped  water. 

Although  the  air  and  water  phases  partially  block  the  pores, 
the  oil  must  replace  one  or  both  of  these  phases  as  it  invades 
a  region.  For  this  experiment  it  appears  that  generally  the  oil 
only  replaces  the  air  as  it  invades  the  soil  leaving  the  water  in 
place.  The  one  exception  is  in  the  30/40  sand  where  the  oil 
finger  causes  the  water  saturation  to  decrease  from  0.08  to  0.05 
cm^/cm\  Here  the  final  water  saturation  value  is  identical  to 
the  final  saturation  value  for  12/20  and  20/30  sands.  This  sat¬ 
uration  is  also  the  value  of  the  so-called  irreducible  water 
saturation  for  all  the  sands,  roughly  where  the  water  conduc¬ 
tivity  becomes  negligible. 

Summary 

Synchrotron  radiation  opens  many  possibilities  to  increase 
the  speed  of  nondestructive  measurements  in  porous  media. 
For  this  particular  experimental  realization,  better  detector 
mounts  and  standard  attenuators  will  make  the  detectors  eas¬ 
ier  to  configure.  Faster  electronics  can  be  installed  to  increase 
the  detectors’  count  rates  and  the  precision  further.  Large 
improvements  are  possible  by  using  custom  monochromator 
crystals.  Miscut  monochromator  crystals  will  inflate  the  beam 


of  simple  oil  infiltrations  illustrate  how  the  initial  water  content 
makes  the  oil  fingers  less  saturated,  wider,  and  slower  moving 
than  fingers  in  dry  sand. 

Appendix 

If  these  cross-binning  probabilities  are  known,  estimates  can 
be  made  for  the  error  that  cross-binning  introduces  into  the 
measured  attenuations  and  attenuation  constants.  It  is  most 
instructive  to  observe  how  cross-binning  changes  the  expected 
attenuations  of  known  quantities  of  the  pure  fluids.  The  mea¬ 
sured  incident  intensity  1^  of  each  X  ray  can  be  related  to  the 
true  incident  intensity  /  of  each  X  ray  on  the  detertor, 

~  ■*"  Ch^l  (A1) 

(A2) 

where  is  the  probability  of  a  j  energy  photon  being  counted 
as  an  i  energy  photon.  Assuming  c^/c^;,  c^lc^  1,  and 
similar  constants  for  the  transmitted  X  rays,  the  difference 
between  the  measured  attenuation  and  true  attenuation  can  be 
expressed  as 


The  calculation  is  identical  for  the  low-energy  attenuations. 
The  true  attenuation  is  linear  with  fluid  concentration,  while 
the  cross-binning  adds  some  nonlinearity  to  the  measured  at¬ 
tenuation.  This  nonlinearity  results  in  an  effective  attenuation 
constant  that  varies  with  fluid  concentration.  To  determine  the 
magnitude  of  this  nonlinearity,  we  differentiate  the  above  ex¬ 
pression  with  respect  to  the  fluid  concentration.  Substituting  in 
the  dependence  of  the  transmitted  intensities  on  the  water 
concentration  yields  the  effective  attenuation  constant 

TdO) 


in  the  vertical  direction,  thus  diminishing  the  variances  caused 
by  coarse  soils.  One  intriguing  possibility  is  the  use  of  a  hori¬ 
zontally  bent  monochromator.  This  could  produce  focused 
dual  energy  X  ray  beams  which  will  converge  at  the  sample 
position  but  be  separated  in  space  elsewhere  so  that  ion  cham¬ 
bers  can  be  used  for  detection.  This  would  eliminate  counting 
statistics  and  make  attenuation  measurements  arbitrarily  fast. 

In  summary,  the  power  of  synchrotron  radiation  has  been 
used  to  increase  the  speed  of  dual  radiation  techniques  in 
order  to  observe  unstable  three-phase  phenomena.  The  dual 
energies  were  provided  by  the  natural  harmonic  content  of  the 
monochromated  beam,  which  can  be  tuned  to  the  most  appro¬ 
priate  energies.  The  advantages  of  the  synchrotron  method 
over  dual  gamma  are  the  tunability  of  the  energies  and  ex¬ 
traordinarily  higher  intensities.  Tunability  makes  the  method 
easily  adaptable  to  maximize  the  accuracy  for  many  different 
experimental  requirements;  for  example,  the  energies  can  be 
chosen  to  be  above  and  below  any  absorption  edge.  The  high 
fluxes  allow  scintillation  detectors  to  operate  at  their  maximum 
efficiency.  Also  different  X  ray  optic  configurations  may  pro¬ 
duce  even  faster  nondestructive  measurements.  A  limitation  of 
synchrotron  X  rays  is  the  variability  of  the  source,  and  thus  the 
incident  radiation  must  be  measured,  slightly  increasing  the 
experimental  error.  The  method  is  most  applicable  to  unstable 


•  [1  +  -  Uy^H)e.x  +  Oie..x)-l  (A4) 

where  r(0)  are  the  transmitted  intensities  with  no  fluid  in  the 
chamber.  The  first  linear  term  gives  the  effective  attenuation 
constant,  and  the  second  nonlinear  term  determines  how  the 
effective  attenuation  constant  changes  with  fluid  concentra¬ 
tion.  The  percentage  change  in  the  measured  attenuation  con¬ 
stant  between  an  empty  and  water-saturated  chamber  is 

^UwH  ”  ^^y^wX  /  M.e\ 

TTM  ~  77?  fjM  •  (A5) 

Using  the  observed  attenuation  parameters  and  cross  bin¬ 
ning  ratios  for  20  and  40  keV  radiation,  the  maximum  devia¬ 
tion  of  the  attenuation  constants  was  1%.  For  25  and  37.5  keV 
radiation  the  cross  binning  ratios  were  greater,  resulting  in  a 
maximum  deviation  of  the  attenuation  constants  of  2%. 
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Abstract.  The  understanding  of  simple  laminar  flow  in  tubes  has  often  been  used  to  interpret  the 
more  complicated  flow  in  porous  media.  A  study  of  the  motion  of  two  immiscible  liquids  in  closed 
tubes  with  relatively  large  diameter  (>  0.3  cm  i.d),  was  conducted  in  order  to  examine  the  influence  of 
wetting  and  nonwetting  liquids  on  the  flow  behavior.  The  results  indicate  that  the  wetting  properties 
of  the  fluids  with  regard  to  the  tube  wall  have  a  major  efffect  on  the  formation  and  motion  of  long 
bubbles.  A  physically  based  model  was  used  to  predict  the  velocity  and  the  conditions  for  no  motion 
of  bubbles  and  drops  in  tubes.  These  results  were  used  to  interpret  the  nature  of  oil  and  water 
flow  in  porous  media.  Experiments  in  which  the  wetting  liquid  was  displaced  by  the  nonwetting, 
or  vice  versa,  were  conducted  by  injecting  the  displacing  liquid  at  a  constant  flux  at  the  center  of  a 
two-dimensional  chamber  saturated  with  the  displaced  liquid.  The  influence  of  wetting-nonwetting 
characteristics  on  the  quantity  of  liquid  displaced,  the  shape  of  the  interface  between  the  two  liquids, 
and  the  interpretation  of  the  no  motion  radius  in  a  closed  tube  to  the  case  of  a  porous  medium  are 
discussed.  It  would  appear  that  the  no  motion  radius  gives  a  good  indication  of  the  minimum  width 
of  a  nonwetting  penetrating  finger  and  the  maximum  width  of  nonwetting  ganglia  left  by  drainage. 

Key  words:  oil  and  water,  NAPL,  bubbles,  slugs,  entrapped  fluid,  immiscible  displacement,  ganglia, 
vadose  zone,  soltrol. 


1.  Introduction 

Contamination  of  soils,  sediments,  and  groundwater  with  non-aqueous  phase  liq¬ 
uids  (NAPLs)  is  one  of  the  most  widespread  environmental  problems.  For  example, 
the  volume  of  contaminants  of  130  Department  of  Energy  facilities  in  the  United 
States  exceeds  2,500  billion  liters  with  a  projected  life  cycle  costs  greated  than  300 
billion  dollars  (DOE,  1995;  Wilkins  etal,  1995). 

In  situ  remediation  is  one  of  the  few  viable  treatment  technologies  that  offer 
the  potential  for  cost-effective  clean  up.  Multiphase  models  (Bear  and  Bachmat, 
1990)  have  been  used  with  limited  success  to  optimize  in  situ  remediation.  One  of 
the  problems  is  that  the  complex  interactions  of  the  system  parameters  is  not  well 
understood  and,  therefore,  cannot  be  modelled  accurately.  This  paper  addresses  the 
interactions  of  wetting  properties  of  the  pore  wall  and  the  liquid.  Two  investigations 
are  carried  out.  First,  the  motion  of  two  immiscible  liquids  in  a  cylindrical  tube 
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with  various  diameters  is  investigated.  Then,  we  apply  these  results  to  explain 
observation  of  fluid  distribution  in  porous  media. 


2.  Flow  in  a  Single  TYibe 

THEORY:  Past  research  for  movement  of  two  fluids  in  closed  tubes  has  been 
mainly  described  by  the  movement  of  a  light  fluid  (gas)  bubble  rising  through  a 
more  dense  fluid  (liquid)  (see  review  of  Fabre  and  Lin6,  1992).  In  cases  where 
the  motion  of  two  immiscible  liquids  is  considered,  the  wetting  properties  of  both 
liquids  and  tube  wall  greatly  influence  the  formation  and  shape  of  the  bubble  or 
drop.  This  influence  has  not  been  studied  systematically  in  the  past.  Zukoski  (1966) 
described  gas  bubbles  movement  in  closed  tubes  in  detail,  and  made  a  few  casual 
remarks  on  the  present  problem  of  two  immiscible  liquids,  mostly  on  the  velocity 
of  propagation  and  very  little  on  wetting  properties. 

Following  early  work  on  gas  bubbles  (Dumitrescu,  1943;  Davies  and  Taylor, 
1950;  Tung  and  Parlange,  1976;  Bendiksen,  1984),  the  stream  function  ip,  of  the 
flow  outside  the  axis-symmetric  bubble  in  a  closed  tube  is  expressed  in  terms  of 
Bessel  functions  of  the  first  kind  and  first  order  Ji  (A:„ra"^ ) 

71=  1 

where  kn  is  the  nth  root,  U  is  the  velocity  of  the  bubble  or  drop,  a  is  the  tube 
radius  and  r  is  the  distance  from  the  axis.  In  these  early  theories,  different  number 
of  coefficients.  An,  and  physical  approaches  were  used  to  predict  both  the  bubble 
shape  {ip  =  0),  and  the  bubble  velocity.  Davies  and  Taylor  (1950)  obtained  a 
solution  by  retaining  only  the  first  term  of  the  series  expansion,  which  would 
result  in  a  velocity  increase  when  surface  tension  is  taken  into  accoxmt,  contrary  to 
experimental  observations.  According  to  Fabre  and  Line  ( 1 992),  Dumitrescu  ( 1 943) 
was  the  first  to  give  the  solution  by  retaining  three  terms  in  the  series  expansion 
of  the  potential  flow  around  a  prescribed  spherical  front  ignoring  surface  tension. 
Bendiksen  (1984)  followed  the  method  of  Dumitrescu  taking  into  account  surface 
tension. 

Tung  and  Parlange  (1976)  calculated  the  coefficients  A\  =  0.2516,^2  = 
-0.0040,  and  A3  =  0.0064  in  the  series  (Equation  (1)),  using  an  expansion  of 
Ip  for  X  and  r  small  near  the  tip  of  the  bubble  and  by  satisfying  the  boundary  con¬ 
dition  of  a  constant  pressure  along  the  surface  of  the  gas  bubble,  and  the  maximum 
velocity  principle  of  Garabedian  (1957)  (Figure  2).  The  predictions  of  Tung  and 
Parlange  (1976)  and  Bendiksen  (1984)  of  the  velocity  of  an  air  bubble  in  a  liq¬ 
uid  both  describe  very  well  the  experimental  results  of  Zukoski  (1966)  and  others 
(Fabre  and  Line,  1992),  but  the  first  model  is  simpler,  and  slightly  better  as  it 
satisfies  Garabedian’s  conditions. 
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Figure  1.  Contour  plot  of  the  stream  function  t/?  (Equation  (2))  as  a  function  of  r,  the  distance 
from  the  center  of  the  tube,  and  x,  the  vertical  coordinate  measured  from  the  tip  of  the  bubble, 
using  the  coefficients  A\  =  0.2516,  Ai  =  —0.0040,  and  A^  =  0.(X)64.  -0  =  0  (thick  line) 
indicates  the  calculated  surface  of  the  bubble,  while  the  measured  and  scaled  bubble  (or  drop) 
shapes  fall  into  the  shaded  area. 


The  velocity  of  the  nonwetting  liquid  bubble  rising  in  a  heavier  wetting  medium, 
or  a  liquid  drop  falling  through  a  lighter  medium  is  calculated  using  Bernoulli’s 
equation  for  the  wetting  liquid  (w)  and  the  nonwetting  bubble  (nw)  along  the 
streamline  separating  the  two  fluids  (indicated  by  ‘0’  of  Figure  1) 

(w)  Pw  +  P^gx  +  =  Cw,  (2a) 

2 

(nw)  Pnw  +  Pn^gx  +  ^  Cnv/,  (2b) 
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Figure  2.  Oil  bubbles  in  water  (red,  top  of  the  figure)  and  water  bubbles  in  the  oil  (gree-blue, 
bottom).  In  the  smaller  diameters  (right)  instability  and  distortions  in  the  interface  are  smaller 
and  occur  further  downstream  from  the  bubble  tip. 
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Table  1.  Nondimensional  parameters  of  flow  in  closed  tubes:  Re: 
bubbles  velocities,  described  in  terms  of  Reynolds  number;  a* :  the 
dimensionless  surface-tension  parameter;  U*\  the  dimensionless 
velocity. 


System 

Diameter  (cm) 

Re 

(T* 

U* 

Glass-water-oil 

1.27 

296 

0.435 

0.214 

1.65 

726 

0.267 

0.354 

2.20 

1436 

0.145 

0.457 

2.54 

1890 

0.109 

0.483 

Acrylic-oil-water 

0.60 

0 

1.947 

0.000 

1.27 

84 

0.435 

0.246 

1.56 

169 

0.287 

0.363 

1.90 

231 

0.193 

0.368 

2.54 

426 

0.109 

0.441 

where  P  is  the  pressure,  q  is  the  tangential  liquid  velocity  along  the  surface  of  the 
bubble,  g  is  the  gravitational  acceleration,  p  is  the  density,  C  is  a  constant,  and  x 
is  measured  along  the  vertical  axis,  positive  upstream  with  zero  at  the  bubble  nose 
(Figure  1).  Subtracting  Equation  (2b)  from  Equation  (2a),  imposing  the  boundary 
conditions  at  the  bubble  nose  (i.e.,  x  -0a.ndq  =  0),  and  using  the  surface  tension, 
cr,  to  describe  the  pressure  differences  at  the  bubble  nose,  results  in 

^pgx  +  2(/?w9w  “  Pnvjqnw)  ~  o'{R  —  Rq  ),  (3) 

where  Ap  is  the  positive  density  difference  of  the  two  immiscible  liquids,  is 
the  curvature  of  the  slug  surface,  and  the  value  of  the  curvature  at  the  bubble 
nose.  (Note  that  we  used  Cnw  -  Cw  =  2a/Ro  in  Equation  (3),  resulting  from 
Equations  (2a)  and  (2b)  applied  at  the  stagnation  point,  x  =  q  =  0.)  As  long  as 
no  significant  circulation  is  created  within  the  liquid  bubble  by  the  outside  liquid 
then  Qnw  -  0  (we  could  just  as  easily  have  considered  the  other  limiting  case  when 
^  q^,  however  anticipating  our  experimental  results  we  found  that  Qnw  —  0, 
and  we  limit  our  theory  to  this  case).  By  introducing  the  coefficients  that  were 
already  calculated  by  Tung  and  Parlange  (1976),  Equation  (3)  can  be  written  in  the 
form 


Apga 


0.272  -  0.472 


a 

Apga?-  ’ 


(4) 


which  describes  the  velocity  at  which  a  bubble  rises  or  a  drop  falls.  Interestingly  it 
shows  exactly  the  same  nondimensional  parameters  as  used  previously  by  Zukoski 
(1966),  based  on  his  experimental  observations. 
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3.  Materials  and  Methods 

Tests  were  carried  out  with  two  immiscible  liquids,  water  and  oil  (Soltrol  220),  and 
two  types  of  tubes  (glass  and  acrylic).  Tbbes  were  1 00  cm  long  and  varied  in  inside 
diameter  from  0.6  to  2.54  cm  (Table  I). 

Initial  experiments  showed  that  in  an  oil-water  system,  the  water  wetted  pref¬ 
erentially  the  glass  and  the  Soltrol  wetted  preferentially  the  acrylic.  Bubbles  and 
drops  were  formed  by  a  sudden  overturn  of  the  tubes.  In  some  cases  flexible  tubes 
with  shutoff  valves  were  used  to  store  the  heavier  liquid  at  the  top  or  the  light  liquid 
on  the  bottom,  and  bubbles  or  drops  were  initiated  after  the  valve  was  opened.  To 
view  the  bubbles/drops  the  wetting  fluid  was  kept  transparent  and  the  nonwetting 
fluid  was  colored.  Blue  dye  FD&C  #1  was  used  for  water  in  the  acrylic  and  Sudan 
IV  (red)  for  Soltrol  in  the  glass  tubes.  Most  consistent  bubbles  and  drops  were 
formed  when  tubes  were  partly  filled  (80-90%)  by  the  wetting  liquid  first,  and 
completed  with  the  nonwetting  liquid.  If  the  nonwetting  liquid  occupied  most  of 
the  tube  the  wetting  liquid  moved  erratically  along  the  walls,  probably  because  the 
non  wetting  liquid  left  some  smear  along  the  walls. 

Bubble/drop  velocity  was  obtained  by  measuring  the  time  required  for  the 
bubble/drop  to  move  30  cm.  Pictures  of  the  bubbles/drops  in  different  diameter 
tubes  were  taken  (Figure  2)  in  order  to  obtain  the  shapes  of  the  bubbles  and  drops. 
The  pictures  were  scanned  into  a  graphic  software  and  were  scaled  to  the  same  size 
by  overlapping  the  lines  of  the  tube  walls. 

Soltrol-water  surface  tension  was  32.84 dyne  cm”  ^  measured  with  a  Model  21 
Fisher  Surface  Tensiomate.  The  Soltrol  viscosity  was  3.64  centipoise  [grcm”^ 
sec”*]  based  on  data  supplied  by  Phillips  Chemicals  Company  data. 


4.  Results  and  Discussion 

Bubbles/drops  formed  in  the  glass  and  acrylic  tubes  with  different  diameters  are 
shown  in  Figure  2.  In  all  cases,  the  wetting  fluid  flowed  along  the  wall  and  the 
nonwetting  fluid  in  the  center.  Thus  the  oil  moved  upward  as  a  bubble  in  the  glass 
tube  (Figures  2a,  b,  c)  and  the  water  moved  downward  as  a  drop  in  the  acrylic 
tube  (Figures  2d,  e,  f).  The  experimental  bubble/drop  nose  shapes,  compared  to  the 
predicted  shape  (Figure  1),  indicate  that  the  zero  stream  line  prediction  (Equation 
(1))  is  accurate  for  r  <  a/2,  and  diverges  increasingly  from  the  actual  shape 
as  r  increases,  in  agreement  with  Tung  and  Parlange  (1976).  We  also  found,  in 
agreement  with  Zukoski  (1966),  that:  (a)  bubble  length  and  initiation  method  did 
not  affect  the  bubble  velocity,  which  was  measured  with  standard  deviations  of  less 
than  5%  (see  also  Figure  3),  and  (b)  an  instability  in  the  two  liquids  interface  was 
observed  mainly  for  the  large  diameter  tubes  (Figure  2,  left)  for  both  the  oil  bubble 
in  water  and  the  water  drop  in  the  oil.  The  instability  did  not  affect  the  velocity  as 
long  as  large-amplitude  waves  were  restricted  to  positions  at  least  one  tube  radius 
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Figure  3.  Dimensionless  velocities,  U* ,  and  standard  deviations  as  a  function  of  the  dimension¬ 
less  surface-tension  parameter,  a* ,  compared  with  the  predicted  velocity  model  (Equation  (3)). 


downstream  of  the  bubble  nose.  With  the  smaller  diameter  the  distortions  were 
smaller  and  occurred  further  downstream  from  the  tip. 

The  bubble/drop  velocities,  described  in  terms  of  Reynolds  number.  Re,  the 
dimensionless  velocity,  U*,  and  the  dimensionless  surface-tension  parameter,  a*, 
are  given  by 


Re  = 


Mw 


U*  =  U 


-(1/2) 


a 


<7 

Apga^  ’ 


(5) 


and  shown  in  Table  I.  U*  was  plotted  as  a  function  of  cr*  and  compared  with 
the  predicted  velocity  model  (Equation  (4))  in  Figure  3.  The  good  agreement 
between  the  predicted  and  measured  results  indicates  that  (a)  the  wall-liquids 
wetting  relationship  determines  whether  a  bubble  or  a  drop  will  be  formed,  but 
while  formed,  it  does  not  affect  its  velocity,  and  (b)  the  assumption  of  no  movement 
inside  the  bubble/drop  is  indeed  the  preferred  assumption  for  modeling  purposes. 
This  assumption  was  also  confirmed  by  tiny  [<  0.25  mm]  dust  particles  floating 
inside  the  bubble/drop,  allowing  us  to  observe  that  no  significant  circulation  is 
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created  within  the  liquid  bubble  by  the  outside  liquid.  The  model  predicted  no 
movement  of  the  bubble  for  cr*  =  0.57  which  corresponds  to  a  =  0.55  cm  for  our 
liquids,  and  indeed,  no  movement  of  the  water  bubble  in  the  oil  was  observed  in 
tubes  with  a  <  0.55  cm. 

There  are  many  other  experiments,  mainly  of  gas  bubbles  in  liquid,  which  could 
be  added  to  Figure  3  (Zukoski,  1 966;  Fabre  and  Line,  1 992).  However,  none  involve 
oil/water,  which  is  of  interest  here,  and  the  case  of  gas  bubbles  is  already  known  to 
satisly  Equation  (4)  for  that  particular  case  with  the  same  accuracy  as  shown  in  the 
present  figure  (Tung  and  Parlange,  1976).  Thus,  these  earlier  experiments  would 
add  nothing  useful. 


5.  Applications  to  Porous  Media 

It  is  of  interest  to  investigate  how  the  model  for  single  tubes  applies  to  a  porous 
media  and  provides  some  qualitative  insight  to  the  observed  flow  patterns.  As  in  the 
single  tube  experiments  we  performed  a  series  of  experiments  with  two  immiscible 
fluids,  and  hydrophobic  and  hydrophilic  solid  phases. 

The  porous  medium  consisted  of  0.8-1. 2 mm  glass  beads.  Glass  is  hydrophilic. 
The  glass  beads  were  made  hydrophobic  by  spraying  them  with  Silicon  lubricant 
spray.  The  beads  were  packed  in  two  two-dimensional  chambers  (27  x  27  x  1  cm). 
One  chamber  had  glass  walls  and  held  the  glass  beads.  The  other  chamber  was 
made  of  acrylic  and  contained  the  hydrophobic  beads. 

Each  chamber  was  saturated  with  either  oil  or  water  (always  noncolored).  After 
placing  the  chamber  horizontally,  the  invading  liquid  was  injected  in  the  center  of 
the  chamber  at  a  rate  of  1  cm^/min.  The  invading  liquid  was  colored  with  FD&C 
blue  #1  for  the  water  and  Sudan  red  IV  for  the  Soltrol.  The  resident  fluid  leaves 
the  chamber  through  six  outlets  at  the  edge  of  the  chamber.  Photographs  of  the 
invading  liquid  were  taken  at  intervals  ranging  from  2  to  5  min. 

The  two  cases  where  water  displaced  the  oil  are  presented  in  Figures  4a  and  4b. 
In  Figure  4a,  the  medium  was  hydrophilic  resulting  in  a  smooth  front  for  the 
intruding  water.  In  Figure  4b,  the  medium  was  hydrophobic  and  the  water  which  is 
nonwetting  forms  a  finger-like  front.  Oil  displacing  water  in  a  hydrophobic  medium 
resulted  in  a  similar  pattern  as  the  smooth  front  in  Figure  4a  and  oil  displacing 
water  in  a  hydrophilic  media  gave  a  finger  structure  as  shown  in  Figure  4b.  When 
water  was  injected  in  Soltrol  in  the  hydrophilic  medium  only  15%  of  the  pore  space 
was  filled  with  oil.  When  water  was  injected  into  oil  in  the  hydrophobic  medium  at 
least  90%  was  filled  with  oil.  Clearly  a  continuous  smooth  invading  front  is  more 
efficient  at  removing  the  resident  fluid. 

If  we  calculate  the  tube  diameter  for  which  the  flow  is  zero,  from  Equation  (4), 
we  find  that  this  is  approximately  5  mm.  The  smallest  finger  width  observed  is 
exactly  of  the  same  order.  This  is  well  above  the  pore  size  of  our  porous  media 
with  the  Soltrol-water  system.  Thus,  nonwetting  fluid  in  different  pores  reconnect 
to  form  pathways  of  a  size  of  around  6  mm  or  more,  i.e.,  independently  of  the  pore 
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Figure  4.  Prints  of  the  colored  fractions  formed  by  injecting  a  constant  flux  of  colored  water  to 
the  oil  saturated  chamber,  using  two  types  of  porous  medium:  (a)  hydrophilic,  (b)  hydrophobic. 
The  numbers  indicate  the  time  (min)  elapsed  from  the  beginning  of  the  injection. 


size  distribution.  However,  it  is  also  clear  that  the  exact  shape  and  dimension  of 
the  invading  fingers  is  affected  by  the  pore  distribution,  so  that  the  6  mm  prediction 
is  only  a  semiquantitative  prediction,  which  should  be  very  useful,  for  instance,  to 
estimate  the  surface  of  contact  between  the  two  fluids  when  the  intruding  fluid  is 
nonwetting. 

After  lifting  the  chamber  the  fluids  always  moved  readily.  Fluid  movement 
eventually  stopped  after  it  had  reached  the  configuration  similar  to  the  one  shown 
in  Figure  5.  Observe  that  the  width,  i.e.,  the  smaller  characteristic  dimension  of  all 
ganglia  is  less  than  6  mm  (only  two  near  the  top  of  Figure  5  are  even  close  to  6  mm). 
Both  results  are  consistent  with  flow  in  a  single  tube.  During  imbibition  the  flow 
requirement  imposes  that  the  pathway  be  larger  than  6  mm.  After  drainage  stops  the 
dimension  of  the  ganglia  in  the  direction  perpendicular  to  the  flow  must  be  less  than 
6  mm.  As  drainage  proceeds  the  exact  size  of  the  ganglia  left  behind  will  depend  on 
the  connections  between  the  entrapped  fluid  and  the  main  pathway.  This  point  was 
made  before,  e.g.,  see  the  excellent  discussion  by  Powers  et  al.  (1992).  Here  we 
suggest  in  addition  that  ganglia  have  a  maximum  size.  It  is  interesting  that  ganglia 
formation  and  movement  has  been  discussed  in  the  past  (e.g.,  see  Payatakes  and 
Dias  (1984)  or  Adler  and  Brenner  (1988))  but  not  in  relation  to  pathway  dimensions 
as  discussed  here. 

Thus  it  seems  that  the  minimum  size  for  movement  in  the  single  tube  experiment 
provides  some  indication  of  the  order  of  magnitude  of  the  width  of  pathways  and 
maximum  size  ganglia  in  a  porous  medium  left  by  drainage  and  minimum  finger  size 
during  imbibition,  which  will  normally  be  greater  than  the  pore  size.  Experiments 
are  currently  underway  to  examine  the  effects  of  surfactants  on  this  ganglion  and 
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Figure  5.  Trapped  water  in  the  hydrophobic  porous  material.  Picture  was  taif^n  after  the 
chamber  from  Figure  lb  was  placed  in  a  vertical  position  and  allowed  to  drain. 


pathway  width,  to  check  whether  Equation  (4)  still  holds,  and  quantify  the  effect 
of  surfactant  in  promoting  soil  remediation. 

6.  Conclusions 

Based  on  Ae  study  of  wetting  and  nonwetting  liquids  motion  in  tubes,  we  expect 
more  efficient  displacement  when  the  wetting  liquid  displaces  the  non  wetting  liquid 
in  a  porous  medium.  Smooth  wetting  fronts  are  expected  if  the  displacing  liquid 
is  the  preferentially  wetting  liquid.  Finger-like  shapes  with  a  width  larger  than  a 
characteristic  length  are  observed  if  the  displacing  liquid  is  nonwetting.  Finally, 
lumps  of  the  nonwetting  fluid  having  a  width  smaller  than  the  same  characteristic 
length  can  be  readily  entrapped  in  the  porous  medium  when  gravity  alone  is  the 
cause  of  motion.  The  characteristic  length  in  both  cases  is  given  by  [a 
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7  Abstract 

8 

9  Predicting  the  occurrence  of  preferential  flow  is  crucial  as  preferential  flow  causes  faster 

10  downward  movement  of  water  and  solutes,  which  may  result  in  groundwater  contamination.  The 

1 1  effects  of  water  repellency  on  preferential  flow  were  studied  by  varying  the  silane  coating  on  a 

12  silica  sand.  Uniform  rain  on  the  non-coated  substrate  resulted  in  a  stable  wetting  front  and  no 

13  preferential  flow,  while  the  coated  sands  produced  an  unstable  wetting  front  and  preferential  flow. 

14  The  occurrence  of  preferential  flow  was  directly  linked  with  the  soil’s  water  repellency,  and  its 

15  pressure-saturation  relationship.  In  addition,  the  front  velocity  and  the  water  entry  pressure 

16  increased  linearly  with  a  higher  degree  of  water  repellency  of  the  sand. 


2 


1  Introduction 

2 

3  Soils  which  have  hydrophobic  properties  (also  called  water  repellent  soils),  can  resist  or  retard 

4  surface  water  infiltration  (Brandt,  1969).  This  affects  plant  growth  and  surface  water  distribution. 

5  These  soils  can  be  found  in  Florida,  Arizona  and  California  (DeBano,  1969),  as  well  as  abroad,  in 

6  New  Zealand,  England,  Australia  (Bond,  1969)  and  the  Netherlands  (Hendrickx  et  al.,  1988; 

7  Dekker  and  Ritsema,  1994).  Water  repellent  soils  may  be  more  prevalent  than  is  commonly 

8  acknowledged  (Bond,  1964;  Gilmour,1968;  DeBano,  1969;  Marshall  and  Holmes,  1979;  Nakaya, 

9  1982;  Dekker  and  Ritsema,  1996).  Some  researchers  suggest  that  repellency  is  the  norm  rather 

10  than  the  exception  (Dekker,  1988;  Dekker  and  Jungerius,  1990;  Wallis  and  Horne,  1992).  The 

1 1  degree  of  repellency  may  vary  between  soils. 

12 

13  Besides  the  retardation  or  resistance  of  surface  water  infiltration,  water  repellent  soils  have  been 

14  associated  with  preferential  flow  paths  (Jamison,  1945;  Bond,  1964;  Gilmour,  1968).  Preferential 

15  flow  paths  create  spatial  variability  in  soil  moisture  affecting  plant  growth  and  surface  distribution 

16  (Dekker  and  Ritsema,  1994).  In  addition,  preferential  flow  allows  much  faster  transport  of  water 

17  and  solutes,  thus  creating  greater  potential  for  groundwater  contamination. 

18 

19  In  this  paper,  the  effects  of  water  repellency  on  preferential  flow  by  varying  the  silane  coating  on  a 

20  silica  sand  is  systematically  studied.  Uniform  rain  on  the  non-coated  substrate  resulted  in  a  stable 

21  wetting  fi'ont  and  no  preferential  flow,  while  the  coated  sands  produced  an  unstable  wetting  front 

22  and  preferential  flow.  A  direct  link  between  the  soil’s  water  repellency,  its  pressure-saturation 

23  relationship,  and  the  occurrence  of  preferential  flow  was  found. 
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1  Materials  and  Methods 

2 

3  Blasting  silica  sand  (W.F.  Saunders  &  Sons,  Inc.)  was  used  as  the  substrate  to  analyze  the  effect 

4  of  repellency.  Table  1  shows  the  screen  analysis  for  this  sand.  Different  degrees  of  water 

5  repellency  were  added  to  the  sand  through  the  following  procedure.  An  ethanol  solution 

6  containing  5%  octadecyltrichlorasilane  (100  ml  OTS  (Fisher  Scientific)  in  2  liter  ethanol)  was 

7  mixed  for  5  hours  with  25  kg  of  blasting  sand  which  produced  an  extremely  water  repellent  sand 

8  after  drying  (Bradford  and  Leij,  1996).  A  basis  sand  batch  was  made  by  mixing  3  kg  of  the 

9  extremely  repellent  sand  with  45  kg  regular  non-repellent  sand.  This  basis  sand  was  than  diluted 

10  with  regular  sand  to  achieve  different  degrees  of  repellencies  (Table  2).  The  most  repellent  sand 

11  tested  was  achieved  by  mixing  the  undiluted  repellent  sand  with  10  kg  of  regular  sand.  In  Table 

12  2,  the  first  number  in  the  notation  for  the  batches  refers  to  the  amount,  in  kg,  used  of  the  basis 

13  repellent  sand  (except  for  the  most  repellent  sand)  and  the  second  number  refers  to  the  amount,  in 

14  kg,  of  the  regular  sand.  The  degree  of  repellency  defines  the  amount  of  highly  repellent  sand  (the 

15  sand  that  was  mixed  with  OTS  in  an  ethanol  solution)  versus  the  total  amount  of  sand  used. 

16 

17  Besides  the  ratio  of  repellent  to  non-repellent  sand,  the  water  drop  penetration  time  (WDPT)  was 

18  used  to  quantify  the  degree  of  repellency.  The  WDPT  test  consists  of  randomly  applying  a  single 

19  water  drop  (0.05  ml)  onto  the  sand  surface  and  measuring  the  amount  of  time  (in  seconds)  it  takes 

20  to  infiltrate  the  soil  (Letey,  1969).  The  WDPT  for  the  different  batches  is  shown  in  Table  2 

21  (Letey,  1969;  King,  1981.) 
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1  A  polycarbonate  slab  chamber  was  built  to  observe  and  quantify  the  preferential  flow  phenomena 

2  (see  Figure  1).  The  interior  dimensions  of  the  chamber  were  42  cm  wide,  61  cm  tall  and  0.8  cm 

3  thick  with  1  cm  polycarbonate  walls.  Seven  high  speed  tensiometers,  to  measure  water  pressures, 

4  were  threaded  into  the  chamber  at  a  height  of  32  cm  and  6  cm  horizontal  spacing.  Each 

5  tensiometer  consisted  of  a  flat  stainless  porous  plate  (0.8  cm  diameter,  2  micrometer  pore  size) 

6  sealed  in  a  hollow  brass  housing  (Selker  et  al.,  1992a).  The  tensiometers  were  inserted  so  the 

7  porous  plate  was  flush  with  the  inside  of  the  chamber,  to  minimize  their  effect  on  packing  and 

8  fluid  flow.  Each  tensiometer  was  filled  with  de-gassed  water  and  attached  through  a  stiff  nylon 

9  tube  and  a  three-way  valve  to  a  pressure  transducer  and  a  vertical  water  tube.  By  turning  the 

10  valve  the  transducer  could  be  connected  to  the  vertical  tube  for  calibration,  or  to  the  tensiometer 

11  for  measurement.  The  output  of  the  transducers  was  read  by  a  data  acquisition  system,  consisting 

12  of  a  personal  computer  with  a  DAS-800  card  (Keithley  Metrabyte),  and  a  multiplexer  to  process 

13  in  total  16  incoming  channels. 

14 

15  The  chamber  was  filled  under  a  hopper  with  three  randomizing  screens  (Glass  et  al.,  1989a). 

16  Before  every  filling,  an  extension  of  15  cm  was  put  on  top  of  the  chamber  to  remove  the  always 

17  poorly  packed  top  section  and  to  also  achieve  an  homogenous  packing  in  the  chamber. 

18 

19  Using  this  chamber,  the  soil  characteristic  curves  were  determined  for  all  the  sands.  The  drying 

20  curves  were  obtained  through  the  following  procedure.  A  peristaltic  pump  pumped  the  water  in  at 

21  0.4  cm/min  until  the  water  reached  a  height  of  45  cm.  After  pausing  for  30  minutes,  the  chamber 

22  was  drained  by  pumping  the  water  out  at  0.4  cm/min  until  the  water  fringe  reached  a  height  of  7 

23  cm.  After  24  hours,  the  front  panel  of  the  chamber  was  removed  and  the  sand  column  was 
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1  segmented  in  1  cm  levels.  Each  segment  was  weighed,  dried  for  24  hours  at  105°  C,  and  weighed 

2  again.  This  procedure  was  practiced  on  the  left  and  right  side  of  the  chamber,  which  resulted  in  2 

3  drying  curves  per  experiment,  those  two  where  then  averaged  to  obtain  the  drying  curve.  The 

4  wetting  curves  were  also  determined  in  the  chamber;  a  constant  head  was  connected  to  the 

5  bottom  of  the  chamber  and  left  connected  for  24  hours  till  equilibrium  was  reached.  The  chamber 

6  was  again  taken  apart  to  section  the  sand  and  determine  the  water  saturation. 

7 

8  For  infiltrations,  the  packed  chamber  was  placed  in  front  of  a  fluorescent  light  bank  and  the 

9  tensiometers  were  attached  to  the  data  acquisition  system.  As  wet  sand  transmits  more  light  than 

10  dry  sand  (Glass  et  al.,  1989b),  light  transmission  was  used  to  observe  the  position  and  velocity  of 

11  the  wetting  front.  The  experiments  were  recorded  by  a  black  and  v  h;  ra  (Panasonic,  TV 

12  Camera  WV  1850)  and  a  video  cassette  recorder  (Panasonic  AG  6720).  The  recorded 

13  experiments  were  than  digitized  with  a  frame  grabber  from  Data  Translaticiis  .iiid  analyzed  with 

14  Global  Labs  Imaging  software. 

15 

16  The  water  application  system  consisted  of  a  peristaltic  pump  driven  point  source  which  was 

17  translated  back  and  forth  across  the  soil  surface  by  a  rotating  cam  shaft.  To  eliminate  the  possible 

18  flow  down  the  side  edges  of  the  chamber,  gutters  were  installed  to  prevent  rain  from  impacting 

19  within  2  cm  of  the  chamber  edges.  Infiltrations  were  conducted  with  an  irrigation  rate  of  0. 16 

20  cm/min.  Experiments  with  infiltration  rates  between  0.058  cm/min  and  0. 19  cm/min  all  produced 

21  similar  results. 
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1  Results 

2 

3  Figure  2  shows  the  averaged  drying  soil  characteristic  (or  pressure-saturation)  curves  for  the  five 

4  different  sands.  The  curves  were  very  similar  over  most  of  the  range  of  water  contents  with 

5  statistical  differences  ( curves  more  than  2  standard  deviations  apart )  only  for  water  saturations 

6  between  10-40  %.  In  this  region,  the  more  repellent  the  sand  the  less  water  it  holds  at  a  particular 

7  negative  pressure. 

8 

9  Figure  3  shows  the  averaged  wetting  soil  characteristic  (or  pressure-saturation)  curves  for  the  five 

10  different  sands.  As  opposed  to  the  drying  curves,  the  wetting  curves  varied  much  more  with  the 

11  degree  of  water  repellency,  with  two  types  of  differences  between  the  curves.  First,  the  water 

12  entry  pressures  increased  with  the  degree  of  repellency.  Second,  the  slope  of  the  repellent  curves 

13  was  flatter  than  that  from  the  pure  sand.  For  the  repellent  sands  the  actual  wetting  slope  may  be 

14  slightly  sharper  than  that  measured.  By  sectioning  the  chamber  every  centimeter,  we  may  have 

15  taken  a  chunk  of  wet  sand  and  a  chunk  of  dry  sand  at  the  same  time.  Thus,  an  abrupt  wetting 

16  front  in  the  water  repellent  sands  will  still  appear  as  having  a  minimum  slope. 

17 

18  Major  differences  between  the  different  water  repellent  sands  were  observed  when  water  was 

19  infiltrated  into  the  sands.  For  the  non-water  repellent  sand,  the  water  moved  downv-’ard  in  a 

20  uniform  stable  front,  while  for  the  repellent  sands,  the  water  moved  downward  in  preferential  flow 

21  paths.  In  the  higher  repellent  sand  the  water  ponded  before  it  moved  downwards  in  the  fingers. 

22  To  gain  insight  into  this  phenomena,  we  measured  front  velocity,  front  width  (for  finger  flow), 

23  and  water  pressures  within  the  front. 
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1 

2  Figure  4  shows  the  time  series  of  the  pressure  at  the  tensiometers  as  the  wetting  front  passed  in 

3  the  non-water  repellent  sand.  The  front  hit  the  tensiometer  at  T=3259  sec,  and  the  pressure  at  the 

4  front  reached  -11.8  cm.  Pressures  before  the  front  were  meaningless  as  the  sand  is  air  dry  and 

5  there  is  no  connection  to  the  water  in  the  soil.  Behind  the  front  the  pressure  remained  -11  cm. 

6  Once  the  infiltration  was  stopped,  the  pressure  decreased  till  -26  cm  and  gradually  went  till  -30 

7  cm. 

8 

9  Figure  5  shows  the  equivalent  time  series  of  the  pressure  in  the  5.6  %  water  repellent  sand.  The 

10  front  hit  the  tensiometer  at  T=  2199  sec,  and  the  pressure  at  the  front  reached  3.6  cm. 

11  Immediately  behind  the  front,  the  pressure  quickly  began  decreasing  before  settling  at  -12  cm. 

12  Once  the  infiltration  was  stopped,  the  pressure  further  decreased  till  -27  cm  and  gradually  moved 

13  till  -30  cm.  This  post  infiltration  pressure  decrease  drops  at  a  rate  which  was  much  higher  than  in 

14  regular  sand. 

15 

16  Figure  6  depicts  the  measured  water  pressure  at  the  wetting  front  (error  bars  i'-dicatc  tire  standard 

17  deviations)  versus  %  repellency  of  the  sand.  The  water  entry  values  measured  from  the  static 

18  wetting  curves  are  also  shown.  These  pressures  are  seen  to  be  closely  related,  and  both  clearly 

19  increased  with  the  degree  of  repellency.  This  increase  was  roughly  linear,  as  a  linear  regression 

20  ■  resulted  in  =  91  and  97%  for  the  wetting  front  pressures  and  water  entry  values,  respectively. 

21  A  positive  water  entry  pressure  was  required  when  the  degree  of  repellency  rose  above  5%,  which 

22  implied  a  ponding  layer  on  top  of  the  sand  grains.  This  ponding  was  also  related  to  the  WDPT 

23  test.  The  WDPT  for  those  5%  repellent  sands  was  >  600  sec,  consequently  there  was  positive 
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1  entry  pressure  necessary  for  infiltration,  thus  we  might  expect  that  a  water  droplet  of  a  couple  of 

2  mm  high  would  never  penetrate  the  repellent  sand  grains. 

3 

4  Figure  7  gives  the  finger  velocity  versus  %  repellency  and  the  fitted  linear  regression  line  (error 

5  bars  indicate  the  standard  deviation  of  the  measured  points).  The  finger  velocity  was  measured, 

6  starting  5  cm  below  the  sand  surface  and  stopped  5  cm  before  reaching  the  bottom.  The  finger 

7  velocity  increased  linearly  with  the  degree  of  repellency  (  =  98%  ) 

8 

9  Figure  8  depicts  the  observed  front  width  at  the  wetting  front  versus  %  repellency.  The  error  bars 

10  indicate  the  standard  deviations.  The  finger  width  was  measured  in  5  cm  intervals  starting  at  the 

11  distribution  layerBecause  the  wetting  front  for  the  non-repellent  sand  is  uniform,  the  width  was 

12  not  depicted.  Widths  at  the  wetting  front  for  the  various  degrees  of  repellencies,  except  the 

13  highest  one;  were  found  to  be  around  3  cm.  The  width  for  the  highest  degree  of  repellency  was 

14  5.8  cm  or  nearly  double  the  width  compared  to  the  other  front  widths. 

15 

16  Discussion 

17 

18  In  this  section,  we  will  discuss  how  different  degrees  of  repellency  affected  the  wetting  and  drying 

19  pressure-saturation  curves  and,  in  turn,  how  these  curves  can  be  related  to  the  observed 

20  infiltration  properties.  Ideally,  this  should  give  us  an  insight  in  the  likelihood  of  preferential  flow 

21  in  soils  from  the  observed  degree  of  wettability,  or  non- wettability. 

22 

23  Repellency  affects  the  wetting  pressure-saturation  curves  of  the  soil  mucli  more  than  the  drying 
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1  pressure-saturation  curves  (Figure's  2  &  3).  The  similarities  in  the  drying  curves  imply  that  once 

2  the  soil  grains  are  wetted,  the  soil  acts  like  a  normal  non-repellent  soil.  Only  when  the  soil  is 

3  drained  to  less  than  40  %  saturation  does  the  water  start  leaving  some  portions  of  the  soil  drier 

4  and,  thus,  the  repellency  starts  affecting  the  water  retention.  Conversely,  on  the  initial  wetting,  the 

5  soil  begins  dry  and  the  water  repellency  reaches  its  full  effect.  Thus,  the  greater  the  water 

6  repellency  the  greater  the  hysteresis  in  the  soil’s  pressure-saturation  relationship. 

7 

8  Importantly,  we  find  that  these  drying  and  wetting  pressure-saturation  relationships  can  be  related 

9  to  the  fingering  properties.  First  of  all.  Figure  3  and  Table  3  show  that  the  water  entry  pressures 

10  measured  from  the  wetting  curve  match  the  water  pressure  observed  from  the  tensiometers 

1 1  directly  behind  the  infiltrating  front.  The  behavior  of  the  infiltrations  differ  from  this  point, 

12  depending  on  the  repellency  of  the  soil.  For  the  non-water  repellent  soil,  the  wetting  front  is 

13  uniform  and  stable,  and  Figure  4  shows  that  the  pressure  behind  the  front  remains  the  initial 

14  wetting  pressure.  This  implies  that  the  water  at  all  positions  behind  the  front  is  driven  by  a  unit 

15  gradient,  and  that  the  conductivity  equals  the  downward  flux.  Thus,  in  the  case  of  an  uniform 

16  wetting  fi'ont,  the  soil  wets  to  the  water  content  such  that  the  unsaturated  conductivity  K(9) 

17  equals  the  downward  flux. 

18 

19  For  the  water  repellent  soil,  the  wetting  front  is  unstable,  and  Figure  5  shows  that  the  pressure  at 

20  the  tip  is  quickly  dropping  from  the  initial  wetting  pressure.  This  drop  in  pressure  for  unstable 

21  fi'onts  was  first  elucidated  by  Raats  (1973).  The  final  pressure  behind  the  wetting  front  is  seen  to 

22  be  roughly  the  pressure  at  the  knee  in  the  drying  curve.  Thus,  in  the  case  of  the  unstable  wetting 

23  front,  the  soil  dries  to  the  water  content  such  that  the  unsaturated  conduc  tivity  K(0)  equals  the 
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1  downward  flux.  Since  the  flow  is  through  preferential  paths,  the  flux  in  each  wetting  path  is 

2  greater  than  the  overall  applied  flux.  Interestingly,  the  water  pressure  behind  the  unstable  front  is 

3  less  than  the  water  pressure  behind  the  stable  front,  even  though  the  soil  behind  the  unstable  front 

4  is  exhibiting  a  higher  conductivity  (and  is  likely  to  be  more  saturated)  than  the  soil  behind  the 

5  stable  front. 

6 

7  Clearly,  the  pressures  inside  the  infiltrating  fronts  can  be  directly  related  to  the  measured  pressure- 

8  saturation  curves  for  each  soil.  These  curves  can  also  be  related  to  the  stability  of  the  wetting 

9  front,  the  velocity  of  the  fingers,  and  the  width  of  the  wetting  front. 

10 

1 1  The  finger  velocity  increases  linearly  with  the  degree  of  repellency  of  the  sand.  This  is,  at  first, 

12  contradictory;  the  higher  the  repellency,  the  more  the  sand  resists  the  infiltration  of  water.  A 

13  simple  argument  on  how  repellency  enhances  fingering  is  given  as  follows.  Previously,  Weitz  et  al. 

14  (1987)  have  found  that  the  velocity  of  the  fluid/fluid  interface  (the  water  front  in  this  case) 

15  increases  with  the  pressure  drop  at  the  interface.  For  our  case,  as  the  sand  is  more  repellent,  a 

16  higher  water  pressure  must  be  reached  before  water  can  enter  it.  But  once  vater  enters  the  sand 

17  and  coats  the  grains  all  the  sands  have  more  or  less  the  same  drying  properties.  The  flow  then 

18  proceeds  normally  (same  conductivity  and  P-S  relationship)  and,  thus,  the  major  pressure  drop  is 

19  taken  up  at  the  interface  and  the  front  velocity  is  higher.  Finger  velocity  is  important,  as  for  a 

20  constant  flux,  the  greater  the  velocity  the  fewer  the  fingers.  Also,  Selker  et  al.  (1992b)  have 

21  shown  that  the  saturation  profile  within  the  finger  can  be  related  to  the  finger  velocity  and  soil 

22  conductivity. 


11 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 


The  stability  of  the  wetting  front  for  different  repellencies  can  be  understood  heuristically  as 
follows.  Since  Saflfman  and  Taylor  (1958),  it  has  been  known  that  for  gravity-driven 
displacements,  the  front  is  unstable  whenever  the  flux  is  less  than  the  soil’s  unsaturated 
conductivity.; 


q  <  K(0) 


(1) 


When  the  wetting  pressure-saturation  curve  has  little  or  no  slope  (the  repellent  sands),  the  wetting 
front  is  very  abrupt  and  directly  behind  it  the  saturation  is  very  high.  This  leads  to  a  high 
conductivity  and  the  instability  criterion  is  satisfied.  However,  when  the  wetting  pressure- 
saturation  curve  has  an  appreciable  slope  (the  non-repellent  sand),  the  wetting  front  is  not  as 
abrupt  and  directly  behind  it  there  is  a  range  of  saturations.  Thus  we  observe  a  situation  where  the 
soil  wets  only  to  a  saturation  such  that  the  conductivity  equals  the  flux,  and  the  instability  criterion 
is  not  satisfied. 


The  slope  of  the  wetting  curve  has  been  more  quantitatively  related  to  the  finger  width  (and,  in 
turn,  front  stability)  by  Liu  et  al.  (1994)  through: 


d  = 
6f  = 


ti  +  1.5  ’ 


(2) 
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where; 
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the  slope  of  the  soil  moisture  characteristic  curve 

ri  =  is  the  exponent  in  the  Brooks  and  Corey  (1964)  conductivity  saturation 
relationship: 

ri  =  f  +  3  (3) 

with  X  a  positive  index  related  to  the  pore  size  distribution  of  the  soil. 

This  predicts  that  the  flatter  the  main  wetting  curve,  the  smaller  the  finger  and  the  increased 
likelihood  of  instability.  This  trend  is  replicated  in  our  data.  Table  4  shows  the  measured  front 
width  versus  the  predicted  front  width  (with  n=  3.3  and  9f  =  100).  These  predictions  do  not 
match  the  measured  finger  widths  very  well.  This  might  be  due  to  the  extreme  small  slope  of  the 
wetting  curve.  Aa  mentioned  earlier,  the  wetting  curves  were  determined  by  sectioning  every 
centimeter,  which  might  be  large  to  accurately  define  the  wetting  front  Nevertheless,  the 
following  trend  can  still  be  deducted:  repellent  sands,  which  finger,  have  wetting  curves  much 
flatter  than  the  non-repellent  sands,  which  do  not  finger.  The  large  finger  width  for  the  most 
repellent  sand  contradicts  the  trend,  however,  the  observed  finger  was  possibly  two  smaller 
fingers  traveling  together  as  the  wet  portion  of  the  tip  was  observed  to  be  splitting  into  two. 
However,  the  finger  reached  the  bottom  before  anything  further  could  be  noticed. 

In  summary,  we  have  shown  that  by  making  identical  soils  water  repellent,  unstable  water  fronts 
and  preferential  flow  patterns  are  created  where  none  existed  before.  Water  repellency  directly 
affects  the  soil's  pressure-saturation  relationship,  and  these  changes  are  represented  in  the 
preferential  flow  characteristics.  Finger  velocity  and  water  entry  pressure  are  both  linearly  related 
to  the  water  repellency  of  the  sand,  but  the  finger  widths  show  a  less  definitive  trend.  This  implies 
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1  for  field  soils,  that  water  repellency,  or  lack  thereof,  plays  a  controlling  roll  in  the  creation  of 

2  preferential  flow,  and  groundwater  pollution  such  as  was  shown  already  experimentally  by 

3  Hendrickx  et  al.  (1993.) 
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1  Figure  captions 

2 

3  Figure  1 ;  A  schematic  of  the  experimental  chamber.  Water  was  infiltrated  from  a  cam-driven  point 

4  source  at  the  top.  Pressure  readings  were  made  by  the  row  of  seven  tensiometers. 

5 

6  Figure  2;  The  averaged  drying  curves  for  identical  sands  with  different  degrees  of  repellency.  The 

7  drying  curves  statistically  differ  only  between  10-40%  saturation.  The  arrow  indicates  the  trend 

8  for  the  more  repellent  sands. 

9 

10  Figure  3:  The  averaged  wetting  curves  for  identical  sands  with  different  degrees  of  repellency. 

11 

12  Figure  4:  A  time  series  of  the  measured  tensiometer  pressure  in  regular  sand. 

13 

14  Figure  5;  A  time  series  of  the  measured  tensiometer  pressure  in  5.6%  water  repellent  sand.  Notice 

15  the  sharp  drop  in  pressure  after  the  finger  hit  the  tensiometer,  this  is  the  typical  property  for 

16  unstable  wetting  fronts  (Raats,  1973.) 

17 

18  Figure  6;  The  water  entry  pressure  obtained  from  the  wetting  pressure-saturation  curves  and  the 

19  tensiometer  readings  directly  behind  the  front  versus  %  repellency.  The  error  bars  indicate  the 

20  standard  deviation  for  the  tensiometer  measurements.  The  line  is  a  linear  regression  on  the 

2 1  tensiometer  measurements  and  has  a  of  9 1  % . 
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1  Figure  7;  The  front  velocity  is  depicted  versus  the  %  water  repellency.  The  error  bars  indicate  the 

2  standard  deviations.  The  line  is  a  linear  regression  and  has  a  of  98%. 

3 

4  Figure  8;  The  finger  widths  versus  %  repellency.  The  error  bars  indicate  the  standard  deviations. 
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1  Table  1 :  The  screen  analysis  (by  weight)  of  the  blasting  silica  sand  provided  by 

2  manufacturer,  W.F.  Saunders  &  Sons  Inc. 


Screen  ( lines  /  inch  ) 

%  Retained 

Cumulated 

20 

Trace 

Trace 

30 

0.8 

0.8 

40 

14.1 

14.9 

50 

38.6 

53.5 

70 

28.4 

81.9 

100 

11.0 

92.9 

Through  100 

7.1 

100.0 

1 


1  Table  2:  The  WDPT  for  the  different  batches  (Letey,  1969;  and  King,  1981 .) 


Water  Repellent  Ratios  WDPT  (Water  Drop  Penetration  Time) 


Batch  type 

%  Water  Repellency 

WDPT 

Description 

Regular 

0 

<0.5 

Wettable  (<5  sec) 

10+10 

3.13% 

40 

Slightly  water 
repellent  (5-60  sec) 

10+2.5 

5.00% 

2400 

Severely  water 
repellent  (600-3600 
sec) 

6+0.6 

5.60% 

>3600 

Extremely  water 
repellent  (>3600sec) 

1+10 

9.00% 

>3600 

Extremely  water 
repellent  (>3600sec) 

2 


Table  3:  Water  entry  pressures  measured  with  tensiometer  and  water  entry  values  from  the 
wetting  curves. 


%  Water  Repellency 

Water  Entry  Pressure 
Data  From  The  Fingertip 
(cm) 

Water  Entry  Pressure 
Data  From  Wetting 
Curves 
(cm) 

data 

standard 

deviation 

Regular 

-9.3 

1.3 

-9.2 

3.13% 

-2.2 

0.7 

-3.0 

5.00% 

-1.5 

1.0 

-1.0 

5.60% 

3.6 

NA 

2.0 

9.00% 

5.6 

0.9 

5.0 

3 


1  Table  4:  Measured  versus  predicted  finger  widths. 

2 

%  Water  Repellency  Measured  Widths  (cm)  Predicted  Widths  (cm) 


Regular 

NA 

16.4 

3.13% 

2.1 

3.6 

5.00% 

2.8 

2.6 

5.60% 

2.7 

3.1 

9.00% 


7.1 


1.3 


0  20 


80 


%  Saturation 


0 


2 
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10  Abstract 

11 

12  Gravity  driven  preferential  flow  (fingering)  can  greatly  affect  how  one  fluid  displaces  another  in 

13  the  subsurface.  Using  synchrotron  x-rays  and  miniature  tensiometers,  we  have  studied  the  internal 

14  properties  of  these  preferential  flow  paths  when  water,  with  and  without  surfactants,  infiltrating 

15  an  oil  saturated  porous  media.  Water  pressure  and  saturation  decrease  behind  the  fi'ont,  similar  to 

16  fingers  in  air-dry  sand,  with  quantitative  differences  for  different  surfactants.  Using  a  formulation 

17  originally  described  for  water  infiltrating  air-dry  sand,  we  find  that  the  finger  widths,  pressure 

18  drops  within  the  fingers,  finger  tip  lengths,  and  finger  splitting  dynamics  can  all  be  related  to  the 

19  interfacial  tension  between  the  fluids. 

20 
21 

22  *  Corresponding  Author 


3 

1  Concurrently,  soil  scientists  have  been  interested  in  preferential  flow  during  infiltration  events, 

2  when  the  flow  rate  is  slow  and  the  fingering  is  driven  by  density  effects,  (hillel?)  In  the  case  of 

3  water  infiltrating  air  dry  soil,  Parlange  and  Hill  (1976)  accounted  for  the  capillary  effects  of  the 

4  porous  media  by  treating  the  wetting  front  as  diffuse.  The  resultant  expression  for  the  finger 

5  dimensions  includes  the  sorptivity  S  of  the  soil  which  is  a  function  of  the  unsaturated  conductivity 

6  and  the  capillary  pressure-saturation  relation.  Recently  this  has  been  extended  to  water  infiltrating 

7  oil-saturated  soils  by  Chandler  et  al.  (1997)  with  the  assumption  that  the  flow  is  slow  enough  to 

8  eliminate  viscous  effects. 

9 

10  Most  experimental  studies  have  focused  on  fingering  in  Hele-Shaw  cells  or  at  the  outside  wall  of 

11  thin  sands  due  to  the  difficulty  in  “seeing”  inside  a  complete  porous  media.  Viscous  fingering  has 

12  been  studied  as  a  function  of  front  velocity,  interfacial  tension,  and  wettability  in  Hele  Shaw  cells 

13  (Homsy  1987)  and  at  the  wall  of  thin  glass  bead  packs  (Stokes  et  al.,  1986).  Of  particular  interest 

14  is  the  use  of  surfactants  or  detergents  in  the  aqueous  phase,  which  are  used  for  enhanced  oil 

15  recovery  or  remediation.  By  altering  the  interfacial  tension  with  a  surfactant,  the  preferential  flow 

16  properties  change  dramatically.  There  has  been  a  wealth  of  work  in  the  very  low  interfacial 

17  tension  limit,  as  fingers  themselves  are  observed  to  be  unstable  and  split.  (Bensimon  1986,  Homsy 

18  1987) 

19 

20  Inside  real  porous  media  the  fluid  saturation  at  and  behind  the  interface  can  show  wide  variations 

21  and  affect  the  scale  of  the  preferential  flow.  (Glass  19??,  Selker  et  al.,  1993b).  In  this  paper,  we 

22  have  used  synchrotron  x-ray  attenuation  to  measure  accurately  and  quickly  fluid  saturations 
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during  fingering  within  a  two-dimensional  slab  of  a  three-dimensional  porous  media.  Along  with 
fast  measurements  of  the  fluids  pressures,  a  better  understanding  of  the  dynamics  of  fingering  can 
be  resolved.  By  using  sands  already  well  characterized  for  fingers  through  dry  sand,  we  find  that 
the  fingers  through  oil  are  qualitatively  similar  and  the  quantitative  differences  can  be  directly 
related  to  the  interfacial  tension  and  the  density  differences  of  the  fluids.  Thus,  the  same  physics  of 
unstable  fluid  displacements  can  be  used  to  describe  and  predict  gravity  driven  fingered  flow  for 
many  fluid  pairs. 


Theory 

We  wish  to  obtain  expressions  for  finger  width  and  the  fluid  energy  state  (pressure)  inside  the 
finger  for  the  case  of  water  with  a  surfactant  infiltrating  a  oil  saturated  column.  The  capillary 
pressure  (Pw/nw,  also  known  as  the  suction  in  soil  science),  which  is  the  pressure  difference 
between  the  non-wetting  (Pnw)  and  wetting  (P*)  phase,  is  related  to  the  pore  radius  (r)  through  the 
Laplace  equation. 


w/nw 


=  P  -  P  = 


w  /  nw  ^0^(4*  w/  nw) 


w 


(2) 


where  the  subscript  w/nw  refers  to  each  particular  wetting  and  non-wetting  fluid  pair,  a  is  the 
interfacial  tension  of  the  pair,  and  <j)  is  the  contact  angle.  Thus  in  a  fixed  porous  medium,  the 
capillary  pressure-saturation  (P-S)  relationship  for  a  fluid  pair  will  scale  with  the  interfacial 
tension  between  the  fluid  pair,  as  long  as  the  contact  angles  are  similar.  In  particular,  the  two- 


3 


1  phase  water-oil  P-S  relationship  can  be  related  to  the  two-phase  water-air  P-S  relationship 

2  through, 

3 

^  Pw/o(®vv)~  ^P\v/a(®w)-  (3) 

*^w/a 

5 

6  In  general,  the  higher  the  capillary  pressure,  the  less  saturation  of  the  wetting  phase.  Also  the 

7  relation  is  always  hysteritic  depending  if  the  wetting  fluid  is  entering  (wetting)  or  leaving  (drying) 

8  a  particular  portion  of  the  porous  media. 

9 

10  From  Parlange  and  Hill  (1976)  and  Chandler  et  al.  (1997)  the  scaling  of  finger  width  into  the 

1 1  water  and  oil  system  is  then  given  by, 

12 


13 


14 


*^oil 


*air 


^vv/o 

^w/a 


Po  u(ef-9i) 


1- 


K 


w 


(4) 


15  where  doii  is  the  finger  width  for  the  water  oil  system,  dair  is  the  finger  width  for  the  water-air 

16  system,  and  K*  =  kg/p.  is  the  water  conductivity  behind  the  wetting  front.  This  expression  has 

17  been  derived  in  the  low  velocity  limit,  where  the  effects  of  viscosity  of  the  displaced  fluid  are 

18  assumed  to  be  negligible. 
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We  also  wish  to  understand  the  fluid  state  within  the  preferential  flow  path.  For  water  infiltrating 
an  air  dry  soil,  Selker  (1992b)  provided  a  description  of  the  water  saturation  behind  a  moving 
gravity  driven  finger.  We  follow  their  arguments  for  the  case  of  water  infiltrating  an  oil-saturated 
soil. 


Starting  with  the  one-dimensional  Richard’s  equation 


dt 


dz 


dz 


(5) 


9 

10  where  z  is  the  vertical  distance  (positive  downward),  K  is  the  unsaturated  conductivity  and  \|/w 

11  and  vj/o  are  the  total  water  and  oil  potentials  in  units  of  pressure  head  (cm  of  water),  respectively. 

12  This  assumes  that  when  preferential  flow  occurs  that  the  oil  moves  under  much  smaller  pressure 

13  gradients  than  the  water,  as  the  oil  can  move  through  the  whole  matrix  and  the  water  is  confined 

14  in  the  flow  path.  This  has  been  observed  to  be  true  for  water  and  soltrol,  but  must  break  down 

15  for  oils  of  high  viscosity. 

16 

17  For  water  infiltrating  a  less  dense  oil  the  potential  difference  is  the  sum  of  the  gravitational  and 

18  capillary  potentials, 

19 


20 


Vw  -  Vo  = 


Pw  Pq 
I  Pw 


z  +  Pw/o(0w)- 


(6) 


21 


7 


1  It  is  experimentally  observed  that  the  saturation  profile  within  the  finger  translates  downward  at  a 

2  constant  velocity  v,  and  thus  the  water  saturation  can  be  expressed  as  0(z-vt).  Defining  Ti=z-vt, 

3  substituting  Eq.  (6)  into  Eq.  (5)  and  integrating  converts  the  partial  differential  equation  into 

4 


5 


ve^  =K(e^) 


1- 


Po 

Pw 


^w/o^ 

dq  J 


(7) 


6 

7  which  relates  the  water  saturation,  the  finger  velocity,  the  unsaturated  conductivity,  and  the 

8  capillary  pressure  inside  the  finger.  This  expression  can  be  used  to  calculate  the  extent  of  the 

9  finger  tip,  and  the  unsaturated  conductivity  among  other  things. 

10 

11  At  a  particular  time  t,  the  capillary  pressure  within  the  finger  path  is  given  by 

12 

13  ^^w/o  _  j  _  Po  _  /g\ 

dz  K(0^) 

14 

15  Since  a  column  of  water  at  rest  has  dP/dz  =  1,  the  above  equation  predicts  that  the  water 

16  saturation  within  the  finger  will  appear  as  a  hanging  water  column  scaled  by  the  right  hand  side  of 

17  Eq  (8).  For  well-sorted  soils,  the  lower  portion  of  a  water  column  (or  equivalently  the  pressure  - 

18  saturation  relation)  has  a  region  of  constant  saturation  between  the  water  entry  pressure  and  the 

19  non-wetting  fluid  entry  pressure.  This  region  of  constant  has  been  observed  in  fingers  in  air-dry 

20  soil.  Since  the  saturation  is  constant,  by  integrating  Eq.  (8)  we  obtain  a  prediction  for  the  length 

21  of  the  tip  region  for  water  fingers  through  oil-saturated  sand. 
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Pw/o(oe)-Pw/o(we.) 


1- 


vG 


w 


Pw 


where  Ztip  is  the  expected  length  of  the  finger  tip,  and  Pw/o(w.e)  and  Pw/o(o.e.)  are  the  water-entry 
and  oil-entry  pressures,  respectively.  Assuming  the  pressure-saturation  scaling  in  Eq.  (3),  Ztip  can 
be  calculated  fi’om  the  known  water  and  air  entry  values  in  the  water  and  air  system,  the  densities, 
the  finger  velocity  and  saturation,  and  the  conductivity  of  the  soil. 

If  the  fluid  pressures  are  also  measured,  it  is  possible  to  obtain  measurements  for  the  unsaturated 
conductivity  using  Eq.  (7). 


Materials  and  Methods 

Experiments  to  observe  and  measure  finger  properties  were  performed  in  two  different  two- 
dimensional  chambers  uniformly  packed  with  sand.  The  sample  chambers  both  had  0.94  cm  thick 
polycarbonate  walls,  interior  thickness’  of  0.94  cm,  and  heights  of  55  cm,  but  different  widths  of 
30  cm  and  51  cm.  Attached  to  the  bottom  end  of  each  chamber  was  a  manifold  of  5  fluid  ports, 
and  the  top  end  was  left  open  to  the  atmosphere.  Chambers  were  packed  as  follows.  The  bottom 
of  the  chamber  was  filled  to  a  height  of  17  cm  with  Soltrol  220  (Phillips  Company),  a  non-toxic 
light  NAPL.  Clean,  dry,  sieved,  quartz  sand  (Unimin  Corporation)  was  dropped  through  the  top 
using  a  dual  randomizing  funnel  in  a  single  pour  to  nunimize  layering.  During  packing  and  for  one 


9 


1  minute  afterwards,  the  chamber  was  vibrated  to  enhance  settling.  Grain  sizes  of  12/20  and  20/30 

2  sand  were  used,  where  the  numbers  correspond  to  the  sieve  sizes  the  sand  passes  through  and 

3  does  not  pass  through,  respectively.  This  produced  a  tight  oil-saturated  pack  with  minimal 

4  entrapped  air  and  a  porosity  of  37%.  Water  fingers  in  these  sand  packs  had  uniform  widths  and 

5  velocities  throughout  the  chamber. 

6 

7  The  aqueous  solutions  used  in  the  experiments  consisted  of  either  water  or  water  with  a  dissolved 

8  surfactant.  Four  commonly  available  surfactants  were  used,  Alfonic  810-4.5  Ethoxylate 

9  (CONDEA  Vista  Company)  and  Surfynol  485  (Air  Products  and  Chemical,  Inc.)  are  ionic 

10  surfactants,  and  Neodol  (R)  25-7  (Shell  Chemical  Company)  and  Poly  Sodium  Vinyl  Sulfonate 

11  (Air  Products  and  Chemical,  Inc.)  are  non-ionic.  The  names  have  been  abbreviated  throughout  the 

12  text  as  alfonic,  surfynol,  neodol,  and  sulfonate,  respectively.  The  surfactants  were  mixed  with 

13  water  to  a  1  %  by  weight  solution.  Interfacial  tensions  between  each  solution  and  soltrol  220 

14  were  measured  using  a  ring  (du  Nouy  ring  method)  tensiomat.  For  all  the  aqueous  solutions, 

15  0.005  %  of  FD&C  blue  #  1  were  added  to  aid  in  the  visual  inspection  of  water  positions. 

16 

17  Infiltrations  of  the  aqueous  solution  were  performed  by  applying  the  aqueous  solution  to  the  top 

18  of  the  chamber  either  uniformly  and  through  a  point  source.  Uniform  applications  were 

19  performed  in  the  wider  chamber  (width  =  51  cm).  The  application  system  consisted  of  a  peristaltic 

20  pump  driven  point  source,  which  a  rotating  cam  shaft  moved  continuously  across  the  soil  surface. 

21  Uniform  application  rates  were  2.25,  6.0,  and  10.0  cm/hour.  Oil  was  removed  from  the  bottom  of 

22  the  chamber  through  a  constant  head  overflow  system  at  the  height  of  the  soil  surface.  For 
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1  uniform  applications,  the  chamber  was  placed  in  front  of  a  light  source  and  the  infiltrations  were 

2  observed  visually. 

3 

4  Point  source  applications  were  performed  in  the  smaller  chamber  (width  =  30  cm).  In  this  case, 

5  the  water  application  system  was  set  up  to  keep  the  fluid  content  in  the  chamber  constant.  This 

6  consisted  of  a  fixed  peristaltic  pump  which  applied  the  solution  through  a  tube  and  needle  at  2 

7  ml/min.  The  other  end  of  the  pump  tubing  was  connected  to  the  bottom  of  a  stoppered  250  mL 

8  Erlenmeyer  flask  containing  the  aqueous  solution,  which  was  in  turn  connected  to  the  bottom  of 

9  the  chamber,  thus  pulling  an  equal  amount  of  oil  out  of  the  chamber.  Using  this  system  to 

10  withdrawal  oil  from  the  bottom  of  the  chamber  as  water  was  being  applied  to  the  top  maintained  a 

11  constant  fluid  level  and  therefore  kept  the  oil  pressure  constant. 

12 

13  For  the  point  source  infiltrations,  fluid  contents  were  measured  using  synchrotron  x-ray 

14  attenuation.  The  dual  energy  attenuations  were  obtained  at  the  F-2  beam  line  of  the  Cornell  High 

15  Energy  Synchrotron  Source  using  a  set-up  described  in  detail  elsewhere  (DiCarlo  et  al.  1997). 

16  The  chamber  was  mounted  on  a  movable  platform  which  had  a  vertical  and  horizontal  range  of  50 

17  and  25  cm,  respectively,  and  a  position  repeatability  of  0.01  cm.  In  this  experiment  we  used  x-ray 

18  energies  of  25  and  50  KeV.  As  the  system  consisted  of  two  fluids,  only  the  25  KeV  x-ray 

19  attenuation  was  used  to  obtain  the  water  saturation,  with  the  50  KeV  x-ray  was  used  to  calibrate 

20  out  any  possible  soil  motions.  For  two  fluid  measurements  the  water  saturation  is  related  to  the 

21  attenuation  through. 
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2 

3  where  A  is  the  measured  attenuation,  Ao  is  the  attenuation  through  the  oil  saturated  soil,  x  is  the 

4  thickness  of  the  chamber,  and  |.i„  and  |.io  are  the  measured  25  KeV  attenuation  constants  for  the 

5  water  and  oil  phases,  respectively. 

6 

7  Fluid  attenuation  constants  were  obtained  using  a  polycarbonate  chamber  with  5  sequential  cells 

8  of  thickness  0.56  cm.  (Lenhard,  et  al.,  1988)  The  attenuation  was  measured  for  the  empty 

9  chamber  and  after  each  cell  was  ftlled  with  either  fluid.  The  measured  attenuations  were  linear 

10  with  fluid  thickness,  and  the  attenuation  constant  was  given  by  the  slope. 

11 

12  Water  and  oil  pressures  were  obtained  using  4  miniature  tensiometers  (Selker  et  al.,  1992a)  placed 

13  in  a  line  40  cm  above  the  bottom  of  the  chamber.  The  water  tensiometers  were  2  cm  left  of  the 

14  horizontal  center  and  on  the  center  (on  the  front  of  the  chamber) ,  while  the  oil  tensiometers  were 

15  2  cm  right  of  the  horizontal  center  and  on  the  center  (on  the  rear  of  the  chamber).  Each 

16  tensiometer  consisted  of  a  hollow  brass  housing,  with  one  end  connected  to  the  soil  through  a  0.8 

17  cm  diameter  porous  plate,  and  the  other  end  connected  to  a  pressure  transducer  through  a  rigid 

18  plastic  tube  and  a  three-way  valve.  The  other  port  of  each  three-way  valve  was  connected  to  a 

19  vertical  reservoir  of  water.  Prior  to  each  experiment  the  valve  was  opened  the  reservoir,  and  the 

20  transducers  were  calibrated  using  different  water  heights  in  the  reservoir.  Water  tensiometers 

21  were  filled  with  de-gassed  water,  and  had  a  sintered  glass  porous  plate  (Ace  Glass)  of  maximum 


22  pore  size  70-100  p,m.  Oil  tensiometers  were  filled  with  soltrol  220  and  had  a  sintered  stainless 
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1  Steel  porous  plate  (Mott  metallurgical)  of  pore  size  20  pm.  As  water  preferentially  wets  the  glass 

2  plate,  and  oil  preferentially  wets  the  stainless  steel  plate,  each  tensiometer  provided  a  connection 

3  to  their  particular  fluid  within  the  sand  pack.  The  capillary  pressures  were  obtained  by  subtracting 

4  the  water  pressure  from  the  oil  pressure. 

5 

6  Finger  widths  were  found  by  analyzing  tracings  from  the  experiments  and  from  horizontal 

7  attenuation  scans  through  the  finger.  Tracings  were  done  on  1  cm  graph  paper  taped  on  the  back 

8  of  the  chamber.  While  the  finger  was  infiltrating,  the  finger  was  then  traced  at  one  minute 

9  intervals.  The  width  was  then  measured  from  the  paper  tracings  every  5  cm  vertical.  Horizontal 

10  attenuation  scans  were  performed  after  the  finger  had  reached  the  bottom  of  the  chamber,  and 

1 1  widths  of  the  fingers  were  taken  to  be  the  full  width  at  half  maximum  of  the  water  saturation  in 

12  the  finger.  Both  procedures  produced  similar  results  and  the  reported  widths  are  an  average  of 

13  the  two  procedures. 

14 

15  Results 

16  Uniform  Water  Application 

17 

18  Due  to  the  difficulties  in  applying  water  uniformly  at  the  synchrotron  and  lack  of  control  of  the 

19  finger  positions,  only  qualitative  results  were  obtained  for  the  uniform  applications.  When  the 

20  aqueous  solution  was  applied  uniformly  to  the  top  of  the  oil  saturated  12/20  sand  pack,  for  every 

21  solution  the  water  moved  through  the  soil  pack  in  preferential  paths.  The  experiments  were 
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1  repeated  for  three  infiltration  rates  ranging  from  2.2  to  13  cm/hour.  The  widths  of  the  fingers  and 

2  the  pattern  of  the  fingers  varied  for  each  surfactant  used. 

3 

4  For  pure  water,  the  observed  preferential  flow  moved  through  the  chamber  at  constant  velocity 

5  and  with  constant  widths  between  5  and  8  cm.  The  flow  paths  are  sketched  in  Figure  la.  No 

6  splitting  or  merging  was  observed.  Higher  infiltration  rates  just  produced  more  fingers; 

7  infiltrations  of  5.6,  9.4,  and  13  cm/hour  produced  1,  2,  and  3  fingers  respectively. 

8 

9  For  the  surfactant  solutions,  sulfonate  produced  only  2  fingers  for  an  infiltration  rate  of  2.2 

10  cm/hour  and  produced  5  fingers  for  an  infiltration  rate  of  10  cm/hour.  Similarly  to  the  water 

11  solution,  these  fingers  did  not  merge  and  remained  the  same  width  and  velocity  throughout  the 

12  chamber.  The  surfynol  solution  initially  produced  10  fingers  for  all  flow  rates  with  widths  of 

13  roughly  2  cm.  Only  a  few  fingers  merged  and  6-8  fingers  remained  at  a  depth  of  50  cm  with 

14  roughly  the  same  2  cm  width.  The  low  interfacial  tension  solutions  containing  alfonic  and  neodol 

15  were  much  different.  Each  initially  produced  20  fingers  of  widths  1  cm  or  less  for  all  flow  rates. 

16  Many  of  these  fingers  quickly  merged  and  produced  fewer  and  faster  moving  fingers.  Typically, 

17  the  finger  merging  occurred  twice  for  each  infiltration,  a  primary  merge  at  a  depth  of  10-20  cm, 

18  and  a  secondary  merge  at  a  depth  of 40-50  cm.  Fig  lb  shows  the  sketch  of  the  flow  for  the  neodol 

19  solution. 
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1  Point  Source  Applications 

2 

3  All  point  source  applications  of  the  aqueous  fluids  produced  preferential  flow  paths.  Water, 

4  sulfonate,  and  surfynol  each  produced  one  single  finger  with  uniform  width  and  visually  similar  to 

5  the  fingers  observed  for  the  uniform  applications.  Neodol  and  alfonic,  however,  exhibited 

6  different  behavior  from  that  observed  in  the  uniform  applications.  Neodol  began  as  one  main 

7  finger,  but  small  fingers  quickly  began  to  splinter  off  the  main  finger.  These  other  fingers 

8  remained  small  and  did  not  grow  continuously  during  the  experiment.  Alfonic,  similar  to  neodol, 

9  had  many  fingers  which  splintered  off  the  main  finger,  but  the  fingers  continued  to  grow  and 

10  became  the  same  size  as  the  main  finger,  if  not  larger.  The  entire  chamber  became  filled  with 

11  alfonic  fingers,  and  the  entire  finger  pattern  closely  resembled  a  large  net.  The  fingers  were 

12  extremely  diffuse  and  pervasive,  especially  when  compared  to  the  one  water  finger.  Table  1  lists 

13  the  widths  and  velocities  of  the  main  fingers.  In  general,  water  fingers  were  the  widest  followed 

14  by  sulfonate,  surfynol,  neodol,  and  alfonic. 

15 

16  Pressure  and  saturation  measurements  were  taken  within  the  finger  fcr  each  solution.  Figure  2 

17  shows  the  water  saturation  and  capillary  pressure  between  the  water  and  oil  phases  in  the  initial 

18  finger  versus  time  when  the  invading  aqueous  solution  is  distilled  water.  Both  data  points  were 

19  taken  at  z  =  40  cm  (40  cm  above  the  bottom  of  the  chamber).  The  behavior  is  qualitatively  similar 

20  to  that  seen  for  water  preferentially  flowing  through  air-dry  sand.  The  finger  has  a  tip  with  a 

21  relatively  constant  water  saturation  of  0.30  cm7cm',  behind  which  the  water  saturation  quickly 
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1  drops  to  0. 17  cm^/cm^  where  it  levels  off  The  capillary  pressure  of  the  is  lowest  at  the  finger  tip 

2  (0.5  cm)  and  rises  continuously  after  the  tip  before  leveling  off  at  4  cm. 

3 

4  Figure  3  shows  the  water  saturation  and  capillary  pressure  within  the  finger  when  the  invading 

5  aqueous  solution  contains  1  %  by  weight  surfynol.  The  same  qualitative  features  as  seen  in  the 

6  water  finger  are  present:  highly  saturated  tip,  and  saturations  decreasing  and  pressures  increasing 

7  behind  the  tip.  But,  an  additional  feature  produced  by  surfactants  is  a  slow  increase  in  the  water 

8  saturation  of  roughly  0.03  cm^/cm^  well  behind  the  finger  tip.  This  feature  is  associated  with  an 

9  decrease  in  the  capillary  pressure  after  the  initial  increase.  Other  than  the  slow  increase  at  the 

10  end,  the  saturations  in  the  finger  are  quantitatively  similar  for  the  surfactant  fingers  with  a  tip 

11  saturation  of  approximately)  0.30  cmVcm^  and  saturation  directly  behind  the  tip  of  0.17  cmVcm^. 

12  The  surfactant  fingers  are  quantitatively  different  in  the  actual  length  of  the  nearly  saturated 

13  region,  and  the  pressure  rise  behind  the  tip. 

14 

15  Identical  sets  of  water  saturation  and  pressure  data  were  taken  for  the  other  surfactant  solutions 

16  and  for  a  coarser  sand  pack  of  12/20  sand.  Table  1  lists  the  finger  widths,  finger  velocities,  finger 

17  pressure  rises,  and  the  lengths  of  the  saturated  tips  for  all  of  the  data  sets.  The  finger  properties 

18  show  definite  trends  with  the  interfacial  tension  of  the  fluids.  As  the  interfacial  tension  between 

19  the  solutions  is  lowered  by  the  surfactants,  the  pressure  rises  behind  the  wetting  front  decrease, 

20  the  finger  widths  decrease,  and  the  lengths  of  the  saturated  tips  decrease. 
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1  Discussion 

2 

3  Due  to  the  qualitative  similarity  between  the  fingers  observed  here  and  those  observed  in  water 

4  infiltrating  air-dry  sand,  we  will  analyze  the  data  in  the  framework  developed  for  the  latter  case.  If 

5  this  formulation  is  correct,  the  finger  properties  such  as  the  finger  widths,  the  pressure  drop 

6  behind  the  finger  and  the  length  of  the  finger  tips  should  be  scaled  through  the  above  equations. 

7 

8  Figure  5  shows  the  finger  width  versus  the  interfacial  tension  for  the  12/20  and  20/30  sand.  The 

9  lines  are  the  predicted  scaling  relation  in  Eq.  4  with  the  slope  determined  by  the  widths  of  the 

10  water  fingers  through  air-dry  sand  (2.0  cm  for  12/20  sand  and  2.5  cm  for  20/30  sand).  The  trend 

1 1  in  the  observed  data  follows  the  theoretical  predictions;  with  the  smallest  fingers  occuring  with 

12  solutions  of  alfonic  and  neodol,  which  have  the  smallest  interfacial  tension  values.  At  these  low 

13  interfacial  tensions,  is  also  where  the  predictions  differ  the  most,  with  the  fingers  being  wider  than 

14  predicted.  This  is  not  surprising  as  we  have  rarely  observed  fingers  of  widths  smaller  than  1.5  cm, 

15  and  never  smaller  than  1  cm.  Thus  it  appears  that  the  scaling  is  roughly  true,  except  that  there 

16  appears  to  be  some  minimum  finger  width.  Clearly  the  theory  must  break  down  at  the  pore  scale 

17  and  thus  must  have  a  length  scale  smaller  than  which  it  is  not  appropriate. 

18 

19  The  simultaneous  measurements  of  pressure  and  saturation  allow  us  to  quickly  obtain  the 

20  pressure-saturation  relationship  for  this  water  oil  system.  Figure  5  shows  the  measured 

21  relationship  from  the  water  infiltration  and  the  re-infiltration  into  12/20  sand.  As  the  wetting 

22  process  is  very  abrupt,  both  curves  are  on  drying  branches  of  the  P-S  relationship. 
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1 

2  Once  the  fingers  are  formed,  hysteresis  in  the  above  pressure-saturation  relation  plays  a  large  part 

3  in  the  stability  of  the  fingers.  The  water  does  not  enter  the  oil  saturated  sand  until  the  capillary 

4  pressure  decreases  to  0.5  cm,  at  which  the  water  enters  at  a  saturation  of  0.30  cmVcm''  (point  A). 

5  Behind  the  finger  tip  the  capillary  pressure  increases  and  the  water  saturation  decreases  (point  B). 

6  Since  the  capillary  pressure  behind  the  finger  tip  quickly  becomes  much  greater  than  the  water 

7  entry  pressure,  the  water  is  confined  within  the  finger  and  cannot  spread.  The  finger  remains  at 

8  point  B  until  the  water  is  turned  off  after  which  the  capillary  pressure  increases  and  water 

9  saturation  decreases  further  (point  C).  On  a  reinfiltration  the  old  finger  core  wets  quickly  to  point 

10  D,  where  it  remains  until  the  water  is  turned  off  after  which  it  dries  to  point  E.  Again,  on 

1 1  reinfiltration  the  capillary  pressure  remains  greater  than  the  water  entry  pressure  and  the  water 

12  finger  does  not  expand. 

13 

14  These  pressure  drops  differ  for  each  surfactant  used.  Figure  6  shows  the  maximum  pressure  drop 

15  inside  the  finger  (Point  A  to  Point  B)  versus  the  interfacial  tension  of  the  fluids.  The  line  shown  is 

16  the  predicted  dependence  on  interfacial  tension.  The  observed  data  correlates  well  over  the  variety 

17  of  surfactants  used  and  for  the  water  into  air  dry  sand  case. 

18 

19  The  pressure  data  can  give  us  insight  into  the  finger  patterns  for  different  fluid  pairs.  For  viscous 

20  fingering,  it  has  been  known  for  some  time  that  at  low  interfacial  tensions,  the  finger  pattern 

21  becomes  very  unstable  or  “chaotic”  with  large  amounts  of  finger  splitting.  (Bensimon  et  al., 

22  1986).  For  gravity  driven  fingers,  we  also  observe  finger  splitting  at  the  lowest  interfacial 
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1  tensions.  This  splitting  can  be  understood  in  terms  of  the  water  entry  pressure  of  the  soil  and  the 

2  pressure  rise  inside  the  finger  tips.  For  large  Interfacial  tensions,  the  pressure  rise  within  the  finger 

3  is  large,  and  the  pressure  behind  the  tip  of  the  finger  becomes  greater  than  the  water-entry 
pressure.  Therefore  there  will  not  be  enough  pressure  for  the  finger  to  expand,  and  the  width  will 

5  not  increase.  This  is  what  is  observed  for  infiltrations  of  suriynol,  sulfonate,  and  water.  For  small 
interfacial  tensions  (  e  g.  neodol  and  alfontc),  however,  the  pressure  rise  within  the  finger  is  small, 

7  and  the  pressure  behind  the  tip  of  the  finger  remains  close  to  or  at  the  water-entry  pressure. 

8  Therefore  there  will  be  enough  pressure  for  the  finger  to  expand,  and  splitting  and  secondary 

9  fingers  can  occur.  This  is  what  is  observed  for  infiltrations  of  neodol  and  alfonic, 

10 

1 1  What  ts  the  role  of  viscosity  in  these  gravity  driven  fingers?  In  measurements  of  the  oil  pressure 

12  directly  outside  the  finger  path,  we  observe  no  change  in  the  oil  pressure  as  the  finger  passes 

13  (Chandler  et  al.  1997).  Thus  for  soltrol  and  water  with  a  viscosity  ratio  of  pVti.  =  4  ,  the  major 

14  pressure  changes  are  in  the  water  phase  and  the  assumption  that  the  oil  acts  inviscidly  (oil 

1 5  movement  with  little  or  no  pressure  gradients)  during  preferential  flow  appears  to  be  valid. 

14  Although  there  are  no  observed  effects  of  oil  viscosity  in  the  bulk  oil  phase,  other  effects  can  be 

17  attnbuted  to  oil  viscosity.  The  most  prominent  of  these  viscous  effects  is  that  the  finger  is  not 

18  completely  saturated  with  water,  unlike  fingers  in  the  water  and  air  system.  This  is  most  likely  due 

19  to  the  oil  not  being  able  to  move  out  of  all  the  pores  fast  enough  during  the  infiltration  because  of 

20  its  viscosity.  The  lower  saturation  in  the  tip  also  affects  the  behavior  of  the  finger  behind  the  tip, 

21  as  the  lower  tip  saturation  yields  a  water  conductivity  of  roughly  half  of  rhe  saturated 

22  conductivity. 
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1 

2  Behind  the  finger  tip,  the  water  saturation  drops  and  oil  re-enters  the  finger  path.  Viscosity  may 

3  affect  the  time  scale  over  which  this  proceeds,  and  thus  the  tip  length  of  the  finger  tip.  Figure  7 

4  plots  the  measured  length  of  the  finger  tip  versus  that  predicted  in  equation  (9).  The  predictions 

5  were  made  by  assuming  that  the  pressure  difference  between  the  oil  entry  and  water  entry  values 

6  scales  with  interfacial  tension,  the  water  saturation  was  0.30  cm'Vcm",  and  the  conductivity  at  this 

7  saturation  was  measured  separately  to  be  0.45  cm/sec  for  12/20  sand  and  0. 15  cm/sec  for  20/30 

8  sand.  The  trend  is  followed  well  by  the  data,  but  the  observed  finger  widths  are  generally  higher 

9  than  the  predicted  widths.  This  is  what  would  be  expected  if  the  displaced  fluid  has  an  appreciable 

10  viscosity,  as  the  viscosity  will  limit  the  rate  at  which  oil  can  return  into  the  finger  path  thus  making 

11  the  finger  longer.  We  have  only  used  an  oil  of  one  viscosity,  it  would  be  interesting  to  see  how  the 

12  tip  length 

13 

14  Summary 

15 

16  In  summary,  when  water  is  gravitionally  driven  through  an  oil-saturated  sand  pack,  fingers  occur 

17  which  resemble  those  seen  when  water  infiltrates  a  dry  sand  pack.  As  opposed  to  viscous  fingers 

18  in  Hele-Shaw  cells,  behind  the  fingers  tips  oil  re-enters  the  finger  paths,  and  the  capillary  pressure 

19  -  saturation  relationship  controls  the  spread  of  the  fingers.  These  phenomena  are  described  well 

20  by  the  theory  developed  for  fingers  from  water  infiltration  with  important  scaling  by  the  density 

21  differences,  and  the  interfacial  tensions,  with  the  viscosity  of  the  oil  playing  a  mino?’  role. 


22 


20 


1  Possible  future  work  involves  looking  into  observing  the  crossover  from  gravity  to  viscosity 

2  driven  fingers  in  porous  media.  Even  with  the  oil  being  more  4  times  more  viscous  than  the 

3  water,  the  pressure  drops  occur  inside  the  water  finger.  Future  work  consists  of  answering 

4  questions  such  as;  Is  the  re-entry  of  the  fluid  behind  the  finger  a  signature  of  gravity  driven 

5  fingers?  Is  the  lateral  expansion  of  fingers  depend  on  if  viscous  or  gravitational  forces  are 

6  dominating?  And  what  is  the  role  of  pressure  saturation  relation  in  viscous  fingering,  as  it  is 

7  usually  assumed  that  the  interface  is  abrupt? 

8 
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1  Table  1: 

2 


12/20  sand 
Fluid 

Interfacial 

Finger 

Tension 

Width 

(dynes/cm) 

(cm) 

Water 

24 

4.8 

Sulfonate 

21 

2.8 

Surfynol 

9 

2.4 

Neodol 

4 

2.2 

Alfonic 

3 

1.6 

20/30  sand 
Fluid 

Interfacial 

Finger 

Tension 

Width 

(dynes/cm) 

(cm) 

Water 

24 

5.1 

Sulfonate 

21 

4.2 

Surfynol 

9 

5.0 

Neodol 

4 

2.1 

Alfonic 

3 

1.9 

3 

4 


Finger 

Pressure 

Finger 

Tip  Length 

Drop 

Velocity 

(cm) 

(cm  of  H2O) 

(cm/sec) 

10.2 

4.0 

0.03 

8.0 

2.6 

0.055 

8.0 

1.5 

0.06 

4.0 

<0.2 

0.09 

3.5 

<0.2 

0.09 

Finger 

Pressure 

Finger 

Tip  Length 

Drop 

Velocity 

(cm) 

(cm  of  H2O) 

(cm/sec) 

25 

3.5 

0.022 

28 

3.5 

0.025 

13 

1.8 

0.03 

6 

0.5 

0.043 

6 

<0.5 

0.045 

5  Summary  of  the  internal  finger  properties  for  each  infiltrating  fluid  and  for  both  sand  sieve  sizes. 
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1  Figure  Captions: 

2 

3  Fig.  1.  Sketch  of  fingering  patterns  when  a  aqueous  solution  is  infiltrated  uniformly  on  an  oil- 

4  saturated  sand,  a)  water,  b)  water  with  neodol  as  a  surfactant. 

5 

6  Fig.  2.  Water  saturation  and  capillary  pressure  inside  a  water  finger  infiltrating  oil  saturated  12/20 

7  sand.  The  highest  saturation  and  capillary  pressures  are  at  the  tip  of  the  finger, 

s 

9  Fig  3.  Water  saturation  and  capillary  pressure  inside  a  surfynol  solution  finger  infiltrating  oil 

10  saturated  12/20  sand.  The  pressure  changes  behind  the  tip  are  smaller  than  for  the  water  case, 

11  and  the  saturation  of  the  surfynol  solution  rises  slightly  after  the  finger  tip. 

12 

13  Fig.  4.  Capillary  pressure  -  saturation  relationship  obtained  from  the  water  finger  into  oil- 

14  saturated  12/20  sand.  The  drying  from  A  to  B  takes  place  in  the  initial  finger,  drying  to  point  C 

15  occurs  when  the  infiltration  is  stopped.  Point  D  is  reached  on  re-infiltration,  and  point  E  after  the 

16  re-infiltration  is  stopped. 

17 

18  Fig.  5.  Finger  widths  versus  interfacial  tension  for  the  five  solutions  and  two  sands  used.  The 

19  solid  and  dotted  lines  are  the  expected  scaling  from  fingers  in  air-dry  12/20  and  20/30  sand, 

20  respectively. 

21 
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1  Fig.  6.  Pressure  drops  behind  the  finger  tip  versus  interfacial  tension  for  the  five  solutions  and 

2  two  sands  used.  The  solid  and  dotted  lines  are  the  expected  scaling  from  fingers  in  air-dry  12/20 

3  and  20/30  sand,  respectively. 

4 

5  Fig.  7.  Measured  finger  tip  length  versus  the  predicted  tip  length  from  Eq.  8.  Tip  lengths  are 

6  slightly  higher  than  predicted  by  the  1:1  solid  line. 
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